Search

In a previous note [Tasso, Phys. Plasmas 2, 1789 (1995)], dissipative magnetized Trkal flows were proved to be Lyapunov stable nonlinearly under the conditions of sufficiently small Reynolds number and of an upper bound of the "force-free" current. An alternative Lyapunov functional containing the cross helicity and introduced recently for conservative magnetized Beltrami flows [Yoshida , J. Math. Phys. 44, 2168 (2003)] is extended here to the Trkal flows mentioned above. Stability is now proven under limitations on maximal velocity and current without introducing an upper bound on the Reynolds number. (C) 2003 American Institute of Physics. Physics of Plasmas

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.

It is found that the ideal magnetohydrodynamic equilibrium of an axisymmetric gravitating magnetically confined plasma with incompressible flows is governed by a second-order elliptic differential equation for the poloidal magnetic flux function containing five flux functions coupled with a Poisson equation for the gravitation potential, and an algebraic relation for the pressure. This set of equations is amenable to analytic solutions. As an application, the magnetic-dipole static axisymmetric equilibria with vanishing poloidal plasma currents derived recently by Krasheninnikov, Catto, and Hazeltine [Phys. Rev. Lett. {\bf 82}, 2689 (1999)] are extended to plasmas with finite poloidal currents, subject to gravitating forces from a massive body (a star or black hole) and inertial forces due to incompressible sheared flows. Explicit solutions are obtained in two regimes: (a) in the low-energy regime $\beta_0\approx \gamma_0\approx \delta_0 \approx\epsilon_0\ll 1$, where $\beta_0$, $\gamma_0$, $\delta_0$, and $\epsilon_0$ are related to the thermal, poloidal-current, flow and gravitating energies normalized to the poloidal-magnetic-field energy, respectively, and (b) in the high-energy regime $\beta_0\approx \gamma_0\approx \delta_0 \approx\epsilon_0\gg 1$. It turns out that in the high-energy regime all four forces, pressure-gradient, toroidal-magnetic-field, inertial, and gravitating contribute equally to the formation of magnetic surfaces very extended and localized about the symmetry plane such that the resulting equilibria resemble the accretion disks in astrophysics., Comment: 12 pages, latex, to be published in Geophys. Astrophys. Fluid Dynamics

From the equilibrium solution of a toroidal plasma with toroidal velocity it is possible to define a relaxed energy state. This state corresponds to the minimum total energy at a given current and under certain assumptions retains the current and pressure profiles of the static case. For large toroidal radius it produces unique current and pressure profiles which however are completely different to those of the static case.

