Towards a scalable fully-implicit fully-coupled resistive MHD formulation with stabilized FE methods

Journal of Computational Physics

Volumn 229, Issue 20, 2010, Pages 7649-7671

  Shadid, J N ^a   Pawlowski, R P ^a   Banks, J W ^b   Chacon L ^c   Lin, P T ^a   Tuminaro, R S ^a  

^a SANDIA NATIONAL LABORATORIES   (United States)

^b LAWRENCE LIVERMORE NATIONAL LABORATORY   (United States)

^c OAK RIDGE NATIONAL LABORATORY   (United States)

Author keywords

Direct to steady state; Fully implicit; Large scale parallel; Multilevel preconditioner; Newton Krylov; Resistive MHD; Stabilized finite element

Indexed keywords

ALGEBRA; FINITE ELEMENT METHOD; LINEAR SYSTEMS; PUMPS;

DIRECT-TO-STEADY-STATE; FULLY-COUPLED; FULLY-IMPLICIT; LARGE-SCALE PARALLEL; LARGE-SCALES; MULTILEVEL PRECONDITIONERS; NEWTON-KRYLOV; RESISTIVE MHD; STABILIZED FINITE ELEMENT; STEADY STATE;

MAGNETOHYDRODYNAMICS;

EID: 77955511672     PISSN: 00219991     EISSN: 10902716     Source Type: Journal    
DOI: 10.1016/j.jcp.2010.06.018     Document Type: Article

References (99)

