BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

G. N. Throumoulopoulos

G. N. Throumoulopoulos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.plasm-phphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

19works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Neural network tokamak equilibria with incompressible flows

We present several numerical solutions to a generalized Grad-Shafranov equation (GGSE), which governs axisymmetric plasma equilibria with incompressible flows of arbitrary direction, using fully connected, feed-forward, deep neural networks, also known as multi-layer perceptrons. Such artificial neural network (ANNs) are trained to approximate tokamak-relevant equilibria upon minimizing the GGSE mean squared residual in the plasma volume and the poloidal flux function on the plasma boundary. Solutions for the Solovev and the general linearizing ansatz for the free functions involved in the GGSE are obtained and benchmarked against known analytic solutions. We also construct a non-linear equilibrium incorporating characteristics relevant to the high confinement mode. In our numerical experiments it was observed that changing the radial distribution of the training points has a surprisingly small effect on the accuracy of the trained solution. In particular it is shown that localizing the training points at the plasma edge results in ANN solutions that describe quite accurately the entire magnetic configuration, thus demonstrating the interpolation capabilities of the ANNs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Energy-Casimir, dynamically accessible, and Lagrangian stability of extended magnetohydrodynamic equilibria

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical equilibria with pressure anisotropy and incompressible plasma rotation parallel to the magnetic field

It is believed that plasma rotation can affect the transitions to the advanced confinement regimes in tokamaks. In addition, in order to achieve fusion temperatures modern tokamaks rely on auxiliary heating methods. These methods generate pressure anisotropy in the plasma. For incompressible rotation with pressure anisotropy the equilibrium is governed by a Generalized Grad-Shafranov (GGS) equation and a decoupled Bernoulli-type equation for the effective pressure, $\bar{p}=(p_\parallel+p_\perp)/2$, where $p_\parallel$ ($p_\perp$) is the pressure tensor element parallel (perpendicular) to the magnetic field. In the case of plasma rotation parallel to the magnetic field the GGS equation can be transformed to one equation identical in form with the GS equation. In this study by making use of the aforementioned property of the GGS equation for parallel plasma rotation we have constructed ITER-like numerical equilibria by extending HELENA, an equilibrium fixed-boundary solver and examined the impact of rotation and anisotropy on certain equilibrium quantities. The main conclusions are that the addition of pressure anisotropy to rotation allows the profile shaping of the equilibrium quantities in much more extent thus favouring the confinement and allows extension of the parametric space of the Mach number corresponding to higher values. Furthermore, the impact of pressure anisotropy in the equilibrium quantities is stronger than that of the rotation, for most of the quantities examined. For the pressure components the impact of the pressure anisotropy is the same regardless of whether the power is deposited parallel or perpendicular to the magnetic surfaces, thus implying that there is no preferable heating direction, while for the current density, the heating parallel to the magnetic surfaces seems to be beneficial for the current-gradient driven instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the linear stability of anisotropic pressure equilibria with field-aligned incompressible flow

We derive a sufficient condition for the linear stability of plasma equilibria with incompressible flow parallel to the magnetic field, $\bf B$, constant mass density and anisotropic pressure such that the quantity $σ_d= μ_0(P_\parallel -P_\perp)/B^2$, where $P_\parallel$ ($P_\perp$) is the pressure tensor element parallel (perpendicular) to $\bf B$, remains constant. This condition is applicable to any steady state without geometrical restriction. The condition, generalising the respective condition for MHD equilibria with isotropic pressure and constant density derived in Throumoulopoulos & Tasso (2007), involves physically interpretable terms related to the magnetic shear, the flow shear and the variation of total pressure perpendicular to the magnetic surfaces. On the basis of this condition we prove that if a given equilibrium is linearly stable, then the ones resulting from the application of Bogoyavlenskij symmetry transformations are linearly stable too, provided that a parameter involved in those transformations is positive. In addition, we examine the impact of pressure anisotropy, flow, and torsion of a helical magnetic axis, for a specific class of analytic equilibria. In this case we find that the pressure anisotropy and the flow may have either stabilising or destabilising effects. Also, helical configurations with small torsion and large pitch seem to have more favorable stability properties.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A discrete Nambu bracket for 2D extended Magnetohydrodynamics

