×

PHOENIX: MHD spectral code for rotating laboratory and gravitating astrophysical plasmas. (English) Zbl 1310.76189

Summary: The new PHOENIX code is discussed together with a sample of many new results that are obtained concerning magnetohydrodynamic (MHD) spectra of axisymmetric plasmas where flow and gravity are consistently taken into account. PHOENIX, developed from the CASTOR code [W. Kerner et al., J. Comput. Phys. 142, No. 2, 271–303, Art. No. CP985910 (1998; Zbl 0921.76100)], incorporates purely toroidal, or both toroidal and poloidal flow and external gravitational fields to compute the entire ideal or resistive MHD spectrum for general tokamak or accretion disk configurations. These equilibria are computed by means of FINESSE [A. J. C. Beliën et al., J. Comput. Phys. 182, No. 1, 91–117 (2002; Zbl 1021.76026)], which discriminates between the different elliptic flow regimes that may occur. PHOENIX makes use of a finite element method in combination with a spectral method for the discretization. This leads to a large generalized eigenvalue problem, which is solved by means of Jacobi-Davidson algorithm [G. L. G. Sleijpen and H. A. van der Vorst, SIAM J. Matrix Anal. Appl. 17, No. 2, 401–425 (1996; Zbl 0860.65023)].
PHOENIX is compared with CASTOR, PEST-1 and ERATO for an internal mode of Soloviev equilibria. Furthermore, the resistive internal kink mode has been computed to demonstrate that the code can accurately handle small values for the resistivity. A new reference test case for a Soloviev-like equilibrium with toroidal flow shows that, on a particular unstable mode, the flow has a quantifiable stabilizing effect regardless of the direction of the flow. PHOENIX reproduces the Toroidal Flow induced Alfvén Eigenmode (TFAE, [B. van der Holst et al., “New Alfvén continuum gaps and global modes induced by toroidal flow”, Phys. Rev. Lett. 84, No. 13, 2865–2868 (2000; doi:10.1103/PhysRevLett.84.2865)]) where finite resistivity in combination with equilibrium flow effects causes resonant damping. Localized ideal gap modes are presented for tokamak plasmas with toroidal and poloidal flow. Finally, we demonstrate the ability to spectrally diagnose magnetized accretion disk equilibria where gravity acts together with either purely toroidal flow or both toroidal and poloidal flow. These cases show that the MHD continua can be unstable or overstable due to the presence of a gravitational field together with equilibrium flow-driven dynamics [J. P. Goedbloed et al., “Unstable continuous spectra of transonic axisymmetric plasmas”, Phys. Plasmas 11, No. 1, 28–54 (2004; doi:10.1063/1.1631812)].

MSC:

