Nonlinear numerical studies of macroscopic modes in a variety of magnetic fusion experiments are made possible by the flexible high-order accurate spatial representation and semi-implicit time advance in the NIMROD simulation code ͓A. H. Glasser et al., Plasma Phys. Controlled Fusion 41, A747 ͑1999͔͒. Simulation of a resistive magnetohydrodynamics mode in a shaped toroidal tokamak equilibrium demonstrates computation with disparate time scales, simulations of discharge 87009 in the DIII-D tokamak ͓J. L. Luxon et al., Plasma Physics and Controlled Nuclear Fusion Research 1986 ͑International Atomic Energy Agency, Vienna, 1987͒, Vol. I, p. 159͔ confirm an analytic scaling for the temporal evolution of an ideal mode subject to plasma- increasing beyond marginality, and a spherical torus simulation demonstrates nonlinear free-boundary capabilities. A comparison of numerical results on magnetic relaxation finds the nϭ1 mode and flux amplification in spheromaks to be very closely related to the mϭ1 dynamo modes and magnetic reversal in reversed-field pinch configurations. Advances in local and nonlocal closure relations developed for modeling kinetic effects in fluid simulation are also described.