skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A second-order semi-implicit δf method for hybrid simulation

Abstract

A second-order accurate semi-implicit Lorentz force ions, fluid electrons δf hybrid model has been developed using a “current closure” scheme. The model assumes quasi-neutrality and is fully electromagnetic. The implicit field solver improves numerical accuracy by separating the equilibrium terms in the presence of small perturbations. The equilibrium part of the generalized Ohm’s law is solved by direct matrix inversion along the direction of gradients for every Fourier mode in the other two directions, while the nonlinear part is solved iteratively. The simulation has been benchmarked on Alfvén waves, ion sound waves and whistler waves against analytical dispersion relation in a slab. In particular, the first-order and second-order schemes are compared by studying the numerical damping of whistler waves. The full evolution of the resistive tearing mode using the Harris sheet equilibrium is also investigated. The linear growth rate and mode structure are compared with the resistive MHD theory. Important tearing mode nonlinear phenomena such as the Rutherford regime and saturation are demonstrated. We also presented systematic study of Rutherford growth rates and saturation island width, which is consistent with previous MHD studies.

Authors:
; ; ;
Publication Date:
OSTI Identifier:
22233609
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 245; Other Information: Copyright (c) 2013 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICAL METHODS AND COMPUTING; COMPARATIVE EVALUATIONS; DISPERSION RELATIONS; ELECTRONS; EQUILIBRIUM; HYBRIDIZATION; IONS; ITERATIVE METHODS; LORENTZ FORCE; MAGNETOHYDRODYNAMICS; NONLINEAR PROBLEMS; PERTURBATION THEORY; SIMULATION; TEARING INSTABILITY; WHISTLERS

Citation Formats

Cheng, Jianhua, E-mail: Jianhua.Cheng@Colorado.Edu, Parker, Scott E., Chen, Yang, and Uzdensky, Dmitri A.. A second-order semi-implicit δf method for hybrid simulation. United States: N. p., 2013. Web. doi:10.1016/J.JCP.2013.03.017.
Cheng, Jianhua, E-mail: Jianhua.Cheng@Colorado.Edu, Parker, Scott E., Chen, Yang, & Uzdensky, Dmitri A.. A second-order semi-implicit δf method for hybrid simulation. United States. doi:10.1016/J.JCP.2013.03.017.
Cheng, Jianhua, E-mail: Jianhua.Cheng@Colorado.Edu, Parker, Scott E., Chen, Yang, and Uzdensky, Dmitri A.. Mon . "A second-order semi-implicit δf method for hybrid simulation". United States. doi:10.1016/J.JCP.2013.03.017.
@article{osti_22233609,
title = {A second-order semi-implicit δf method for hybrid simulation},
author = {Cheng, Jianhua, E-mail: Jianhua.Cheng@Colorado.Edu and Parker, Scott E. and Chen, Yang and Uzdensky, Dmitri A.},
abstractNote = {A second-order accurate semi-implicit Lorentz force ions, fluid electrons δf hybrid model has been developed using a “current closure” scheme. The model assumes quasi-neutrality and is fully electromagnetic. The implicit field solver improves numerical accuracy by separating the equilibrium terms in the presence of small perturbations. The equilibrium part of the generalized Ohm’s law is solved by direct matrix inversion along the direction of gradients for every Fourier mode in the other two directions, while the nonlinear part is solved iteratively. The simulation has been benchmarked on Alfvén waves, ion sound waves and whistler waves against analytical dispersion relation in a slab. In particular, the first-order and second-order schemes are compared by studying the numerical damping of whistler waves. The full evolution of the resistive tearing mode using the Harris sheet equilibrium is also investigated. The linear growth rate and mode structure are compared with the resistive MHD theory. Important tearing mode nonlinear phenomena such as the Rutherford regime and saturation are demonstrated. We also presented systematic study of Rutherford growth rates and saturation island width, which is consistent with previous MHD studies.},
doi = {10.1016/J.JCP.2013.03.017},
journal = {Journal of Computational Physics},
number = ,
volume = 245,
place = {United States},
year = {Mon Jul 15 00:00:00 EDT 2013},
month = {Mon Jul 15 00:00:00 EDT 2013}
}
  • Cited by 1
  • An iterative implicit-explicit hybrid scheme is proposed for hyperbolic systems of conservation laws. Each wave in a system may be implicitly, or explicitly, or partially implicitly and partially explicitly treated depending on its associated Courant number in each numerical cell, and the scheme is able to smoothly switch between implicit and explicit calculations. The scheme is of Godunov-type in both explicit and implicit regimes, is in a strict conservation form, and is accurate to second-order in both space and time for all Courant numbers. The computer code for the scheme is easy to vectorize. Multicolors proposed in this paper maymore » reduce the number of iterations required to reach a converged solution by several orders for a large time step. The feature of the scheme is shown through numerical examples. 38 refs., 12 figs.« less
  • A semi-implicit method for solving the full compressible resistive magnetohydrodynamic equations in three dimensions is presented. The method is designed for use in the modeling of fusion plasmas in magnetic confinement devices. The method is unconditionally stable with respect to the fast compressional modes. The time step is limited instead by the slower shear Alfven motion. The computing time required for one time step is essentially the same as that for explicit methods. The code is applied to resistive instabilities in cylindrical tokamak equilibria.
  • Time domain simulation of plasmas in the rf time scale range is difficult because the time-scale is long compared to the electron plasma wave period, and in addition, the various cutoff and resonance behaviors within the plasma insure that any explicit finite-difference scheme would be numerically unstable. We resolve this dilemma with a new algorithm [1] based upon a semi-implicit method (e.g., explicit Maxwell, implicit plasma), such that all linear plasma dispersion behavior are faithfully reproduced at the available temporal and spatial resolution, despite the fact that the simulation time-step may exceed the electron gyro and electron plasma time scalesmore » by orders of magnitude. This new algorithm is now available in the Vorpal parallel computing simulation framework, which can provide complex boundary modeling of edge geometry such as antenna, limiters, and other 3-D structures. The result is a unified model joining the power delivery system to the bulk plasma. We report on benchmarking of the new capability for several classical benchmarks, including tunneling through low density edge plasma, cyclotron resonance, and mode-conversion to ICW (ion cyclotron wave). We also report on initial efforts to include 3-D edge geometry, including wave launcher with realistic edge density profiles.« less
  • An iterative implicit-explicit hybrid scheme is proposed for one-dimensional ideal magnetohydrodynamical (MHD) equations. The scheme is accurate to second order in both space and time for all Courant numbers. Each kind of MHD wave, both compressible and incompressible, may be either implicitly or explicitly, or partially implicitly and partially explicitly treated depending on their associated Courant numbers in each numerical cell, and the scheme is able to smoothly switch between explicit and implicit calculations. The scheme is of Godunov type in both explicit and implicit regimes, in which the flux needed in a Godunov scheme is calculated from Riemann problems,more » is in a strictly conservation form, and is able to resolve MHD discontinuities. Only a single level of iterations is involved in the scheme, and the iterations solve both the implicit relations arising from upstream centered differences for all wave families and the nonlinearity of MHD equations. Multicolors proposed in the paper may significantly reduce the number of iterations required to reach a converged solution. Compared with the largest Courant number in a simulation, only a small number of iterations are needed in the scheme. The feature of the scheme has been shown through numerical examples. It is easy to vectorize the computer code of the scheme.« less