It is shown that the magnetohydrodynamic equilibrium states of an axisymmetric toroidal plasma with finite resistivity and flows parallel to the magnetic field are governed by a second-order partial differential equation for the poloidal magnetic flux function $\psi$ coupled with a Bernoulli type equation for the plasma density (which are identical in form to the corresponding ideal MHD equilibrium equations) along with the relation $\Delta^\star \psi=V_c \sigma$. (Here, $\Delta^\star$ is the Grad-Schl\"{u}ter-Shafranov operator, $\sigma$ is the conductivity and $V_c$ is the constant toroidal-loop voltage divided by $2 \pi $). In particular, for incompressible flows the above mentioned partial differential equation becomes elliptic and decouples from the Bernoulli equation [H. Tasso and G. N. Throumoulopoulos, Phys. Plasmas {\bf 5}, 2378 (1998)]. For a conductivity of the form $\sigma=\sigma(R, \psi)$ ($R$ is the distance from the axis of symmetry) several classes of analytic equilibria with incompressible flows can be constructed having qualitatively plausible $\sigma$ profiles, i.e. profiles with $\sigma$ taking a maximum value close to the magnetic axis and a minimum value on the plasma surface. For $\sigma=\sigma(\psi)$ consideration of the relation $\Delta^\star\psi = V_c \sigma(\psi)$ in the vicinity of the magnetic axis leads therein to a proof of the non-existence of either compressible or incompressible equilibria. This result can be extended to the more general case of non-parallel flows lying within the magnetic surfaces., Comment: Post script file, 17 pages, no figures, to be published in J. Plasma Physics

The impact of reversed magnetic shear (s < 0) along with a radial electric field E-r on negative-energy perturbations (NEP's) is investigated for the case of cylindrical tokamak-like equilibria. NEP's can imply instabilities leading to anomalous transport [G. N. Throumoulopoulos and D. Pfirsch, Phys. Rev. E 56, 5979 (1997) and references therein]. For equilibria with E-r corresponding to \e(nu)phi\much less than beta T-nu and eta(nu)< 4/3 (phi is the electrostatic potential, beta equivalent to 8 pi P/B-2, eta(nu)equivalent to partial derivative lnT(nu)/partial derivative lnN(nu) and nu denotes the particle species) s < 0 results in a reduction of the fraction of particles associated with NEP's (active particles) as compared with equilibria with s > 0, by making the pressure profile steeper. The reduction is stronger the closer the minimum of the safety factor (q(min)) to the plasma center and the lower the negative value of s. For equilibria with \e(nu)phi\much less than beta T-nu and eta(nu)> 4/3 the fraction of active particles decreases as the magnetic tension B-theta(2)/r becomes larger. For E-r< 0 corresponding to \e(nu)phi\approximate to beta T-nu the electric field leads to a reduction of the fraction of active ions for eta(i)< 4/3 and of the fraction of active electrons for eta(e)> 4/3 in addition to that caused by s < 0 and B-theta(2)/r, respectively. For T(nu)greater than or equal to\e(nu)phi\much greater than beta T-nu, which corresponds to experimental observations, the reduction of active particles caused by E-r< 0 is insensitive to the sign and the value of s, the position of q(min), and the value of B-theta(2)/r. In addition, it is shown that the experimentally evidenced facts that (a) s < 0 is associated with a steepness of the pressure profile and (b) the position of q(min) is displaced towards the plasma edge as the toroidal current is increased are consistent with equilibrium considerations. It is emphasized that these results cannot yet give a complete picture, since the present paper still neglects toroidal effects. (C) 1999 American Institute of Physics. [S1070-664X(99)03208-5]. Physics of Plasmas

A recent study on axisymmetric ideal magnetohydrodynamic equilibria with incompressible flows [H. Tasso and G. N. Throumoulopoulos, Phys. Plasmas {\bf 5}, 2378 (1998)] is extended to the generic case of helically symmetric equilibria with incompressible flows. It is shown that the equilibrium states of the system under consideration are governed by an elliptic partial differential equation for the helical magnetic flux function $\psi$ containing five surface quantities along with a relation for the pressure. The above mentioned equation can be transformed to one possessing differential part identical in form to the corresponding static equilibrium equation, which is amenable to several classes of analytic solutions. In particular, equilibria with electric fields perpendicular to the magnetic surfaces and non-constant-Mach-number flows are constructed. Unlike the case in axisymmetric equilibria with isothermal magnetic surfaces, helically symmetric $T=T(\psi)$ equilibria are over-determined, i.e., in this case the equilibrium equations reduce to a set of eight ordinary differential equations with seven surface quantities. In addition, it is proved the non-existence of incompressible helically symmetric equilibria with (a) purely helical flows (b) non-parallel flows with isothermal magnetic surfaces and the magnetic field modulus being a surface quantity (omnigenous equilibria)., Comment: Latex file, 13 pages, accepted in J. Plasma Phys

It is shown that the ideal magnetohydrodynamic (MHD) equilibrium states of an axisymmetric plasma with incompressible flows are governed by an elliptic partial differential equation for the poloidal magnetic flux function psi containing five surface quantities along with a relation for the pressure. Exact equilibria are constructed including those with nonvanishing poloidal and toroidal flows and differentially varying radial electric fields. Unlike the case in cylindrical incompressible equilibria with isothermal magnetic surfaces which should have necessarily circular cross sections [G. N. Throumoulopoulos and H. Tasso, Phys. Plasmas 4, 1492 (1997)], no restriction appears on the shapes of the magnetic surfaces in the corresponding axisymmetric equilibria. The latter equilibria satisfy a set of six ordinary differential equations which for flows parallel to the magnetic field B can be solved semianalytically. In addition, it is proved the nonexistence of incompressible axisymmetric equilibria with (a) purely poloidal flows and (b) nonparallel flows with isothermal magnetic surfaces and \B\ = \B\(psi) (omnigenous equilibria). (C) 1998 American Institute of Physics. [S1070-664X(98)02806-7] Physics of Plasmas

The impact of an equilibrium radial electric field $E $ on negative-energy perturbations (NEPs) (which are potentially dangerous because they can lead to either linear or nonlinear explosive instabilities) in cylindrical equilibria of magnetically confined plasmas is investigated within the framework of Maxwell-drift kinetic theory. It turns out that for wave vectors with a non-vanishing component parallel to the magnetic field the conditions for the existence of NEPs in equilibria with E=0 [G. N. Throumoulopoulos and D. Pfirsch, Phys. Rev. E 53, 2767 (1996)] remain valid, while the condition for the existence of perpendicular NEPs, which are found to be the most important perturbations, is modified. For $|e_i\phi|\approx T_i$ ($\phi$ is the electrostatic potential) and $T_i/T_e > \beta_c\approx P/(B^2/8\pi)$ ($P$ is the total plasma pressure), a case which is of operational interest in magnetic confinement systems, the existence of perpendicular NEPs depends on $e_\nu E$, where $e_\nu$ is the charge of the particle species $\nu$. In this case the electric field can reduce the NEPs activity in the edge region of tokamaklike and stellaratorlike equilibria with identical parabolic pressure profiles, the reduction of electron NEPs being more pronounced than that of ion NEPs., Comment: 30 pages, latex

The equilibrium of a resistive axisymmetric plasma with purely toroidal flow surrounded by a conductor is investigated within the framework of the nonlinear magnetohydrodynamic theory. It is proved that a) the poloidal current density vanishes and b) apart from an idealized case the pressure profile should vanish on the plasma boundary. For the cases of isothermal magnetic surfaces, isentropic magnetic surfaces and magnetic surfaces with constant density, the equilibrium states obey to an elliptic partial differential equation for the poloidal magnetic flux function, which is identical in form to the corresponding equation governing ideal equilibria. The conductivity, which can not be either uniform or a surface quantity, results however in a restriction of the possible classes of equilibrium solutions, e.g., for the cases considered the only possible equilibria with Spitzer conductivity are of cylindrical shape., Comment: latex, 16 pages, 1 figure available on request to the author (gthroum@cc.uoi.gr), to be published in J. Plasma Physics

