Discrete Calderon's projections on parallelepipeds and their application to computing exterior magnetic fields for FRC plasmas

Confining dense plasma in a field reversed configuration (FRC) is considered a promising approach to fusion. Numerical simulation of this process requires setting artificial boundary conditions (ABCs) for the magnetic field because whereas the plasma itself occupies a bounded region (within the FRC coils), the field extends from this region all the way to infinity. If the plasma is modeled using single fluid magnetohydrodynamics (MHD), then the exterior magnetic field can be considered quasi-static. This field has a scalar potential governed by the Laplace equation. The quasi-static ABC for the magnetic field is obtained using the method of difference potentials, in the form of a discrete Calderon boundary equation with projection on the artificial boundary shaped as a parallelepiped. The Calderon projection itself is computed by convolution with the discrete fundamental solution on the three-dimensional Cartesian grid.

[1]  Uri Shumlak,et al.  A discontinuous Galerkin method for the full two-fluid plasma model , 2005, Comput. Phys. Commun..

[2]  Stefaan Poedts,et al.  Principles of Magnetohydrodynamics: With Applications to Laboratory and Astrophysical Plasmas , 2004 .

[3]  Rony Keppens,et al.  Advanced Magnetohydrodynamics: With Applications to Laboratory and Astrophysical Plasmas , 2010 .

[4]  Donald B. Melrose,et al.  Electromagnetic processes in dispersive media : a treatment based on the dielectric tensor , 1991 .

[5]  Semyon Tsynkov,et al.  Active Shielding and Control of Noise , 2001, SIAM J. Appl. Math..

[6]  N. K. Kulman,et al.  Method of difference potentials and its applications , 2001 .

[7]  T. Hussey,et al.  Three-dimensional magnetic field enhancement in a liner implosion system , 1995 .

[8]  J. Freidberg,et al.  Plasma Physics and Fusion Energy , 2007 .

[9]  Viktor S. Ryaben’kii,et al.  Boundary equations with projections , 1985 .

[10]  P. Morse,et al.  Methods of theoretical physics , 1955 .

[11]  Semyon Tsynkov,et al.  The Method of Difference Potentials for the Helmholtz Equation Using Compact High Order Schemes , 2012, Journal of Scientific Computing.

[12]  I. Tsukerman.  A Singularity-Free Boundary Equation Method for Wave Scattering , 2010, IEEE Transactions on Antennas and Propagation.

[13]  Harold Grad,et al.  HYDROMAGNETIC EQUILIBRIA AND FORCE-FREE FIELDS , 1958 .

[14]  Semyon Tsynkov,et al.  External Boundary Conditions for Three-Dimensional Problems of Computational Aerodynamics , 1999, SIAM J. Sci. Comput..

[15]  David R. Owen,et al.  an Introduction to Finite Element Method , 1979 .

[16]  Nathan Ida,et al.  Introduction to the Finite Element Method , 1997 .

[17]  Robert T. Seeley,et al.  SINGULAR INTEGRALS AND BOUNDARY VALUE PROBLEMS. , 1966 .

[18]  Uri Shumlak,et al.  Approximate Riemann solver for the two-fluid plasma model , 2003 .

[19]  Uri Shumlak,et al.  Advanced physics calculations using a multi-fluid plasma model , 2011, Comput. Phys. Commun..

[20]  Semyon Tsynkov,et al.  On the Definition of Surface Potentials for Finite-Difference Operators , 2003, J. Sci. Comput..

[21]  Sergei Utyuzhnikov.  Generalized Calderon - Ryaben'kii's Potentials , 2008 .

[22]  V. Shafranov.  Plasma Equilibrium in a Magnetic Field , 1966 .

[23]  I. Furno,et al.  A high-density field reversed configuration plasma for magnetized target fusion , 2004, IEEE Transactions on Plasma Science.

[24]  A. Tikhonov,et al.  Equations of Mathematical Physics , 1964 .

[25]  R. E. Peterkin,et al.  Transport of Magnetic Flux in an Arbitrary Coordinate ALE Code , 1998 .

[26]  S. Tsynkov.  Numerical solution of problems on unbounded domains. a review , 1998 .

[27]  Carlos E. Kenig,et al.  Boundary value problems for elliptic equations , 1991 .

[28]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .