The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and study the turbulent dynamo mechanism in the laboratory. Here we describe the design of these experiments through simulation campaigns using FLASH, a highly capable radiation magnetohydrodynamics code that we have developed, and large-scale three-dimensional simulations on the Mira supercomputer at Argonne National Laboratory. The simulation results indicate that the experimental platform may be capable of reaching a turbulent plasma state and study dynamo amplification. We validate and compare our numerical results with a small subset of experimental data using synthetic diagnostics.

Self-gravity computation by multipole expansion is a common approach in problems such as core-collapse and Type Ia supernovae, where single large condensations of mass must be treated. The standard formulation of multipole self-gravity suffers from two significant sources of error, which we correct in the formulation presented in this article. The first source of error is due to the numerical approximation that effectively places grid cell mass at the central point of the cell, then computes the gravitational potential at that point, resulting in a convergence failure of the multipole expansion. We describe a new scheme that avoids this problem by computing gravitational potential at cell faces. The second source of error is due to sub-optimal choice of location for the expansion center, which results in angular power at high multipole $l$ values in the gravitational field, requiring a high --- and expensive --- value of multipole cutoff \lmax. By introducing a global measure of angular power in the gravitational field, we show that the optimal coordinate for the expansion is the square-density-weighted mean location. We subject our new multipole self-gravity algorithm to two rigorous test problems: MacLaurin spheroids for which exact analytic solutions are known, and core-collapse supernovae. We show that key observables of the core-collapse simulations, particularly shock expansion, proto-neutron star motion, and momentum conservation, are extremely sensitive to the accuracy of the multipole gravity, and the accuracy of their computation is greatly improved by our reformulated solver.