Source author record

M. V. Umansky

M. V. Umansky appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.plasm-ph physics.comp-ph

Catalog footprint

What is connected

7works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2021arXiv

INGRID: an interactive grid generator for 2D edge plasma modeling

A fusion boundary-plasma domain is defined by axisymmetric magnetic surfaces where the geometry is often complicated by the presence of one or more X-points; and modeling boundary plasmas usually relies on computational grids that account for the magnetic field geometry. The new grid generator INGRID (Interactive Grid Generator) presented here is a Python-based code for calculating grids for fusion boundary plasma modeling, for a variety of configurations with one or two X-points in the domain. Based on a given geometry of the magnetic field, INGRID first calculates a skeleton grid which consists of a small number of quadrilateral patches; then it puts a subgrid on each of the patches, and joins them in a global grid. This domain partitioning strategy makes possible a uniform treatment of various configurations with one or two X-points in the domain. This includes single-null, double-null, and other configurations with two X-points in the domain. The INGRID design allows generating grids either interactively, via a parameter-file driven GUI, or using a non-interactive script-controlled workflow. Results of testing demonstrate that INGRID is a flexible, robust, and user-friendly grid-generation tool for fusion boundary-plasma modeling.

preprint2014arXiv

BOUT++: Recent and current developments

BOUT++ is a 3D nonlinear finite-difference plasma simulation code, capable of solving quite general systems of PDEs, but targeted particularly on studies of the edge region of tokamak plasmas. BOUT++ is publicly available, and has been adopted by a growing number of researchers worldwide. Here we present improvements which have been made to the code since its original release, both in terms of structure and its capabilities. Some recent applications of these methods are reviewed, and areas of active development are discussed. We also present algorithms and tools which have been developed to enable creation of inputs from analytic expressions and experimental data, and for processing and visualisation of output results. This includes a new tool Hypnotoad for the creation of meshes from experimental equilibria. Algorithms have been implemented in BOUT++ to solve a range of linear algebraic problems encountered in the simulation of reduced MHD and gyro-fluid models: A preconditioning scheme is presented which enables the plasma potential to be calculated efficiently using iterative methods supplied by the PETSc library, without invoking the Boussinesq approximation. Scaling studies are also performed of a linear solver used as part of physics-based preconditioning to accelerate the convergence of implicit time-integration schemes.

preprint2013arXiv

Energy dynamics in a simulation of LAPD turbulence

Energy dynamics calculations in a 3D fluid simulation of drift wave turbulence in the linear Large Plasma Device (LAPD) [W. Gekelman et al., Rev. Sci. Inst. 62, 2875 (1991)] illuminate processes that drive and dissipate the turbulence. These calculations reveal that a nonlinear instability dominates the injection of energy into the turbulence by overtaking the linear drift wave instability that dominates when fluctuations about the equilibrium are small. The nonlinear instability drives flute-like ($k_\parallel = 0$) density fluctuations using free energy from the background density gradient. Through nonlinear axial wavenumber transfer to $k_\parallel \ne 0$ fluctuations, the nonlinear instability accesses the adiabatic response, which provides the requisite energy transfer channel from density to potential fluctuations as well as the phase shift that causes instability. The turbulence characteristics in the simulations agree remarkably well with experiment. When the nonlinear instability is artificially removed from the system through suppressing $k_\parallel=0$ modes, the turbulence develops a coherent frequency spectrum which is inconsistent with experimental data.

preprint2013arXiv

Nonlinear instability in simulations of Large Plasma Device turbulence

