Paper Info
Title
Nested Iteration and First-Order System Least Squares for Incompressible, Resistive Magnetohydrodynamics
Abstract
This paper develops a nested iteration algorithm to solve time-dependent nonlinear systems of partial differential equations. For each time step, Newton's method is used to form approximate solutions from a sequence of nested spaces, where the resolution of the approximations increases as the algorithm progresses. Nested iteration results in most of the iterations being performed on coarser grids, where minimal work is needed to reduce error to the level of discretization error. The approximate solution on a given coarse grid is interpolated to a refined grid and is used as an initial guess for the problem posed there. The approximation is then already close enough to the solution on the current grid that a minimal amount of work is needed to solve the refined problem due to the rapid convergence of Newton's method near a solution. The paper develops an algorithm that attempts to optimize accuracy-per-computational-cost on each grid, so that essentially no unnecessary work is done on any grid. The nested iteration algorithm is then applied to a reduced two-dimensional model of the incompressible, resistive magnetohydrodynamic (MHD) equations. Using this algorithm on the MHD equations in the context of a first-order system least squares finite element discretization and algebraic multigrid to solve the linearized systems, instabilities in a model tokamak fusion reactor are simulated. Numerical results show that this highly complex nonlinear problem is solved in an equivalent of 30-80 fine-grid relaxation sweeps per time step.
Year
DOI
Venue
2010
10.1137/090766905
SIAM J. Scientific Computing
Keywords
Field
DocType
magnetohydrodynamics,refined grid,resistive magnetohydrodynamics,nested iteration,algebraic multigrid,coarser grid,nested iteration result,nested iteration algorithm,approximate solution,coarse grid,complex nonlinear problem,fosls,minimal work,time step,first-order system,current grid,aspect ratio,nonlinear system,adaptive mesh refinement,charged particles,finite element,plasma physics,iterative algorithm,partial differential equation,mathematics,steady state,linear system
Discretization,Mathematical optimization,Nonlinear system,Mathematical analysis,Finite element method,Numerical analysis,Mathematics,Grid,Domain decomposition methods,Multigrid method,Newton's method
Journal
Volume
Issue
ISSN
32
3
1064-8275
Citations 
PageRank 
References 
9
0.73
11
Authors
5
Name
Order
Citations
PageRank
J. H. Adler15610.02
T. A. Manteuffel227838.19
S. F. McCormick326638.47
J. Ruge429333.76
G. D. Sanders5101.09