×

zbMATH — the first resource for mathematics

Calculation of the vacuum Green’s function valid even for high toroidal mode numbers in tokamaks. (English) Zbl 1129.78014
Summary: The present evaluation of the Green’s function used for the magnetic scalar potential in vacuum calculations for axisymmetric geometry has been found to be deficient even for moderately high, \(n\), the toroidal mode number. This is relevant to the edge localized peeling-ballooning modes calculated by gato, pest and other mhd stability codes. The deficiency was due to the loss of numerical precision arising from the upward recursion relation used for generating the functions from the values at \(n = 0\) from the complete elliptic integrals of the first and second kinds. To ameliorate this, a direct integration of the integral representation of the function is crafted to achieve the necessary high accuracy for moderately high mode numbers, with due consideration to the singular behavior of the integrand involved. At higher mode numbers the loss of numerical precision due to cancellations from the oscillatory behavior of the integrand is further avoided by judiciously deforming the integration contour into the complex plane to obtain a new integral representation for the Green’s function. Near machine precision, roughly 12–16 digits, can be achieved by using a combination of these techniques. The relation to the associated Legendre functions, as well as a novel integral representation of these are also described.

MSC:
78A55 Technical applications of optics and electromagnetic theory
76W05 Magnetohydrodynamics and electrohydrodynamics
65D30 Numerical integration
65D15 Algorithms for approximation of functions
Software:
ERATO
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Troyon, F.; Bernard, L.C.; Gruber, R., Comput. phys. commun., 19, 161, (1980)
[2] Gruber, R.; Troyon, F.; Berger, D.; Bernard, L.C.; Rousset, S.; Schreiber, R.; Kerner, W.; Schneider, W.; Roberts, K.V., Comput. phys. commun., 21, 323, (1981)
[3] Gruber, R.; Semenzato, S.; Troyon, F.; Tsunematsu, T.; Kerner, W.; Merkel, P.; Schneider, W., Comput. phys. commun., 24, 363, (1981)
[4] Grimm, R.C.; Greene, J.M.; Johnson, J.L., Computation of the magnetohydrodynamic spectrum in axisymmetric toroidal confinement systems, (), 253-280
[5] Chance, M.S., Phys. plasmas, 4, 2161, (1997)
[6] Wilson, H.R.; Snyder, P.B.; Miller, R.L.; Huysmans, G.T.A., Phys. plasmas, 1277, (2002)
[7] Snyder, P.B.; Wilson, H.R., Phys. plasmas, 2037, (2002)
[8] M. Abramovitz, I.A. Stegun, Handbook of Mathematical Functions, ninth ed., Number 55 in Applied Mathematics Series, National Bureau of Standards, Superintendent of Documents, U.S. Government Printing Office, Washington, D.C. 20402, 1970.
[9] Press, W.H.; Flannery, B.P.; Teukolsky, S.A.; Vetterling, W.T., Numerical recipes: the art of scientific computing, (1988), Cambridge University Press · Zbl 0661.65001
[10] Bulirsch, R., Numer. math., 7, 78, (1965)
[11] Dahlquist, G.; Björck, Å., Numerical methods, (1974), Prentice-Hall Inc. Englewood Cliffs, New Jersey, pp. 297-301 (chapter 7)
[12] Carlson, B.C., SIAM J. math. anal., 9, 525, (1978)
[13] Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P., Numerical recipes in Fortran 77: the art of scientific and parallel computing, (1996), Cambridge University Press · Zbl 0878.68049
[14] Dwight, H.B., Tables of integrals and other mathematical data, (1966), The Macmillan Company New York · Zbl 0154.18410
[15] Whittaker, E.T.; Watson, G.N., A course of modern analysis, (1963), Cambridge University Press England · Zbl 0108.26903
[16] Erdélyi, A., Higher transcendental functions, The Bateman manuscript project, vol. 1, (1953), McGraw-Hill Book Co. Inc. New York, Toronto, Canada, London, England, (chapter 3)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.