It is proved that (a) the solutions of the ideal magnetohydrodynamic equation, which describe the equlibrium states of a cylindrical plasma with purely poloidal flow and arbitrary cross sectional shape [G. N. Throumoulopoulos and G. Pantis, Plasma Phys. and Contr. Fusion 38, 1817 (1996)] are also valid for incompressible equlibrium flows with the axial velocity component being a free surface quantity and (b) for the case of isothermal incompressible equilibria the magnetic surfaces have necessarily circular cross section., Comment: 7 pages, latex, no figures

The conditions for the existence of negative-energy perturbations (which could be nonlinearly unstable and cause anomalous transport) are investigated in the framework of linearized collisionless Maxwell-drift kinetic theory for the case of equilibria of magnetically confined, circularly cylindrical plasmas and vanishing initial field perturbations. For wave vectors with a non-vanishing component parallel to the magnetic field, the plane equilibrium conditions (derived by Throumoulopoulos and Pfirsch [Phys Rev. E {\bf 49}, 3290 (1994)]) are shown to remain valid, while the condition for perpendicular perturbations (which are found to be the most important modes) is modified. Consequently, besides the tokamak equilibrium regime in which the existence of negative-energy perturbations is related to the threshold value of 2/3 of the quantity $\eta_\nu = \frac {\partial \ln T_\nu} {\partial \ln N_\nu}$, a new regime appears, not present in plane equilibria, in which negative-energy perturbations exist for {\em any} value of $\eta_\nu$. For various analytic cold-ion tokamak equilibria a substantial fraction of thermal electrons are associated with negative-energy perturbations (active particles). In particular, for linearly stable equilibria of a paramagnetic plasma with flat electron temperature profile ($\eta_e=0$), the entire velocity space is occupied by active electrons. The part of the velocity space occupied by active particles increases from the center to the plasma edge and is larger in a paramagnetic plasma than in a diamagnetic plasma with the same pressure profile. It is also shown that, unlike in plane equilibria, negative-energy perturbations exist in force-free reversed-field pinch equilibria with a substantial fraction of active particles., Comment: 31 pages, latex

