Magnetohydrodynamic Modeling of Three Van Allen Probes Storms in 2012 and 2013

Coronal mass ejection (CME)-shock compression of the dayside magnetopause has been observed to cause both prompt enhancement of radiation belt electron flux due to inward radial transport of electrons conserving their first adiabatic invariant and prompt losses which at times entirely eliminate the outer zone. Recent numerical studies suggest that enhanced ultra-low frequency (ULF) wave activity is necessary to explain electron losses deeper inside the magnetosphere than magnetopause incursion following CME-shock arrival. A combination of radial transport and magnetopause shadowing can account for losses observed at radial distances into L=4.5, well within the computed magnetopause location. We compare ULF wave power from the Electric Field and Waves (EFW) electric field instrument on the Van Allen Probes for the 8 October 2013 storm with ULF wave power simulated using the Lyon–Fedder–Mobarry (LFM) global magnetohydrodynamic (MHD) magnetospheric simulation code coupled to the Rice Convection Model (RCM). Two other storms with strong magnetopause compression, 8–9 October 2012 and 17–18 March 2013, are also examined. We show that the global MHD model captures the azimuthal magnetosonic impulse propagation speed and amplitude observed by the Van Allen Probes which is responsible for prompt acceleration at MeV energies reported for the 8 October 2013 storm. The simulation also captures the ULF wave power in the azimuthal component of the electric field, responsible for acceleration and radial transport of electrons, at frequencies comparable to the electron drift period. This electric field impulse has been shown to explain observations in related studies (Foster et al., 2015) of electron acceleration and drift phase bunching by the Energetic Particle, Composition, and Thermal Plasma Suite (ECT) instrument on the Van Allen Probes

Magnetohydrodynamic activity inside a sphere

We present a computational method to solve the magnetohydrodynamic equations in spherical geometry. The technique is fully nonlinear and wholly spectral, and uses an expansion basis that is adapted to the geometry: Chandrasekhar-Kendall vector eigenfunctions of the curl. The resulting lower spatial resolution is somewhat offset by being able to build all the boundary conditions into each of the orthogonal expansion functions and by the disappearance of any difficulties caused by singularities at the center of the sphere. The results reported here are for mechanically and magnetically isolated spheres, although different boundary conditions could be studied by adapting the same method. The intent is to be able to study the nonlinear dynamical evolution of those aspects that are peculiar to the spherical geometry at only moderate Reynolds numbers. The code is parallelized, and will preserve to high accuracy the ideal magnetohydrodynamic (MHD) invariants of the system (global energy, magnetic helicity, cross helicity). Examples of results for selective decay and mechanically-driven dynamo simulations are discussed. In the dynamo cases, spontaneous flips of the dipole orientation are observed.Comment: 15 pages, 19 figures. Improved figures, in press in Physics of Fluid

Orthogonal, solenoidal, three-dimensional vector fields for no-slip boundary conditions

Viscous fluid dynamical calculations require no-slip boundary conditions. Numerical calculations of turbulence, as well as theoretical turbulence closure techniques, often depend upon a spectral decomposition of the flow fields. However, such calculations have been limited to two-dimensional situations. Here we present a method that yields orthogonal decompositions of incompressible, three-dimensional flow fields and apply it to periodic cylindrical and spherical no-slip boundaries.Comment: 16 pages, 2 three-part figure

Local estimates for entropy densities in coupled map lattices

We present a method to derive an upper bound for the entropy density of coupled map lattices with local interactions from local observations. To do this, we use an embedding technique being a combination of time delay and spatial embedding. This embedding allows us to identify the local character of the equations of motion. Based on this method we present an approximate estimate of the entropy density by the correlation integral.Comment: 4 pages, 5 figures include

Hydrodynamic and magnetohydrodynamic computations inside a rotating sphere

