×

Additive Schwarz-based fully coupled implicit methods for resistive Hall magnetohydrodynamic problems. (English) Zbl 1343.76049

Summary: A parallel, fully coupled, nonlinearly implicit Newton-Krylov-Schwarz algorithm is proposed for the numerical simulation of a magnetic reconnection problem described by a system of resistive Hall magnetohydrodynamics equations in slab geometry. A key component of the algorithm is a restricted additive Schwarz preconditioner defined for problems with doubly periodic boundary conditions. We show numerically that with such a preconditioned nonlinearly implicit method the time step size is no longer constrained by the CFL number or the convergence of the Newton solver. We report the parallel performance of the algorithm and software on machines with thousands of processors.

MSC:

76M25 Other numerical methods (fluid mechanics) (MSC2010)
76W05 Magnetohydrodynamics and electrohydrodynamics
76X05 Ionized gas flow in electromagnetic fields; plasmic flow

Software:

PETSc; MPI
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] S. Balay, K. Buschelman, V. Eijkhout, W.D. Gropp, D. Kaushik, M. Knepley, L.C. McInnes, B.F. Smith, H. Zhang, PETSc Users Manual, ANL-95/11 - Revision 2.3.1, Argonne National Laboratory, 2006.; S. Balay, K. Buschelman, V. Eijkhout, W.D. Gropp, D. Kaushik, M. Knepley, L.C. McInnes, B.F. Smith, H. Zhang, PETSc Users Manual, ANL-95/11 - Revision 2.3.1, Argonne National Laboratory, 2006.
[2] Biskamp, D., Nonlinear Magnetohydrodynamics (1993), Cambridge University Press: Cambridge University Press New York
[3] Braginskii, S. I., Reviews of Plasma Physics, vol.1 (1965), Consultants Bureau: Consultants Bureau New York
[4] Cafaro, E.; Grasso, D.; Pegoraro, F.; Porcelli, F.; Saluzzi, A., Invariants and geometric structures in nonlinear Hamiltonian magnetic reconnection, Phys. Rev. Lett., 80, 20, 4430-4433 (1998)
[5] Cai, X.-C., Additive Schwarz algorithms for parabolic convection-diffusion equations, Numer. Math., 60, 41-62 (1990) · Zbl 0737.65078
[6] Cai, X.-C.; Gropp, W. D.; Keyes, D. E.; Melvin, R. G.; Young, D. P., Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation, SIAM J. Sci. Comput., 19, 246-265 (1998) · Zbl 0917.76035
[7] Cai, X.-C.; Gropp, W. D.; Keyes, D. E.; Tidriri, M. D., Newton-Krylov-Schwarz methods in CFD, (Rannacher, R., Proceedings of the International Workshop on the Navier-Stokes Equations. Proceedings of the International Workshop on the Navier-Stokes Equations, Notes in Numerical Fluid Mechanics (1994), Vieweg Verlag: Vieweg Verlag Braunschweig) · Zbl 0876.76059
[8] Cai, X.-C.; Keyes, D. E., Nonlinearly preconditioned inexact Newton algorithms, SIAM J. Sci. Comput., 24, 183-200 (2002) · Zbl 1015.65058
[9] Cai, X.-C.; Sarkis, M., A restricted additive Schwarz preconditioner for general sparse linear systems, SIAM J. Sci. Comput., 21, 792-797 (1999) · Zbl 0944.65031
[10] Chacón, L.; Knoll, D. A., A 2D high-β Hall MHD implicit nonlinear solver, J. Comput. Phys., 188, 573-592 (2003) · Zbl 1127.76375
[11] Chacón, L.; Knoll, D. A.; Finn, J. M., An implicit, nonlinear reduced resistive MHD solver, J. Comput. Phys., 178, 15-36 (2002) · Zbl 1139.76328
[12] Coleman, T. F.; Moré, J. J., Estimation of sparse Jacobian matrices and graph coloring problems, SIAM J. Numer. Anal., 20, 187-209 (1983) · Zbl 0527.65033
[13] Dennis, J. E.; Schnabel, R. B., Numerical Methods for Unconstrained Optimization and Nonlinear Equations (1996), SIAM: SIAM Philadelphia · Zbl 0847.65038
[14] Edwards, A.; Campbell, D.; Engelhardt, W.; Fahrbach, H.; Gill, R.; Granetz, R.; Tsuji, S.; Tubbing, B.; Weller, A.; Wesson, J.; Zasche, D., Rapid collapse of a plasma sawtooth oscillation in the JET tokamak, Phys. Rev. Lett., 57, 2, 210-213 (1986)
[15] Fitzpatrick, R., Scaling of forced magnetic reconnection in the Hall-magnetohydrodynamical Taylor problem with arbitrary guide field, Phys. Plasmas, 11, 8, 3961-3968 (2004)
[16] Fitzpatrick, R., Scaling of forced magnetic reconnection in the Hall-magnetohydrodynamic Taylor problem, Phys. Plasmas, 11, 3, 937-946 (2004)
[17] Germaschewski, K.; Bhattacharjee, A.; Grauer, R.; Keyes, D.; Smith, B., Using Krylov-Schwarz methods in an adaptive mesh refinement environment, (Proceedings of the Chicago Workshop on adaptive mesh refinement methods. Proceedings of the Chicago Workshop on adaptive mesh refinement methods, Lecture Notes in Computational Science and Engineering, vol. 41 (2003), Springer), 115-124 · Zbl 1246.65228
[18] Goldston, R. J.; Rutherford, P. H., Introduction to Plasmas Physics (1995), IOP Publishing
[19] Grasso, D.; Pegoraro, F.; Porcelli, F.; Califano, F., Hamiltonian magnetic reconnection, Plasma Phys. Control Fusion, 41, 1497-1515 (1999)
[20] Gropp, W.; Lusk, E.; Skjellum, A., Using MPI: Portable Parallel Programming with the Message-passing Interface (1994), The MIT Press: The MIT Press Cambridge
[21] Hairer, E.; Nørsett, S.; Wanner, G., Solving Ordinary Differential Equations I (1993), Springer-Verlag: Springer-Verlag New York · Zbl 0789.65048
[22] Hwang, F.-N.; Cai, X.-C., A parallel nonlinear additive Schwarz preconditioned inexact Newton algorithm for incompressible Navier-Stokes equations, J. Comput. Phys., 204, 666-691 (2005) · Zbl 1267.76083
[23] Jothiprasad, G.; Mavriplis, D.; Caughey, D., Higher-order time integration schemes for the unsteady Navier-Stokes equations on unstructured meshes, J. Comput. Phys., 191, 542-566 (2003) · Zbl 1134.76428
[24] Ladyzhenskaya, O. A., The Boundary Value Problems of Mathematical Physics (1985), Springer-Verlag: Springer-Verlag New York · Zbl 0588.35003
[25] Lewis, J.; Rehm, R., The numerical solution of a nonseparable elliptic partial differential equation by preconditioned conjugate gradients, J. Res. Nat. Bur. Standards, 85, 5, 367-390 (1980) · Zbl 0444.65070
[26] Ottaviani, M.; Porcelli, F., Nonlinear collisionless magnetic reconnection, Phys. Rev. Lett., 71, 23, 3802-3805 (1993)
[27] Ottaviani, M.; Porcelli, F., Fast nonlinear magnetic reconnection, Phys. Plasmas, 2, 11, 4104-4117 (1995)
[28] Reynolds, D. R.; Samtaney, R.; Woodward, C. S., A fully implicit numerical method for single-fluid resistive magnetohydrodynamics, J. Comput. Phys., 219, 144-162 (2006) · Zbl 1103.76036
[29] Saad, Y., Iterative Methods for Sparse Linear Systems (2003), SIAM · Zbl 1002.65042
[30] Smith, B.; Bjørstad, P.; Gropp, W., Domain Decomposition. Parallel Multilevel Methods for Elliptic Partial Differential Equations (1996), Cambridge University Press: Cambridge University Press New York · Zbl 0857.65126
[31] Strauss, H. R.; Longcope, D. W., An adaptive finite element method for magnetohydrodynamics, J. Comput. Phys., 147, 318-336 (1998) · Zbl 0936.76032
[32] Toselli, A.; Widlund, O., Domain Decomposition Methods - Algorithms and Theory (2005), Springer-Verlag: Springer-Verlag Berlin · Zbl 1069.65138
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.