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

Title: Plasma transport in an Eulerian AMR code

Abstract

A plasma transport model has been implemented in an Eulerian AMR radiation-hydrodynamics code, xRage, which includes plasma viscosity in the momentum tensor, viscous dissipation in the energy equations, and binary species mixing with consistent species mass and energy fluxes driven by concentration gradients, ion and electron baro-diffusion terms and temperature gradient forces. The physics basis, computational issues, numeric options, and results from several test problems are discussed. The transport coefficients are found to be relatively insensitive to the kinetic correction factors when the concentrations are expressed with the molar fractions and the ion mass differences are large. The contributions to flow dynamics from plasma viscosity and mass diffusion were found to increase significantly as scale lengths decrease in an inertial confinement fusion relevant Kelvin-Helmholtz instability mix layer. The mixing scale lengths in the test case are on the order of 100 μm and smaller for viscous effects to appear and 10 μm or less for significant ion species diffusion, evident over durations on the order of nanoseconds. The temperature gradient driven mass flux is seen to deplete a high Z tracer ion at the ion shock front. The plasma transport model provides the generation of the atomic mix per unitmore » of interfacial area between two species with no free parameters. The evolution of the total atomic mix then depends also on an accurate resolution or estimate of the interfacial area between the species mixing by plasma transport. High resolution simulations or a more Lagrangian-like treatment of species interfaces may be required to distinguish plasma transport and numerical diffusion in an Eulerian computation of complex and dynamically evolving mix regions.« less

Authors:
ORCiD logo [1];  [1];  [1];  [1]; ORCiD logo [1]; ORCiD logo [1]
  1. 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) Advanced Simulation and Computing (ASC) Program; USDOE
OSTI Identifier:
1361797
Alternate Identifier(s):
OSTI ID: 1357127; OSTI ID: 1420619
Report Number(s):
LA-UR-16-29421
Journal ID: ISSN 1070-664X
Grant/Contract Number:
AC52-06NA25396
Resource Type:
Journal Article: Published Article
Journal Name:
Physics of Plasmas
Additional Journal Information:
Journal Volume: 24; Journal Issue: 4; 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; plasma transport Eulerian code; plasma transport; Eulerian code

Citation Formats

