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

Title: Conservative discretization of the Landau collision integral

Abstract

Here we describe a density, momentum-, and energy-conserving discretization of the nonlinear Landau collision integral. The method is suitable for both the finite-element and discontinuous Galerkin methods and does not require structured meshes. The conservation laws for the discretization are proven algebraically and demonstrated numerically for an axially symmetric nonlinear relaxation problem using a finite-element implementation.

Authors:
 [1];  [2]
  1. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
  2. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Publication Date:
Research Org.:
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21)
OSTI Identifier:
1379746
Alternate Identifier(s):
OSTI ID: 1373960
Grant/Contract Number:
AC02-05CH11231; AC02-09CH11466
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Physics of Plasmas
Additional Journal Information:
Journal Volume: 24; Journal Issue: 3; 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; Tensor methods; Polynomials; Galerkin methods; Plasma collisions; Conservation of momentum

Citation Formats

Hirvijoki, E., and Adams, M. F.. Conservative discretization of the Landau collision integral. United States: N. p., 2017. Web. doi:10.1063/1.4979122.
Hirvijoki, E., & Adams, M. F.. Conservative discretization of the Landau collision integral. United States. doi:10.1063/1.4979122.
Hirvijoki, E., and Adams, M. F.. Tue . "Conservative discretization of the Landau collision integral". United States. doi:10.1063/1.4979122. https://www.osti.gov/servlets/purl/1379746.
@article{osti_1379746,
title = {Conservative discretization of the Landau collision integral},
author = {Hirvijoki, E. and Adams, M. F.},
abstractNote = {Here we describe a density, momentum-, and energy-conserving discretization of the nonlinear Landau collision integral. The method is suitable for both the finite-element and discontinuous Galerkin methods and does not require structured meshes. The conservation laws for the discretization are proven algebraically and demonstrated numerically for an axially symmetric nonlinear relaxation problem using a finite-element implementation.},
doi = {10.1063/1.4979122},
journal = {Physics of Plasmas},
number = 3,
volume = 24,
place = {United States},
year = {Tue Mar 28 00:00:00 EDT 2017},
month = {Tue Mar 28 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

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

Save / Share:
  • Cited by 2
  • The problem of expansion of the Boltzmann collision integral for arbitrary interaction potential is considered in the form of a power series in the momentum transferred in a collision. Such a series is termed a Landau series. An exact (without divergence) expansion of the collision integral in a Landau series is performed, the explicit form of the general term of this series is found, and its properties analyzed. (JFP)
  • We present fast numerical algorithms to solve the nonlinear Fokker-Planck-Landau equation in 3D velocity space. The discretization of the collision operator preserves the properties required by the physical nature of the Fokker-Planck-Landau equation, such as the conservation of mass, momentum, and energy, the decay of the entropy, and the fact that the steady states are Maxwellians. At the end of this paper, we give numerical results illustrating the efficiency of these fast algorithms in terms of accuracy and CPU time. 20 refs., 7 figs.
  • Boundary-element methods (BEMs) for solving integral equations numerically have been used in many fields to compute the induced charges at dielectric boundaries. In this paper, we consider a more accurate implementation of BEM in the context of ions in aqueous solution near proteins, but our results are applicable more generally. The ions that modulate protein function are often within a few angstroms of the protein, which leads to the significant accumulation of polarization charge at the protein-solvent interface. Computing the induced charge accurately and quickly poses a numerical challenge in solving a popular integral equation using BEM. In particular, themore » accuracy of simulations can depend strongly on seemingly minor details of how the entries of the BEM matrix are calculated. We demonstrate that when the dielectric interface is discretized into flat tiles, the qualocation method of Tausch et al. [IEEE Trans Comput.-Comput.-Aided Des. 20, 1398 (2001)] to compute the BEM matrix elements is always more accurate than the traditional centroid-collocation method. Qualocation is not more expensive to implement than collocation and can save significant computational time by reducing the number of boundary elements needed to discretize the dielectric interfaces.« less