Several simulations of turbulence in the Large Plasma Device (LAPD) [W. Gekelman et al., Rev. Sci. Inst. 62, 2875 (1991)] are energetically analyzed and compared with each other and with the experiment. The simulations use the same model, but different axial boundary conditions. They employ either periodic, zero-value, zero-derivative, or sheath axial boundaries. The linear stability physics is different between the scenarios because the various boundary conditions allow the drift wave instability to access different axial structures, and the sheath boundary simulation contains a conducting wall mode instability which is just as unstable as the drift waves. Nevertheless, the turbulence in all the simulations is relatively similar because it is primarily driven by a robust nonlinear instability that is the same for all cases. The nonlinear instability preferentially drives $k_\parallel = 0$ potential energy fluctuations, which then three-wave couple to $k_\parallel \ne 0$ potential energy fluctuations in order to access the adiabatic response to transfer their energy to kinetic energy fluctuations. The turbulence self-organizes to drive this nonlinear instability, which destroys the linear eigenmode structures, making the linear instabilities ineffective.

preprint2011arXiv

Modeling of plasma turbulence and transport in the Large Plasma Device

Numerical simulation of plasma turbulence in the Large Plasma Device (LAPD) [Gekelman et al, Rev. Sci. Inst., 62, 2875, 1991] is presented. The model, implemented in the BOUndary Turbulence (BOUT) code [M. Umansky et al, Contrib. Plasma Phys. 180, 887 (2009)], includes 3-D collisional fluid equations for plasma density, electron parallel momentum, and current continuity, and also includes the effects of ion-neutral collisions. In nonlinear simulations using measured LAPD density profiles but assuming constant temperature profile for simplicity, self-consistent evolution of instabilities and nonlinearly-generated zonal flows results in a saturated turbulent state. Comparisons of these simulations with measurements in LAPD plasmas reveal good qualitative and reasonable quantitative agreement, in particular in frequency spectrum, spatial correlation and amplitude probability distribution function of density fluctuations. For comparison with LAPD measurements, the plasma density profile in simulations is maintained either by direct azimuthal averaging on each time step, or by adding particle source/sink function. The inferred source/sink values are consistent with the estimated ionization source and parallel losses in LAPD. These simulations lay the groundwork for more a comprehensive effort to test fluid turbulence simulation against LAPD data.

preprint2010arXiv

Simulation of Edge Localised Modes using BOUT++

The BOUT++ code is used to simulate ELMs in a shifted circle equilibrium. Reduced ideal MHD simulations are first benchmarked against the linear ideal MHD code ELITE, showing good agreement. Diamagnetic drift effects are included finding the expected suppression of high toroidal mode number modes. Nonlinear simulations are performed, making the assumption that the anomalous kinematic electron viscosity is comparable to the anomalous electron thermal diffusivity. This allows simulations with realistically high Lundquist numbers S = 1e8, finding ELM sizes of 5-10% of the pedestal stored thermal energy. Scans show a strong dependence of the ELM size resistivity at low Lundquist numbers, with higher resistivity leading to more violent eruptions. At high Lundquist numbers relevant to high-performance discharges, ELM size is independent of resistivity as hyper-resistivity becomes the dominant dissipative effect.

preprint2008arXiv

BOUT++: a framework for parallel plasma fluid simulations

A new modular code called BOUT++ is presented, which simulates 3D fluid equations in curvilinear coordinates. Although aimed at simulating Edge Localised Modes (ELMs) in tokamak X-point geometry, the code is able to simulate a wide range of fluid models (magnetised and unmagnetised) involving an arbitrary number of scalar and vector fields, in a wide range of geometries. Time evolution is fully implicit, and 3rd-order WENO schemes are implemented. Benchmarks are presented for linear and non-linear problems (the Orszag-Tang vortex) showing good agreement. Performance of the code is tested by scaling with problem size and processor number, showing efficient scaling to thousands of processors. Linear initial-value simulations of ELMs using reduced ideal MHD are presented, and the results compared to the ELITE linear MHD eigenvalue code. The resulting mode-structures and growth-rate are found to be in good agreement (BOUT++ = 0.245, ELITE = 0.239). To our knowledge, this is the first time dissipationless, initial-value simulations of ELMs have been successfully demonstrated.