Vold, E. L., Rauenzahn, R. M., Aldrich, C. H., Molvig, K., Simakov, A. N., and Haines, B. M. Plasma transport in an Eulerian AMR code. United States: N. p., 2017. Web. doi:10.1063/1.4979171.
Vold, E. L., Rauenzahn, R. M., Aldrich, C. H., Molvig, K., Simakov, A. N., & Haines, B. M. Plasma transport in an Eulerian AMR code. United States. doi:10.1063/1.4979171.
Vold, E. L., Rauenzahn, R. M., Aldrich, C. H., Molvig, K., Simakov, A. N., and Haines, B. M. Tue . "Plasma transport in an Eulerian AMR code". United States. doi:10.1063/1.4979171.
@article{osti_1361797,
title = {Plasma transport in an Eulerian AMR code},
author = {Vold, E. L. and Rauenzahn, R. M. and Aldrich, C. H. and Molvig, K. and Simakov, A. N. and Haines, B. M.},
abstractNote = {A plasma transport model has been implemented in an Eulerian AMR radiation-hydrodynamics code, xRage, which includes plasma viscosity in the momentum tensor, viscous dissipation in the energy equations, and binary species mixing with consistent species mass and energy fluxes driven by concentration gradients, ion and electron baro-diffusion terms and temperature gradient forces. The physics basis, computational issues, numeric options, and results from several test problems are discussed. The transport coefficients are found to be relatively insensitive to the kinetic correction factors when the concentrations are expressed with the molar fractions and the ion mass differences are large. The contributions to flow dynamics from plasma viscosity and mass diffusion were found to increase significantly as scale lengths decrease in an inertial confinement fusion relevant Kelvin-Helmholtz instability mix layer. The mixing scale lengths in the test case are on the order of 100 μm and smaller for viscous effects to appear and 10 μm or less for significant ion species diffusion, evident over durations on the order of nanoseconds. The temperature gradient driven mass flux is seen to deplete a high Z tracer ion at the ion shock front. The plasma transport model provides the generation of the atomic mix per unit of interfacial area between two species with no free parameters. The evolution of the total atomic mix then depends also on an accurate resolution or estimate of the interfacial area between the species mixing by plasma transport. High resolution simulations or a more Lagrangian-like treatment of species interfaces may be required to distinguish plasma transport and numerical diffusion in an Eulerian computation of complex and dynamically evolving mix regions.},
doi = {10.1063/1.4979171},
journal = {Physics of Plasmas},
number = 4,
volume = 24,
place = {United States},
year = {Tue Apr 04 00:00:00 EDT 2017},
month = {Tue Apr 04 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1063/1.4979171

Citation Metrics:
Cited by: 2works
Citation information provided by
Web of Science

Save / Share:
  • A plasma transport model has been implemented in an Eulerian AMR radiation-hydrodynamics code, xRage, which includes plasma viscosity in the momentum tensor, viscous dissipation in the energy equations, and binary species mixing with consistent species mass and energy fluxes driven by concentration gradients, ion and electron baro-diffusion terms and temperature gradient forces. The physics basis, computational issues, numeric options, and results from several test problems are discussed. The transport coefficients are found to be relatively insensitive to the kinetic correction factors when the concentrations are expressed with the molar fractions and the ion mass differences are large. The contributions tomore » flow dynamics from plasma viscosity and mass diffusion were found to increase significantly as scale lengths decrease in an inertial confinement fusion relevant Kelvin-Helmholtz instability mix layer. The mixing scale lengths in the test case are on the order of 100 μm and smaller for viscous effects to appear and 10 μm or less for significant ion species diffusion, evident over durations on the order of nanoseconds. The temperature gradient driven mass flux is seen to deplete a high Z tracer ion at the ion shock front. The plasma transport model provides the generation of the atomic mix per unit of interfacial area between two species with no free parameters. The evolution of the total atomic mix then depends also on an accurate resolution or estimate of the interfacial area between the species mixing by plasma transport. High resolution simulations or a more Lagrangian-like treatment of species interfaces may be required to distinguish plasma transport and numerical diffusion in an Eulerian computation of complex and dynamically evolving mix regions.« less
  • Cited by 2
  • A Three-Dimensional Finite Volume Arbitrary Lagrangian-Eulerian simulation code was developed to study different plasma physics problems in 3D+t. The code is based on a complex multi-component species program with transport and radiation terms written and applied to plasma and fusion physics problems. Three different examples are shown: double-base chemical propellant combustion, ignition and propagation of a thermonuclear detonation wave, and, the development of the Kelvin-Helmholtz (KH) instability in local plane slab models of the magnetopause, showing the response of a background equilibrium to the excitation by finite amplitude perturbations generated upstream.
  • The Ranchero Magnetic Flux Compression Generator (FCG) has been used to create current pulses in the 10-100 MA range for driving both “static” low inductance (0.5 nH) loads1 for generator demonstration purposes and high inductance (10-20 nH) imploding liner loads for ultimate use in physics experiments at very high energy density. This report offers a description of these efforts.
  • The Ranchero Magnetic Flux Compression Generator (FCG) has been used to create current pulses in the 10-­100 MA range for driving both “static” low inductance (0.5 nH) loads1 for generator demonstration purposes and high inductance (10-­20 nH) imploding liner loads2 for ultimate use in physics experiments at very high energy density. Simulations of the standard Ranchero generator have recently shown that it had a design issue that could lead to flux trapping in the generator, and a non-­ robust predictability in its use in high energy density experiments. A re-­examination of the design concept for the standard Ranchero generator, promptedmore » by the possible appearance of an aneurism at the output glide plane, has led to a new generation of Ranchero generators designated the RancheroS (for swooped). This generator has removed the problematic output glide plane and replaced it with a region of constantly increasing diameter in the output end of the FCG cavity in which the armature is driven outward under the influence of an additional HE load not present in the original Ranchero. The resultant RancheroS generator, to be tested in LA43S-­L13, probably in early FY17, has a significantly increased initial inductance and may be able to drive a somewhat higher load inductance than the standard Ranchero. This report will use the Eulerian AMR code Roxane to study the ability of the new design to drive static loads, with a goal of providing a database corresponding to the load inductances for which the generator might be used and the anticipated peak currents such loads might produce in physics experiments. Such a database, combined with a simple analytic model of an ideal generator, where d(LI)/dt = 0, and supplemented by earlier estimates of losses in actual use of the standard Ranchero, scaled to estimate the increase in losses due to the longer current carrying perimeter in the RancheroS, can then be used to bound the expectations for the current drive one may apply to any load assembly in future experiments.« less