Date of Award

Spring 1-1-2014

Document Type


Degree Name

Doctor of Philosophy (PhD)


Applied Mathematics

First Advisor

Thomas A. Manteuffel

Second Advisor

Stephen F. McCormick

Third Advisor

John W. Ruge

Fourth Advisor

Harihar Rajaram

Fifth Advisor

Luis Chacon


Efforts are currently being directed towards a fully implicit, electromagnetic, JFNK-based solver, motivating the necessity of developing a fluid-based, electromagnetic, preconditioning strategy [15]. The two-fluid plasma (TFP) model is an ideal approximation to the kinetic Jacobian. The TFP model couples both an ion and an electron fluid with Maxwell’s equations. The fluid equations consist of the conservation of momentum and number density. A Darwin approximation of Maxwell is used to eliminate light waves from the model in order to facilitate coupling to non-relativistic particle models. We analyze the TFP-Darwin system in the context of a stand-alone solver with consideration of preconditioning a kinetic-JFNK approach. The TFP-Darwin system is addressed numerically by use of nested iteration (NI) and a First-Order Systems Least Squares (FOSLS) discretization. An important goal of NI is to produce an approximation that is within the basis of attraction for Newton’s method on a relatively coarse mesh and, thus, on all subsequent meshes. After scaling and modification, the TFP-Darwin model yields a nonlinear, first-order system of equations whose Fréchet derivative is shown to be uniformly H1-elliptic in a neighborhood of the exact solution. H1 ellipticity yields optimal finite element performance and linear systems amenable to solution with Algebraic Multigrid (AMG). To efficiently focus computational resources, an adaptive mesh refinement scheme, based on the accuracy per computational cost, is leveraged. Numerical tests demonstrate the efficacy of the approach, yielding an approximate solution within discretization error in a relatively small number of computational work units.