76X05 Ionized gas flow in electromagnetic fields; plasmic flow
85A30 Hydrodynamic and hydromagnetic problems in astronomy and astrophysics
85-04 Software, source code, etc. for problems pertaining to astronomy and astrophysics
76M25 Other numerical methods (fluid mechanics) (MSC2010)
76U05 General theory of rotating fluids
76W05 Magnetohydrodynamics and electrohydrodynamics
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Beliën, A. J.C.; Botchev, M. A.; Goedbloed, J. P.; van der Holst, B.; Keppens, R., J. Comp. Phys., 182, 91 (2002)
[2] Bondeson, A.; Iacono, R.; Bhattacharjee, A., Phys. Fluids, 30, 2167 (1987)
[3] A. Bondeson, G. Vlad, H. Lütjens, Controlled fusion and plasma heating, in: Proceedings of the 17th European Conference, Amsterdam, The Netherlands (European Physical Society, Petit-Lancy, 1990) B Part II 14, 1990, p. 906.; A. Bondeson, G. Vlad, H. Lütjens, Controlled fusion and plasma heating, in: Proceedings of the 17th European Conference, Amsterdam, The Netherlands (European Physical Society, Petit-Lancy, 1990) B Part II 14, 1990, p. 906.
[4] Bondeson, A.; Vlad, G.; Lütjens, H., Phys. Fluids B, 4, 1889 (1992)
[5] C. Busch, in: Proceedings of 32nd EPS Conference on Plasma Physics, ECA VCol C, vol. 29, 2005, P-1.016.; C. Busch, in: Proceedings of 32nd EPS Conference on Plasma Physics, ECA VCol C, vol. 29, 2005, P-1.016.
[6] Chandra, D.; Sen, A.; Kaw, P.; Bora, M. P.; Kruger, S., Nucl. Fusion, 45, 524 (2005)
[7] Cheng, C. Z.; Chance, M. S., J. Comp. Phys., 71, 124 (1987)
[8] Crombé, K., Phys. Rev. Lett., 95, 155003 (2005)
[9] Goedbloed, J. P., Physica, 53, 501 (1971)
[10] Goedbloed, J. P., Phys. Fluids, 18, 1258 (1975)
[11] Goedbloed, J. P., Space Sci. Rev., 107, 353 (2003)
[12] Goedbloed, J. P., Phys. Scripta T, 107, 159 (2004)
[13] Goedbloed, J. P.; Beliën, A. J.C.; van der Holst, B.; Keppens, R., AIP Conf. Proc., 703, 42 (2004)
[14] Goedbloed, J. P.; Beliën, A. J.C.; van der Holst, B.; Keppens, R., Phys. Plasmas, 11, 28 (2004)
[15] Goedbloed, J. P.; Beliën, A. J.C.; van der Holst, B.; Keppens, R., Phys. Plasmas, 11, 4332 (2004)
[16] Goedbloed, J. P.; Sakanaka, P. H., Phys. Fluids, 17, 908 (1974)
[17] Golub, G. H.; van Loan, C. F., Matrix Computations (1996), The Johns Hopkins University Press: The Johns Hopkins University Press Baltimore, MD · Zbl 0865.65009
[18] H. Grad, H. Rubin, in: Proceedings of the Second UN International Conference on the Peaceful Uses of Atomic Energy, vol. 31, 1958, p. 190.; H. Grad, H. Rubin, in: Proceedings of the Second UN International Conference on the Peaceful Uses of Atomic Energy, vol. 31, 1958, p. 190.
[19] Grimm, R. C.; Greene, J. M.; Johnson, J. L., Methods Comput. Phys., 16, 253 (1976)
[20] Gruber, R.; Rappaz, J., Finite Element Methods in Linear Ideal Magnetohydrodynamics (1985), Springer: Springer New York · Zbl 0573.76001
[21] Gruber, R.; Troyon, F.; Berger, D.; Bernard, L. C.; Rousset, S.; Schreider, R.; Kerner, W.; Schneider, W.; Roberts, K. V., Comp. Phys. Commun., 21, 323 (1981)
[22] Guazzotto, L.; Betti, R.; Manickham, J.; Kaye, S., Phys. Plasmas, 11, 604 (2004)
[23] Hain, K.; Lüst, R., Z. Naturforsch., 13a, 936 (1958)
[24] E. Hameiri, Ph.D. thesis, New York University, 1976.; E. Hameiri, Ph.D. thesis, New York University, 1976.
[25] Hameiri, E., Phys. FLuids, 26, 230 (1983)
[26] Hirshman, S. P., Nucl. Fusion, 18, 917 (1978)
[27] van der Holst, B.; Beliën, A. J.C.; Goedbloed, J. P., Phys. Rev. Lett., 84, 2865 (2000)
[28] van der Holst, B.; Beliën, A. J.C.; Goedbloed, J. P., Phys. Plasmas, 7, 4208 (2000)
[29] Hsu, C. T.; Shaing, K. C.; Gormley, R., Phys. Plasmas, 1, 132 (1996)
[30] G.T.A. Huysmans, J.P. Goedbloed, W. Kerner, in: Proceedings of CP90 Conf. on Comp. Phys. Proc., 1991, p. 371.; G.T.A. Huysmans, J.P. Goedbloed, W. Kerner, in: Proceedings of CP90 Conf. on Comp. Phys. Proc., 1991, p. 371.
[31] Huysmans, G. T.A.; Goedbloed, J. P.; Kerner, W., Phys. Fluids B, 5, 1545 (1993)
[32] Keppens, R.; Casse, F.; Goedbloed, J. P., ApJ, 569, L121 (2002)
[33] Kerner, W.; Goedbloed, J. P.; Huysmans, G. T.A.; Poedts, S.; Schwarz, E., J. Comput. Phys., 142, 271 (1998)
[34] Morris, R. C.; Haines, M. G.; Hastie, R. J., Phys. Plasmas, 3, 4513 (1996)
[35] Nool, M.; van der Ploeg, A., SIAM J. Sci. Comput., 22, 95 (2000)
[36] Pao, Y. P., Nucl. Fusion, 15, 631 (1975)
[37] Poedts, S.; Kerner, W., Phys. Rev. Lett., 66, 2871 (1991)
[38] Poedts, S.; Schwarz, E., J. Comp. Phys., 105, 165 (1993)
[39] Rappaz, J., Numer. Math., 28, 15 (1977)
[40] Saad, Y.; Schultz, M. H., SIAM J. Sci. Stat. Comput., 7, 856 (1986)
[41] Semenzato, S.; Gruber, R.; Zehrfeld, H. P., Comput. Phys. Rep., 1, 389 (1984)
[42] Shafranov, V. D., Sov. Phys. JETP, 6, 545 (1958)
[43] Shaing, K. C.; Hirshman, S. P., Phys. Fluids B, 1, 705 (1989)
[44] Sleijpen, G. L.G.; van der Vorst, H. A., SIAM J. Matrix Anal. Appl., 17, 401 (1996)
[45] Soloviev, L. S., Sov. Phys. JETP, 26, 400 (1968)
[46] Stix, T. H., Phys. Fluids, 16, 1260 (1973)
[47] Strang, G.; Fix, G. J., Analysis of the finite element method (1973), Prentice Hall: Prentice Hall Englewood, NJ · Zbl 0278.65116
[48] Strumberger, E.; Günter, S.; Merkel, P.; Riondato, S.; Schwarz, E.; Tichmann, C.; Zehrfeld, H. P., Nucl. Fusion, 45, 1156 (2005)
[49] Taguchi, M., Phys. Plasmas Control Fusion, 33, 859 (1991)
[50] Zehrfeld, H. P.; Green, B. J., Nucl. Fusion, 12, 569 (1972)
[51] Zelazny, R.; Stankiewicz, R.; Galkowski, A.; Potempski, S., Plasma Phys. Control Fusion, 35, 1215 (1993)
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. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.