In this note we propose a trilinear bracket formulation for the Hamiltonian extended Magnetohydrodynamics (XMHD) model with homogeneous mass density. The corresponding two-dimensional representation is derived by performing spatial reduction on the three-dimensional bracket, upon introducing a symmetric representation for the field variables. Subsequently, the trilinear bracket of the resulting two-dimensional, four-field model is discretized using a finite difference scheme, which results in semi-discrete dynamics that involve the Arakawa Jacobian. Simulations of planar dynamics show that this scheme respects the desired conservation properties to high precision.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Axisymmetric equilibria with pressure anisotropy and plasma flow

A generalised Grad-Shafranov equation that governs the equilibrium of an axisymmetric toroidal plasma with anisotropic pressure and incompressible flow of arbitrary direction is derived. This equation includes six free surface functions and recovers known Grad-Shafranov-like equations in the literature as well as the usual static, isotropic one. The form of the generalised equation indicates that pressure anisotropy and flow act additively on equilibrium. In addition, two sets of analytical solutions, an extended Solovev one with a plasma reaching the separatrix and an extended Hernegger-Maschke one for a plasma surrounded by a fixed boundary possessing an X-point, are constructed, particularly in relevance to the ITER and NSTX tokamaks. Furthermore, the impacts both of pressure anisotropy and plasma flow on these equilibria are examined. It turns out that depending on the maximum value and the shape of an anisotropy function, the anisotropy can act either paramagnetically or diamagnetically. Also, in most of the cases considered both the anisotropy and the flow have stronger effects on NSTX equilibria than on ITER ones.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Equilibria with incompressible flows from symmetry analysis

We identify and study new nonlinear axisymmetric equilibria with incompressible flow of arbitrary direction satisfying a generalized Grad Shafranov equation by extending the symmetry analysis presented in [G. Cicogna and F. Pegoraro, Phys. Plasmas Vol. 22, 022520 (2015)]. In particular, we construct a typical tokamak D-shaped equilibrium with peaked toroidal current density, monotonically varying safety factor and sheared electric field.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Vlasov tokamak equilibria with shearad toroidal flow and anisotropic pressure

By choosing appropriate deformed Maxwellian ion and electron distribution functions depending on the two particle constants of motion, i.e. the energy and toroidal angular momentum, we reduce the Vlasov axisymmetric equilibrium problem for quasineutral plasmas to a transcendental Grad-Shafranov-like equation. This equation is then solved numerically under the Dirichlet boundary condition for an analytically prescribed boundary possessing a lower X-point to construct tokamak equilibria with toroidal sheared ion flow and anisotropic pressure. Depending on the deformation of the distribution functions these steady states can have toroidal current densities either peaked on the magnetic axis or hollow. These two kinds of equilibria may be regarded as a bifurcation in connection with symmetry properties of the distribution functions on the magnetic axis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Tokamak-like Vlasov equilibria

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Vlasov versus reduced kinetic theories for helically symmetric equilibria

A new constant of motion for helically symmetric equilibria in the vicinity of the magnetic axis is obtained in the framework of Vlasov theory. In view of this constant of motion the Vlasov theory is compared with drift kinetic and gyrokinetic theories near axis. It turns out that as in the case of axisymmetric equilibria [H. Tasso and G. N. Throumoulopoulos, Phys. Plasmas 18, 064507 (2011)] the Vlasov current density thereon can differ appreciably from the drift kinetic and gyrokinetic current densities. This indicates 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear translational symmetric equilibria relevant to the L-H transition

