Global MHD Simulations of Neptune's Magnetosphere

A global magnetohydrodynamic (MHD) simulation has been performed in order to investigate the outer boundaries of Neptune's magnetosphere at the time of Voyager 2's flyby in 1989 and to better understand the dynamics of magnetospheres formed by highly inclined planetary dipoles. Using the MHD code Gorgon, we have implemented a precessing dipole to mimic Neptune's tilted magnetic field and rotation axes. By using the solar wind parameters measured by Voyager 2, the simulation is verified by finding good agreement with Voyager 2 magnetometer observations. Overall, there is a large-scale reconfiguration of magnetic topology and plasma distribution. During the “pole-on” magnetospheric configuration, there only exists one tail current sheet, contained between a rarefied lobe region which extends outward from the dayside cusp, and a lobe region attached to the nightside cusp. It is found that the tail current always closes to the magnetopause current system, rather than closing in on itself, as suggested by other models. The bow shock position and shape is found to be dependent on Neptune's daily rotation, with maximum standoff being during the pole-on case. Reconnection is found on the magnetopause but is highly modulated by the interplanetary magnetic field (IMF) and time of day, turning “off” and “on” when the magnetic shear between the IMF and planetary fields is large enough. The simulation shows that the most likely location for reconnection to occur during Voyager 2's flyby was far from the spacecraft trajectory, which may explain the relative lack of associated signatures in the observations

Self-similar solutions for resistive diffusion, Ohmic heating, and Ettingshausen effects in plasmas of arbitrary beta

ABSTRACT Magneto-inertial fusion (MIF) approaches, such as the MagLIF experiment, use magnetic fields in dense plasma to suppress cross-field thermal conduction, attempting to reduce heat loss and trap alpha particles to achieve ignition. However, the magnetic field can introduce other transport effects, some of which are deleterious. An understanding of these processes is thus crucial for accurate modeling of MIF. We generalize past work exploiting self-similar solutions to describe transport processes in planar geometry and compare the model to the radiation-magnetohydrodynamics (MHDs) code Chimera. We solve the 1D extended MHD equations under pressure balance, making no assumptions about the ratio of magnetic and thermal pressures in the plasma. The resulting ordinary differential equation (ODE) boundary value problem is solved using a shooting method, combining an implicit ODE solver and a Newton–Raphson root finder. We show that the Nernst effect dominates over resistive diffusion in high β plasma, but its significance is reduced as the β decreases. On the other hand, we find that Ettingshausen and Ohmic heating effects are dominant in low β plasma and can be observable in even order unity β plasma, though in the presence of a strong temperature gradient heat conduction remains dominant. We then present a test problem for the Ohmic heating and Ettingshausen effects which will be useful to validate codes modeling these effects. We also observe that the Ettingshausen effect plays a role in preventing temperature separation when Ohmic heating is strong. Neglecting this term may lead to overestimates for the electron temperature at a vacuum–plasma interface, such as at the edge of a z-pinch. The model developed can be used to provide test problems with arbitrary boundary conditions for magnetohydrodynamics codes with the ability to freely switch on terms to compare their individual implementations

Dipole tilt effect on magnetopause reconnection and the steady‐state magnetosphere‐ionosphere system: global MHD simulation

The Earth’s dipole tilt angle changes both diurnally and seasonally and introduces numerous variabilities in the coupled magnetosphere‐ionosphere system. By altering the location and intensity of magnetic reconnection, the dipole tilt influences convection on a global scale. However, due to the nonlinear nature of the system, various other effects like dipole rotation, varying IMF orientation and non‐uniform ionospheric conductance can smear tilt effects arising purely from changes in coupling with the solar wind. To elucidate the underlying tilt angle‐dependence, we perform MHD simulations of the steady‐state magnetosphere‐ionosphere system under purely southward IMF conditions for tilt angles from 0°‐90°. We identify the location of the magnetic separator in each case, and find that an increasing tilt angle shifts the 3‐D X‐line southward on the magnetopause due to changes in magnetic shear angle. The separator is highly unsteady above 50° tilt angle, characteristic of regular FTE generation on the magnetopause. The reconnection rate drops as the tilt angle becomes large, but remains continuous across the dayside such that the magnetosphere is open even for 90°. These trends map down to the ionosphere, with the polar cap contracting as the tilt angle increases, and region‐I field‐aligned current (FAC) migrating to higher latitudes with changing morphology. The tilt introduces a north‐south asymmetry in magnetospheric convection, thus driving more FAC in the northern (sunward‐facing) hemisphere for large tilt angles than in the south independent of conductance. These results highlight the strong sensitivity to onset time in the potential impact of a severe space weather event

Synthetic nuclear diagnostics for inferring plasma properties of inertial confinement fusion implosions