1. Goedbloed H., Poedts S. Principles of Magnetohydrodynamics with Applications to Laboratory and Astrophysical Plasmas 2004, Cambridge University Press.
2. Ascher U.M., Petzold L.R. Computer methods for ordinary differential equations and differential-algebraic equations. SIAM 1998.
3. Dai W., Woodward P.R. A simple finite difference scheme for multi-dimensional magnetohydrodynamic equations. J. Comput. Phys. 1998, 142:331.
4. Ryu D., Miniati F., Jones T.W., Frank A. A divergence-free upwind code for multi-dimensional magnetohydrodynamics flows. Astrophys. J. 1998, 509:244.
5. Balsara D.S. Divergence-free adaptive mesh refinement for magnetohydrodynamics. J. Comput. Phys. 2001, 174:614-648.
6. Tóth G. The ∇·B=0 constraint in shock-capturing magnetohydrodynamics codes. J. Comput. Phys. 2000, 161:605-652.
7. Tóth G., Keppens R., Botchev M.A. Implicit and semi-implicit schemes in the Versatile Advection Code: numerical tests. Astron. Astrophys. 1998, 332:1159-1170.
8. Keppens R., Tóth G., Botchev M.A., Ploeg A.V.D. Implicit and semi-implicit schemes: algorithms. Int. J. Numer. Meth. Fl. 1999, 30:335-352.
9. Aydemir A.Y., Barnes D.C. An implicit algorithm for compressible three-dimensional magnetohydrodynamic calculations. J. Comput. Phys. 1985, 59(1):108-119.
10. Park W., Breslau J., Chen J., Fu G.Y., Jardin S.C., Klasky S., Menard J., Pletzer A., Stratton B.C., Stutman D., Strauss H.R., Sugiyama L.E. Nonlinear simulation studies of Tokamaks and STS. Nucl. Fusion 2003, 43(6):483-489.
11. Jardin S.C., Breslau J.A. Implicit solution of the four-field extended-magnetohydrodynamic equations using high-order high-continuity finite elements. Phys. Plasmas 2005, 12(5):056101.
12. Harned D.S., Kerner W. Semi-implicit method for three-dimensional compressible magnetohydrodynamic simulation. J. Comput. Phys. 1985, 60:62-75.
13. Schnack D.D., Barnes D.C., Harned D.S., Caramana E.J. Semi-implicit magnetohydrodynamic calculations. J. Comput. Phys. 1987, 70:330-354.
14. Harned D.S., Mikic Z. Accurate semi-implicit treatment of the Hall effect in magnetohydrodynamic computations. J. Comput. Phys. 1989, 83:1-15.
15. Sovinec C.R., Glasser A.H., Gianakon T.A., Barnes D.C., Nebel R.A., Kruger S.E., Schnack D.D., Plimpton S.J., Tarditi A., Chu M.S. Nonlinear magnetohydrodynamics simulation using high-order finite elements. J. Comput. Phys. 2004, 195(1):355-386. The NIMROD team.
16. Hujeirat A. IRMHD: an implicit radiative and magnetohydrodynamical solver for self-gravitating systems. Mon. Not. Roy. Astron. Soc. 1998, 298:310-320.
17. A.C. Robinson, et al., ALEGRA: An arbitrary Lagrangian-Eulerian multimaterial, multiphysics code, in: AIAA 2008-1235 46th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, 2008.
18. G. Toth, R. Keppens, Versatile advection code, http://www.phys.uu.nl/~toth/.
19. Knoll D.A., Chacón L., Margolin L., Mousseau V.A. On balanced approximations for the time integration of multiple time-scale systems. J. Comput. Phys. 2003, 185:583-611.
20. Hujeirat A., Rannacher R. On the efficiency and robustness of implicit methods in computational astrophysics. New Astron. Rev. 2001, 45:425-447.
21. Reynolds D.R., Samtaney R., Woodward C.S. A fully implicit numerical method for single-fluid resistive magnetohydrodynamics. J. Comput. Phys. 2006, 219(1):144-162.
22. Ovtchinnikov S., Dobrian F., Cai X.-C., Keyes D. Additive Schwarz-based fully coupled implicit methods for resistive Hall magnetohydrodynamic problems. J. Comput. Phys. 2007, 225:1919-1936.
23. Chacón L., Knoll D.A., Finn J.M. Implicit, nonlinear reduced resistive MHD nonlinear solver. J. Comput. Phys. 2002, 178(1):15-36.
24. Chacón L., Knoll D.A. A 2D high-β Hall MHD implicit nonlinear solver. J. Comput. Phys. 2003, 188(2):573-592.
25. Chacón L. An optimal, parallel, fully implicit Newton-Krylov solver for three-dimensional visco-resistive magnetohydrodynamics. Phys. Plasmas 2008, 15:056103.
26. Chacón L. Scalable solvers for 3D magnetohydrodynamics. J. Phys: Conf. Ser. 2008, 125:012041.
27. Reynolds D.R., Samtaney R., Woodward C.S. Operator-based preconditioning of stiff hyperbolic systems. SIAM J. Sci. Comput. 2010, 32:150-170.
28. Salah N.B., Soulaimani A., Habashi W.G., Fortin M. A conservative stabilized finite element method for the magnetohydrodynamics equations. Int. J. Num. Meth. Fuilds 1999, 29:535-554.
29. Brezzi F. On existence, uniqueness and approximation of saddle-point problems arising from Lagrange multipliers. RAIRO Model. Math. Anal. Numer. 1974, 21:129-151.
30. Gunzburger M. Finite Element Methods for Viscous Incompressible Flows 1989, Academic Press, Boston.
31. Barth T., Bochev P., Gunzburger M., Shadid J. A taxonomy of consistently stabilized finite element methods for the Stokes problem. SIAM J. Sci. Comput. 2004, 25(5):1585-1607.
32. Brezzi F., Pitkaranta J. On the stabilization of finite element approximations of the stokes problem, Efficient Solutions of Elliptic Systems. Vieweg Notes on Numerical Fluid Mechanics 1984, 10:11-19. Vieweg, Wiesbaden.
33. Salah N.B., Soulaimani A., Habashi W.G. A finite element method for magnetohydrodynamics. Comput. Meth. Appl. M. 2001, 190:5867-5892.
34. Codina R., Hernandez-Silva N. Stabilized finite element approximation of the stationary magneto-hydrodynamics equations. Comput. Mech. 2006, 38(4-5):344-355.
35. Gerbeau J.-F. A stabilized finite element method for the incompressible magnetohydrodynamic equations. Numer. Math. 2000, 87:83-111.
36. Lankalapalli S., Flaherty J., Shephard M., Strauss H. An adaptive finite element method for magnetohydrodynamics. J. Comput. Phys. 2007, 225:363-381.
37. Brown P.N., Saad Y. Convergence theory of nonlinear Newton-Krylov algorithms. SIAM J. Optimiz. 1994, 4:297-330.
38. Eisenstat S., Walker H. Globally convergent inexact Newton methods. SIAM J. Optimiz. 1994, 4:393-422.
39. Quarteroni A., Valli A. Domain Decomposition Methods for Partial Differential Equations 1999, Oxford University Press, Oxford.
40. Sala M., Shadid J.N., Tuminaro R.S. An improved convergence bound for aggregation-based domain decomposition preconditioners. SIAM J. Matrix Anal. A. 2006, 27(3):744-756.
41. Lin P.T., Sala M., Shadid J.N., Tuminaro R.S. Performance of fully-coupled algebraic multilevel domain decomposition preconditioners for incompressible flow and transport. Int. J. Numer. Meth. Eng. 2006, 67(9):208-225.
42. Strauss H.R. Nonlinear, 3-dimensional magnetohydrodynamics of noncircular tokamaks. Phys. Fluids 1976, 19(1):134-140.
43. Hazeltine R.D., Kotschenreuther M., Morrison P.J. A four-field model for tokamak plasma dynamics. Phys. Fluids 1985, 28(8):2466-2477.
44. Drake J.F., Antonsen T.M. Nonlinear reduced fluid equations for toroidal plasmas. Phys. Fluids 1984, 27(4):898-908.
45. Moreau R. Magnetohydrodynamics 1990, Kluwer, Dordrecht.
46. Davidson P.A. An Introduction to Magnetohydrodynamics 2001, Cambridge University Press.
47. Donea J., Huerta A. Finite Element Methods for Flow Problems 2002, John Wiley.
48. Dedner A., Kemm F., Kroner D., Munz C.-D., Schnitzer T., Wesenberg M. Hyperbolic divergence cleaning for the MHD equations. J. Comput. Phys. 2002, 175:645-673.
49. Chacón L. A non-staggered, conservative, ∇·B=0, finite-volume scheme for 3D implicit extended magnetohydrodynamics in curvilinear geometries. Comput. Phys. Commun. 2004, 163:143-171.
50. Jackson J.D. Classical Electrodynamics 1975, John Wiley & Sons. 2nd ed.
51. Brooks A.N., Hughes T. Streamline upwind/petrov-galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations. Comput. Meth. Appl. M. 1982, 32:199-259.
52. Hughes T., Franca L., Balestra M. A new finite element formulation for computational fluid dynamics: V. Circumventing the Babuska-Brezzi condition: A stable Petrov-Galerkin formulation of the Stokes problem accommodating equal-order interpolations. Comput. Meth. Appl. M. 1986, 59:85-99.
53. Hughes T., Mallet M. A new finite element formulation for computational fluid dynamics: III. The generalized streamline operator for multidimentional advective-diffusive systems. Comput. Meth. Appl. M. 1986, 58:305-328.
54. F. Shakib, Finite element analysis of the compressible Euler and Navier-Stokes equations, Ph.D. Thesis, Division of Applied Mathematics, Stanford University, 1989.
55. Hughes T. Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods. Comput. Meth. Appl. M. 1995, 127:387-401.
56. Dennis J.E., Schnabel R.B. Numerical methods for unconstrained optimization and nonlinear equations. Series in Automatic Computation 1983, Prentice-Hall, Englewood Cliffs, NJ.
57. Pawlowski R.P., Shadid J.N., Simonis J.P., Walker H.F. Globalaization techniques for Newton-Krylov methods and applications to the fully-coupled solution of the Navier-Stokes equations. SIAM Rev. 2006, 48:700-721.
58. Dembo R.S., Eisenstat S.C., Steihaug T. Inexact Newton methods. SIAM J. Numer. Anal. 1982, 19:400-408.
59. Eisenstat S.C., Walker H.F. Choosing the forcing terms in an inexact Newton method. SIAM J. Sci. Comput. 1996, 17:16-32.
60. Heroux M., Bartlett R., Howle V., Hoekstra R., Hu J., Kolda T., Lehoucq R., Long K., Pawlowski R., Phipps E., Salinger A., Thornquist H., Tuminaro R., Willenbring J., Williams A. An Overview of Trilinos Project. ACM Trans. Math. Software 2005, 31(3):397-423.
61. Saad Y. Iterative Methods for Sparse Linear Systems 2003, SIAM.
62. Axelsson O. Iterative Solution Methods 1994, Cambridge University Press, New York.
63. Tuminaro R., Tong C., Shadid J., Devine K.D., Day D. On a multilevel preconditioning module for unstructured mesh Krylov solvers: two-level Schwarz. Comm. Num. Method. Eng. 2002, 18:383-389.
64. Shadid J., Tuminaro R., Devine K., Henningan G., Lin P. Performance of fully-coupled domain decomposition preconditioners for finite element transport/reaction simulations. J. Comput. Phys. 2005, 205(1):24-47.
65. Trottenberg U., Oosterlee C., Schüller A. Multigrid 2001, Academic Press, London.
66. Shadid J.N., Salinger A.G., Pawlowski R.P., Lin P.T., Hennigan G.L., Tuminaro R.S., Lehoucq R.B. Stabilized FE computational analysis of nonlinear steady state transport/reaction systems. Comput. Meth. Appl. M. 2006, 195:1846-1871.
67. Ruge J., Stüben K. Algebraic multigrid (AMG). Frontiers in Applied Mathematics 1987, vol. 3:73-130. SIAM, Philadelphia, PA. S.F. McCormick (Ed.).
68. G. Karypis, V. Kumar, Parallel multilevel k-way partitioning scheme for irregular graphs, in: ACM/IEEE Proceedings of SC96: High Performance Networking and Computing, 1996. Article No. 35.
69. Vaněk P., Mandel J., Brezina M. Algebraic multigrid by smoothed aggregation for second and fourth order elliptic problems. Computing 1996, 56:179-196.
70. Vaněk P., Brezina M., Mandel J. Convergence of algebraic multigrid based on smoothed aggregation. Numer. Math. 2001, 88:559-579.
71. Sala M., Tuminaro R. A new Petrov-Galerkin smoothed aggregation preconditioner for nonsymmetric linear systems. SIAM J. Sci. Comput. 2008, 31(1):143-166.
72. M. Gee, C. Siefert, J. Hu, R. Tuminaro, M. Sala, ML 5.0 smoothed aggregation user's guide, Technical Report. SAND2006-2649, Sandia National Laboratories, Albuquerque, NM 87185, 2006.
73. Davis T. Direct methods for sparse linear systems 2006, SIAM, Philadelphia, PA.
74. M. Sala, Amesos 2.0 reference guide, Technical Report. SAND2004-4820, Sandia National Laboratories, September 2004.
75. Lin P.T., Sala M., Shadid J.N., Tuminaro R.S. Performance of a geometric and an algebraic multilevel preconditioner for incompressible flow and transport. Computational Mechanics: WCCM VI in conjunction with APCOM'04 2004, Tsinghus University Press/Springer-Verlag.
76. Lin P.T., Shadid J.N., Sala M., Tuminaro R., Hennigan G.L., Hoekstra R.J. Performance of a parallel algebraic multilevel preconditioner for stabilized finite element semiconductor device modeling. J. Comput. Phys. 2009, 228(17):6079-6616.
77. Moffat H.K. Magnetic Field Generation in Electrically Conducting Fluids 1983, Cambridge University Press.
78. Codina R., Hernadez-Silva N. Stabilized finite element approximation of the stationary magnetohydrodynamics equations. Comput. Mech. 2006, 38:344-355.
79. Ainsworth M., Oden J.T. A Posteriori Error Estimation in Finite Element Analysis 2000, John Wiley.
80. Batchelor G.K. An Introduction to Fluid Mechanics 1967, Cambridge University Press.
81. Fadeev V.M., Kvartskhava I.F., Komarov N.N. Self-focusing of local plasma currents. Nucl. Fusion 1965, 5(3):202-209.
82. Hughes M., Pericleous K.A., Cross M. The numerical modeling of dc electromagnetic pump and brake flow. Appl. Math. Model. 1995, 19:713-723.
83. Chandrasekhar S. Hydrodynamic and Hydromagnetic Stability 1961, Oxford University Press.
84. Biskamp D. Magnetic Reconnection in Plasmas 2000, Cambridge University Press, Cambridge, UK.
85. Knoll D.A., Chacón L. Coalescence of magnetic islands, sloshing, and the pressure problem. Phys. Plasmas 2006, 13(3):32307-32311.
86. Priest E., Forbes T. Magnetic Reconnection: MHD Theory and Applications 2006, Cambridge University Press.
87. Shadid J. A fully-coupled Newton-Krylov solution method for parallel unstructured finite element fluid flow, heat and mass transfer simulations. Int. J. CFD 1999, 12:199-211.
88. Bell J.B., Colella P., Glaz H.M. A second-order projection method for the incompressible Navier-Stokes equations. J. Comput. Phys. 1989, 85:257.
89. Chorin A.J. A numerical method for solving incompressible viscous problems. J. Comput. Phys. 1967, 2:12.
90. Gresho P.M. On the theory of semi-implicit projection methods for viscous incompressible flow and its implementation via a finite element method that also introduces a nearly consistent mass matrix. part 1: Theory. Int. J. Num. Meth. Fluids 1990, 11:587-620.
91. Oran E.S., Boris J.P. Numerical Simulation of Reactive Flow 2001, Cambridge Universty Press, Cambridge.
92. Guermond J.L., Quartapelle L. On stability of convergence of projection methods based on pressure poisson equation. Int. J. Numer. Meth. Fl. 1998, 26:1039-1053.
93. Karniadakis G.E., Israeli M., Orszag S.A. High-order splitting methods for the incompressible Navier-Stokes equations. J. Comput. Phys. 1991, 97:414-443.
94. Perot J.B. An analysis of the fractional step method. J. Comput. Phys. 1993, 108:51-58.
95. Strikwerda J.C., Lee Y.S. The accuracy of the fractional step method. SIAM J. Numer. Anal. 1999, 37(1):37-47.
96. Patankar S.V. Numerical heat transfer and fluid flow 1980, Hemisphere Public Corporation, New York.
97. Issa R. Solution of the implicitly discretized fluid flows equations by operator splitting. J. Comp. Phys 1986, 62(1):40-65.
98. Deng G.B., Piquet J., Queutey P., Visonneau M.A. A new fully coupled solution of the Navier-Stokes equations. Int. J. Numer. Meth. Fl. 1994, 19:605-639.
99. Deng G.B., Piquet J., Vasseur X., Visonneau M.A. A new fully coupled method for computing turbulent flows. Comput. Fluids 2001, 30:445-472.