Nonlinear z-independent solutions to a generalized Grad-Shafranov equation (GSE) with up to quartic flux terms in the free functions and incompressible plasma flow non parallel to the magnetic field are constructed quasi-analytically. Through an ansatz the GSE is transformed to a set of three ordinary differential equations and a constraint for three functions of the coordinate x, in cartesian coordinates (x,y), which then are solved numerically. Equilibrium configurations for certain values of the integration constants are displayed. Examination of their characteristics in connection with the impact of nonlinearity and sheared flow indicates that these equilibria are consistent with the L-H transition phenomenology. For flows parallel to the magnetic field one equilibrium corresponding to the H-state is potentially stable in the sense that a sufficient condition for linear stability is satisfied in an appreciable part of the plasma while another solution corresponding to the L-state does not satisfy the condition. The results indicate that the sheared flow in conjunction with the equilibrium nonlinearity play a stabilizing role.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On Lyapunov boundary control of unstable magnetohydrodymanic plasmas

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A comparison of Vlasov with drift kinetic and gyrokinetic theories

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lyapunov stability of flowing MHD plasmas surrounded by resistive walls

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.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Axisymmetric equilibria with anisotropic resistivity and toroidal flow

The equilibrium of an axisymmetric magnetically confined plasma with anisotropic resistivity and toroidal flow is investigated in the framework of magnetohydrodynamics (MHD). The stationary states are determined by an elliptic differential equation for the poloidal magnetic flux function, a Bernoulli equation for the pressure and two relations for the resistivity components parallel and perpendicular to the magnetic field. The equilibrium characteristics of a tokamak with rectangular cross-section are studied by means of eigenfunctions in connection with exact solutions for the cases of "compressible" flows with constant temperature, but varying density on magnetic surfaces and incompressible ones with constant density, but varying temperature thereon. Those eigenfunctions can describe either single or multiple toroidal configurations. In the former case the equilibrium has the following characteristics: (i) the profiles of the resistivity's components on the poloidal cross-section are roughly collisional (ii) the electric field perpendicular to the magnetic surfaces possesses two local extrema within the plasma and vanishes on the boundary and (iii) the toroidal current density is peaked close to the magnetic axis and vanishes on the boundary. The impact of the flow and the aspect ratio on the aforementioned quantities is evaluated for both "compressible" and incompressible flows.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Exact magnetohydrodynamic equilibria with flow and effects on the Shafranov shift

Exact solutions of the equation governing the equilibrium magetohydrodynamic states of an axisymmetric plasma with incompressible flows of arbitrary direction [H. Tasso and G.N.Throumoulopoulos, Phys. Pasmas {\bf 5}, 2378 (1998)] are constructed for toroidal current density profiles peaked on the magnetic axis in connection with the ansatz $S=-ku$, where $S=d/d u [\varrho (dΦ/du)^2]$ ($k$ is a parameter, $u$ labels the magnetic surfaces; $\varrho(u)$ and $Φ(u)$ are the density and the electrostatic potential, respectively). They pertain to either unbounded plasmas of astrophysical concern or bounded plasmas of arbitrary aspect ratio. For $k=0$, a case which includes flows parallel to the magnetic field, the solutions are expressed in terms of Kummer functions while for $k\neq 0$ in terms of Airy functions. On the basis of a tokamak solution with $k\neq 0$ describing a plasma surrounded by a perfectly conducted boundary of rectangular cross-section it turns out that the Shafranov shift is a decreasing function which can vanish for a positive value of $k$. This value is larger the smaller the aspect ratio of the configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Axisymmetric equilibria of a gravitating plasma with incompressible flows

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 $β_0\approx γ_0\approx δ_0 \approxε_0\ll 1$, where $β_0$, $γ_0$, $δ_0$, and $ε_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 $β_0\approx γ_0\approx δ_0 \approxε_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.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

On resistive magnetohydrodynamic equilibria of an axisymmetric toroidal plasma with flow

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 $ψ$ 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 $Δ^\star ψ=V_c σ$. (Here, $Δ^\star$ is the Grad-Schlüter-Shafranov operator, $σ$ is the conductivity and $V_c$ is the constant toroidal-loop voltage divided by $2 π$). 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 $σ=σ(R, ψ)$ ($R$ is the distance from the axis of symmetry) several classes of analytic equilibria with incompressible flows can be constructed having qualitatively plausible $σ$ profiles, i.e. profiles with $σ$ taking a maximum value close to the magnetic axis and a minimum value on the plasma surface. For $σ=σ(ψ)$ consideration of the relation $Δ^\starψ= V_c σ(ψ)$ 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.

0 citations0 reviewsNo code linkedOpen access