Chaotic magnetic field reversals in turbulent dynamos

We present direct numerical simulations of reversals of the magnetic field generated by swirling flows in a spherical domain. In agreement with a recent model, we observe that coupling dipolar and quadrupolar magnetic modes by an asymmetric forcing of the flow generates field reversals. In addition, we show that this mechanism strongly depends on the value of the magnetic Prandtl number.Comment: 4 pages, 5 figure

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed

Mode analysis of numerical geodynamo models

It has been suggested in Hoyng (2009) that dynamo action can be analysed by expansion of the magnetic field into dynamo modes and statistical evaluation of the mode coefficients. We here validate this method by analysing a numerical geodynamo model and comparing the numerically derived mean mode coefficients with the theoretical predictions. The model belongs to the class of kinematically stable dynamos with a dominating axisymmetric, antisymmetric with respect to the equator and non-periodic fundamental dynamo mode. The analysis requires a number of steps: the computation of the so-called dynamo coefficients, the derivation of the temporally and azimuthally averaged dynamo eigenmodes and the decomposition of the magnetic field of the numerical geodynamo model into the eigenmodes. For the determination of the theoretical mode excitation levels the turbulent velocity field needs to be projected on the dynamo eigenmodes. We compare the theoretically and numerically derived mean mode coefficients and find reasonably good agreement for most of the modes. Some deviation might be attributable to the approximation involved in the theory. Since the dynamo eigenmodes are not self-adjoint a spectral interpretation of the eigenmodes is not possible

A spherical shell numerical dynamo benchmark with pseudo vacuum magnetic boundary conditions

It is frequently considered that many planetary magnetic fields originate as a result of convection within planetary cores. Buoyancy forces responsible for driving the convection generate a fluid flow that is able to induce magnetic fields; numerous sophisticated computer codes are able to simulate the dynamic behaviour of such systems. This paper reports the results of a community activity aimed at comparing numerical results of several different types of computer codes that are capable of solving the equations of momentum transfer, magnetic field generation and heat transfer in the setting of a spherical shell, namely a sphere containing an inner core. The electrically conducting fluid is incompressible and rapidly rotating and the forcing of the flow is thermal convection under the Boussinesq approximation. We follow the original specifications and results reported in Harder & Hansen to construct a specific benchmark in which the boundaries of the fluid are taken to be impenetrable, non-slip and isothermal, with the added boundary condition for the magnetic field <b>B</b> that the field must be entirely radial there; this type of boundary condition for <b>B</b> is frequently referred to as ‘pseudo-vacuum’. This latter condition should be compared with the more frequently used insulating boundary condition. This benchmark is so-defined in order that computer codes based on local methods, such as finite element, finite volume or finite differences, can handle the boundary condition with ease. The defined benchmark, governed by specific choices of the Roberts, magnetic Rossby, Rayleigh and Ekman numbers, possesses a simple solution that is steady in an azimuthally drifting frame of reference, thus allowing easy comparison among results. Results from a variety of types of code are reported, including codes that are fully spectral (based on spherical harmonic expansions in angular coordinates and polynomial expansions in radius), mixed spectral and finite difference, finite volume, finite element and also a mixed Fourier-finite element code. There is good agreement among codes

Simulations of galactic dynamos

We review our current understanding of galactic dynamo theory, paying particular attention to numerical simulations both of the mean-field equations and the original three-dimensional equations relevant to describing the magnetic field evolution for a turbulent flow. We emphasize the theoretical difficulties in explaining non-axisymmetric magnetic fields in galaxies and discuss the observational basis for such results in terms of rotation measure analysis. Next, we discuss nonlinear theory, the role of magnetic helicity conservation and magnetic helicity fluxes. This leads to the possibility that galactic magnetic fields may be bi-helical, with opposite signs of helicity and large and small length scales. We discuss their observational signatures and close by discussing the possibilities of explaining the origin of primordial magnetic fields.Comment: 28 pages, 15 figure, to appear in Lecture Notes in Physics "Magnetic fields in diffuse media", Eds. E. de Gouveia Dal Pino and A. Lazaria

Continuum Halos in Nearby Galaxies -- an EVLA Survey (CHANG-ES) -- I: Introduction to the Survey

We introduce a new survey to map the radio continuum halos of a sample of 35 edge-on spiral galaxies at 1.5 GHz and 6 GHz in all polarization products. The survey is exploiting the new wide bandwidth capabilities of the Karl G. Jansky Very Large Array (i.e. the Expanded Very Large Array, or EVLA) in a variety of array configurations (B, C, and D) in order to compile the most comprehensive data set yet obtained for the study of radio halo properties. This is the first survey of radio halos to include all polarization products. In this first paper, we outline the scientific motivation of the survey, the specific science goals, and the expected improvements in noise levels and spatial coverage from the survey. Our goals include investigating the physical conditions and origin of halos, characterizing cosmic ray transport and wind speed, measuring Faraday rotation and mapping the magnetic field, probing the in-disk and extraplanar far-infrared - radio continuum relation, and reconciling non-thermal radio emission with high-energy gamma-ray models. The sample size allows us to search for correlations between radio halos and other properties, including environment, star formation rate, and the presence of AGNs. In a companion paper (Paper II) we outline the data reduction steps and present the first results of the survey for the galaxy, NGC 4631.Comment: 17 pages, 1 figure, accepted to the Astronomical Journal, Version 2 changes: added acknowledgement to NRA

The influence of swelling capacity of superdisintegrants in different pH media on the dissolution of hydrochlorothiazide from directly compressed tablets

The purpose of this study was to investigate the efficiency of superdisintegrants in promoting tablet disintegration and drug dissolution under varied media pH. Significant reductions in the rate and extent of water uptake and swelling were observed for both sodium starch glycolate (Primojel) and croscarmellose sodium (Ac-Di-Sol) in an acidic medium (0.1 N HCl) but not for crospovidone NF (Polyplasdone XL10), a nonionic polymer. When Primojel and Ac-Di-Sol were incorporated in model formulations, a significant increase in tablet disintegration time was observed for slowly disintegrating tablets (lactose-based tablets) but not for the rapidly disintegrating tablets (dicalcium phosphate-based tablets). The dissolution rate of the model drug, hydrochlorothiazide, was found highly dependent on both tablet disintegration efficiency and the solubility of base material(s) in the testing medium. A laser diffraction particle size analyzer proved to be an effective tool for determining the intrinsic swelling of disintegrant particles in different media. Water uptake and swelling were confirmed as 2 important functions of superdisintegrants. The reduced water uptake and swelling capacity of disintegrants containing ionizable substituents in an acidic medium can potentially jeopardize their efficiency in promoting tablet disintegration and the drug dissolution rate