Numerical solutions of the incompressible magnetohydrodynamic (MHD) equations are reported for the interior of a rotating, perfectly-conducting, rigid spherical shell that is insulator-coated on the inside. A previously-reported spectral method is used which relies on a Galerkin expansion in Chandrasekhar-Kendall vector eigenfunctions of the curl. The new ingredient in this set of computations is the rigid rotation of the sphere. After a few purely hydrodynamic examples are sampled (spin down, Ekman pumping, inertial waves), attention is focused on selective decay and the MHD dynamo problem. In dynamo runs, prescribed mechanical forcing excites a persistent velocity field, usually turbulent at modest Reynolds numbers, which in turn amplifies a small seed magnetic field that is introduced. A wide variety of dynamo activity is observed, all at unit magnetic Prandtl number. The code lacks the resolution to probe high Reynolds numbers, but nevertheless interesting dynamo regimes turn out to be plentiful in those parts of parameter space in which the code is accurate. The key control parameters seem to be mechanical and magnetic Reynolds numbers, the Rossby and Ekman numbers (which in our computations are varied mostly by varying the rate of rotation of the sphere) and the amount of mechanical helicity injected. Magnetic energy levels and magnetic dipole behavior are exhibited which fluctuate strongly on a time scale of a few eddy turnover times. These seem to stabilize as the rotation rate is increased until the limit of the code resolution is reached.Comment: 26 pages, 17 figures, submitted to New Journal of Physic

Comparison of the HAARM-3 Fallout Model with Nuclear Safety Pilot Plant (NSPP) Data

Sodium--oxide aerosols have been generated in the Nuclear Safety Pilot Plant (NSPP) experiment under vessel conditions appropriate for the LMFBR secondary containment following an HCDA. The primary goal of the NSPP program is to provide experimental validation for the HAARM-3 aerosol behavioral code. Validation can be determined by direct comparison of code predictions and integral test results using the same initial and boundary conditions. In addition, separate physical models within the code (e.g., for agglomeration, plate-out, and gravitational fallout) may be validated by a comparison of predicted values with measured instantaneous values. A comparison is presented for the gravitational fallout rate which is considered to be the primary depletion mechanism under prototypic secondary containment conditions

Just plain Wronga?: a multimodal critical analysis of online payday loan discourse

Payday loans constitute one of the most rapidly expanding and controversial forms of consumer lending today. Payday lending – the selling of high-interest, short-term credit – has thrived in the wake of the decline of the traditional high street banking system and the reluctance on the part of many mainstream credit services, following the 2007/8 Global Financial Crisis, to lend to low income earners. This study critically examines the website of the industry leader in the UK, Wonga, a payday lender which recently rebranded and relaunched itself (in 2015) after being embroiled in a series of financial scandals. Our analysis centres on the new Wonga website, the gateway to its financial services, and identifies three inter-related discursive strategies through which the lender, in the wake of its financial misconduct, seeks to present itself as a reputable financial service provider, namely by (1) constructing the empowered and responsible borrower, (2) de-stigmatising both its service provision and its prospective customers, the payday borrower, and (3) minimising the consequences and risks associated with payday borrowing. We argue that, collectively, these strategies constitute an artful response by Wonga to the changing legislative and socio-economic contexts in which it and other payday lenders are now required to operate, permitting it to continue marketing and selling its high-interest rate financial services

Evolution of Ultracold, Neutral Plasmas

We present the first large-scale simulations of an ultracold, neutral plasma, produced by photoionization of laser-cooled xenon atoms, from creation to initial expansion, using classical molecular dynamics methods with open boundary conditions. We reproduce many of the experimental findings such as the trapping efficiency of electrons with increased ion number, a minimum electron temperature achieved on approach to the photoionization threshold, and recombination into Rydberg states of anomalously-low principal quantum number. In addition, many of these effects establish themselves very early in the plasma evolution ($\sim$ ns) before present experimental observations begin.Comment: 4 pages, 3 figures, submitted to PR

Electronic structure of the (111) and (-1-1-1) surfaces of cubic BN: A local-density-functional ab initio study

We present ab initio local-density-functional electronic structure calculations for the (111) and (-1-1-1) surfaces of cubic BN. The energetically stable reconstructions, namely the N adatom, N3 triangle models on the (111), the (2x1), boron and nitrogen triangle patterns on the (-1-1-1) surface are investigated. Band structure and properties of the surface states are discussed in detail.Comment: 8 pages, 12 figure