The equilibrium of a cylindrical plasma with purely poloidal mass flow and cross section of arbitrary shape is investigated within the framework of the ideal MHD theory. For the system under consideration it is shown that only incompressible flows are possible and, conscequently, the general two dimensional flow equilibrium equations reduce to a single second-order quasilinear partial differential equation for the poloidal magnetic flux function $\psi$, in which four profile functionals of $\psi$ appear. Apart from a singularity occuring when the modulus of Mach number associated with the Alfv\'en velocity for the poloidal magnetic field is unity, this equation is always elliptic and permits the construction of several classes of analytic solutions. Specific exact equlibria for a plasma confined within a perfectly conducting circular cylindrical boundary and having i) a flat current density and ii) a peaked current density are obtained and studied., Comment: Accepted to Plasma Physics & Controlled Fusion, 14 pages, revtex

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

The cat-eyes steady state solution in the framework of hydrodynamics describing an infinite row of identical vortices is extended to the magnetohydrodynamic equilibrium equation with incompressible flow of arbitrary direction. The extended solution covers a variety of equilibria including one- and two-dimensional generalized force-free and Harris-sheet configurations which are preferable from those usually employed as initial states in reconnection studies. Although the vortex shape is not affected by the magnetic field, the flow in conjunction with the equilibrium nonlinearity has a strong impact on isobaric surfaces by forming pressure islands located within the cat-eyes vortices. More importantly, a magnetic-field-aligned flow of experimental fusion relevance and the flow shear have significant stabilizing effects in the region of pressure islands. The stable region is enhanced by an external axial ("toroidal") magnetic field., 17 pages, 8 figures

Several classes of analytic equilibria of a toroidal axisymmetric plasma with toroidal mass flow are obtained for isothermal as well as isentropic magnetic surfaces under various pressure and toroidal current density profiles. The effects of flow on certain global equilibrium and stability characteristics are examined. It is shown that for a wide family of equilibria the average beta 〈 β〉 does not depend explicitly on flow and falls within the limits of static equilibria. In particular, for the case of field‐reversed configurations 〈 β〉 is totally independent of plasma flow. Calculations based on a specific equilibrium solution for a tokamak with fixed average toroidal magnetic energy show that plasma elongation is an increasing function of the Mach number whereas the safety factor is a decreasing function of it.

We derive exact solutions of a linear form of the Grad-Shafranov (GS) equation, including incompressible equilibrium flow, using ansatz-based similarity reduction methods. The linearity of the equilibrium equation allows linear combinations of solutions in order to obtain axisymmetric MHD equilibria with closed and nested magnetic surfaces which are favorable for the effective confinement of laboratory plasmas. In addition, employing the same reduction methods we obtain analytical solutions for several non-linear forms of the GS equation. In this context analytic force-free solutions in both linear and nonlinear regimes are also derived., 5 pages, 3 figures

Under the assumption of time-dependent net core breeding gain of fusion-fission symbiotic systems, analytical forms for the stockpile fuel trajectories have been derived, both for continuous and discontinuous operation. Also, usable regions for characteristic system parameters such as doubling time, minimum inventory, minimum inventory time, and minimum initial inventory have been defined.It is shown that fuel trajectories defined previously are underestimating the importance of the time dependence of the net core breeding gain.

Vlasov equilibria of axisymmetric plasmas with vacuum toroidal magnetic field can be reduced, up to a selection of ions and electrons distributions functions, to a Grad-Shafranov-like equation. Quasineutrality narrow the choice of the distributions functions. In contrast to two-dimensional translationally symmetric equilibria whose electron distribution function consists of a displaced Maxwellian, the toroidal equilibria need deformed Maxwellians. In order to be able to carry through the calculations, this deformation is produced by means of either a Heaviside step function or an exponential function. The resulting Grad-Shafranov-like equations are established explicitly., 11 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.