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

This content will become publicly available on January 22, 2020

Title: High-order two-fluid plasma solver for direct numerical simulations of plasma flows with full transport phenomena

Abstract

The two-fluid plasma equations for a single ion species, with full transport terms, including temperature and magnetic field dependent ion and electron viscous stresses and heat fluxes, frictional drag force, and ohmic heating terms, have been implemented in the CFDNS code and solved by using sixth-order non-dissipative compact finite differences for plasma flows in several different regimes. In order to be able to fully resolve all the dynamically relevant time and length scales, while maintaining computational feasibility, the assumptions of infinite speed of light and negligible electron inertia have been made. Non-dimensional analysis of the two-fluid plasma equations shows that, by varying the characteristic/background number density, length scale, temperature, and magnetic strength, the corresponding Hall, resistive, and ideal magnetohydrodynamic equations can be recovered as limiting cases. The accuracy and robustness of this two-fluid plasma solver in handling plasma flows in different regimes have been validated against four canonical problems: Alfven and whistler dispersion relations, electromagnetic plasma shock, and magnetic reconnection. For all test cases, by using physical dissipation and diffusion, with negligible numerical dissipation/diffusion, fully converged Direct Numerical Simulation (DNS)-like solutions are obtained when the ion Reynolds number based on the grid size is smaller than a threshold value whichmore » is about 2.3. For the magnetic reconnection problem, the results show that the magnetic flux saturation time and value converge when the ion and magnetic Reynolds numbers are large enough. Thus, the DNS-like results become relevant to practical problems with much larger Reynolds numbers.« less

Authors:
ORCiD logo [1]; ORCiD logo [2]
  1. Texas A & M Univ., Corpus Christi, TX (United States). Dept. of Engineering
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA); LANL Laboratory Directed Research and Development (LDRD) Program
OSTI Identifier:
1503195
Report Number(s):
LA-UR-18-31228
Journal ID: ISSN 1070-664X
Grant/Contract Number:  
89233218CNA000001
Resource Type:
Accepted Manuscript
Journal Name:
Physics of Plasmas
Additional Journal Information:
Journal Volume: 26; Journal Issue: 1; Journal ID: ISSN 1070-664X
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Citation Formats

Li, Z., and Livescu, D. High-order two-fluid plasma solver for direct numerical simulations of plasma flows with full transport phenomena. United States: N. p., 2019. Web. doi:10.1063/1.5082190.
Li, Z., & Livescu, D. High-order two-fluid plasma solver for direct numerical simulations of plasma flows with full transport phenomena. United States. doi:10.1063/1.5082190.
Li, Z., and Livescu, D. Tue . "High-order two-fluid plasma solver for direct numerical simulations of plasma flows with full transport phenomena". United States. doi:10.1063/1.5082190.
@article{osti_1503195,
title = {High-order two-fluid plasma solver for direct numerical simulations of plasma flows with full transport phenomena},
author = {Li, Z. and Livescu, D.},
abstractNote = {The two-fluid plasma equations for a single ion species, with full transport terms, including temperature and magnetic field dependent ion and electron viscous stresses and heat fluxes, frictional drag force, and ohmic heating terms, have been implemented in the CFDNS code and solved by using sixth-order non-dissipative compact finite differences for plasma flows in several different regimes. In order to be able to fully resolve all the dynamically relevant time and length scales, while maintaining computational feasibility, the assumptions of infinite speed of light and negligible electron inertia have been made. Non-dimensional analysis of the two-fluid plasma equations shows that, by varying the characteristic/background number density, length scale, temperature, and magnetic strength, the corresponding Hall, resistive, and ideal magnetohydrodynamic equations can be recovered as limiting cases. The accuracy and robustness of this two-fluid plasma solver in handling plasma flows in different regimes have been validated against four canonical problems: Alfven and whistler dispersion relations, electromagnetic plasma shock, and magnetic reconnection. For all test cases, by using physical dissipation and diffusion, with negligible numerical dissipation/diffusion, fully converged Direct Numerical Simulation (DNS)-like solutions are obtained when the ion Reynolds number based on the grid size is smaller than a threshold value which is about 2.3. For the magnetic reconnection problem, the results show that the magnetic flux saturation time and value converge when the ion and magnetic Reynolds numbers are large enough. Thus, the DNS-like results become relevant to practical problems with much larger Reynolds numbers.},
doi = {10.1063/1.5082190},
journal = {Physics of Plasmas},
number = 1,
volume = 26,
place = {United States},
year = {2019},
month = {1}
}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on January 22, 2020
Publisher's Version of Record

Save / Share:

Works referenced in this record:

Self-Consistent Equations Including Exchange and Correlation Effects
journal, November 1965