A suite of synthetic nuclear diagnostics has been developed to post-process radiation hydrodynamics simulations performed with the code Chimera. These provide experimental observables based on simulated capsule properties and are used to assess alternative experimental and data analysis techniques. These diagnostics include neutron spectroscopy, primary and scattered neutron imaging, neutron activation, γ-ray time histories and carbon γ-ray imaging. Novel features of the neutron spectrum have been analysed to infer plasma parameters. The nT and nD backscatter edges have been shown to provide a shell velocity measurement. Areal density asymmetries created by low mode perturbations have been inferred from the slope of the downscatter spectrum down to 10 MeV. Neutron activation diagnostics showed significant aliasing of high mode areal density asymmetries when observing a capsule implosion with 3D multimode perturbations applied. Carbon γ-ray imaging could be used to image the ablator at a high convergence ratio. Time histories of both the fusion and carbon γ signals showed a greater time difference between peak intensities for the perturbed case when compared to a symmetric simulation

Neutron backscatter edge: A measure of the hydrodynamic properties of the dense DT fuel at stagnation in ICF experiments

The kinematic lower bound for the single scattering of neutrons produced in deuterium-tritium (DT) fusion reactions produces a backscatter edge in the measured neutron spectrum. The energy spectrum of backscattered neutrons is dependent on the scattering ion velocity distribution. As the neutrons preferentially scatter in the densest regions of the capsule, the neutron backscatter edge presents a unique measurement of the hydrodynamic conditions in the dense DT fuel. It is shown that the spectral shape of the edge is determined by the scattering rate weighted fluid velocity and temperature of the dense DT fuel layer during neutron production. In order to fit the neutron spectrum, a model for the various backgrounds around the backscatter edge is developed and tested on synthetic data produced from hydrodynamic simulations of OMEGA implosions. It is determined that the analysis could be utilized on current inertial confinement fusion experiments in order to measure the dense fuel properties

SpK: A fast atomic and microphysics code for the high-energy-density regime

SpK is part of the numerical codebase at Imperial College London used to model high energy density physics (HEDP) experiments. SpK is an efficient atomic and microphysics code used to perform detailed configuration accounting calculations of electronic and ionic stage populations, opacities and emissivities for use in post-processing and radiation hydrodynamics simulations. This is done using screened hydrogenic atomic data supplemented by the NIST energy level database. An extended Saha model solves for chemical equilibrium with extensions for non-ideal physics, such as ionisation potential depression, and non thermal equilibrium corrections. A tree-heap (treap) data structure is used to store spectral data, such as opacity, which is dynamic thus allowing easy insertion of points around spectral lines without a-priori knowledge of the ion stage populations. Results from SpK are compared to other codes and descriptions of radiation transport solutions which use SpK data are given. The treap data structure and SpK’s computational efficiency allows inline post-processing of 3D hydrodynamics simulations with a dynamically evolving spectrum stored in a treap

The structure of 3-D collisional magnetized bow shocks in pulsed-power-driven plasma flows

We investigate three-dimensional (3-D) bow shocks in a highly collisional magnetized aluminium plasma, generated during the ablation phase of an exploding wire array on the MAGPIE facility (1.4 MA, 240 ns). Ablation of plasma from the wire array generates radially diverging, supersonic ( MS∼7 ), super-Alfvénic ( MA>1 ) magnetized flows with frozen-in magnetic flux ( RM≫1 ). These flows collide with an inductive probe placed in the flow, which serves both as the obstacle that generates the magnetized bow shock, and as a diagnostic of the advected magnetic field. Laser interferometry along two orthogonal lines of sight is used to measure the line-integrated electron density. A detached bow shock forms ahead of the probe, with a larger opening angle in the plane parallel to the magnetic field than in the plane normal to it. Since the resistive diffusion length of the plasma is comparable to the probe size, the magnetic field decouples from the ion fluid at the shock front and generates a hydrodynamic shock, whose structure is determined by the sonic Mach number, rather than the magnetosonic Mach number of the flow. The 3-D simulations performed using the resistive magnetohydrodynamic (MHD) code Gorgon confirm this picture, but under-predict the anisotropy observed in the shape of the experimental bow shock, suggesting that non-MHD mechanisms may be important for modifying the shock structure

Impact of asymmetries on fuel performance in inertial confinement fusion

Low-mode asymmetries prevent effective compression, confinement, and heating of the fuel in inertial confinement fusion (ICF) implosions, and their control is essential to achieving ignition. Ion temperatures (Tion) in ICF experiments are inferred from the broadening of primary neutron spectra. Directional motion (flow) of the fuel at burn also impacts broadening and will lead to artificially inflated "Tion" values. Flow due to low-mode asymmetries is expected to give rise to line-of-sight variations in measured Tion. We report on intentionally asymmetrically driven experiments at the OMEGA laser facility designed to test the ability to accurately predict and measure line-of-sight differences in apparent Tion due to low-mode asymmetry-seeded flows. Contrasted to chimera and xrage simulations, the measurements demonstrate how all asymmetry seeds have to be considered to fully capture the flow field in an implosion. In particular, flow induced by the stalk that holds the target is found to interfere with the seeded asymmetry. A substantial stalk-seeded asymmetry in the areal density of the implosion is also observed