Identification of the Sequence of Steps Intrinsic to Spheromak Formation

A planar coaxial electrostatic helicity source is used for studying the relaxation process intrinsic to spheromak formation Experimental observations reveal that spheromak formation involves: (1) breakdown and creation of a number of distinct, arched, filamentary, plasma-filled flux loops that span from cathode to anode gas nozzles, (2) merging of these loops to form a central column, (3) jet-like expansion of the central column, (4) kink instability of the central column, (5) conversion of toroidal flux to poloidal flux by the kink instability. Steps 1 and 3 indicate that spheromak formation involves an MHD pumping of plasma from the gas nozzles into the magnetic flux tube linking the nozzles. In order to measure this pumping, the gas puffing system has been modified to permit simultaneous injection of different gas species into the two ends of the flux tube linking the wall. Gated CCD cameras with narrow-band optical filters are used to track the pumped flows

Consideration of the relationship between Kepler and cyclotron dynamics leading to prediction of a non-MHD gravity-driven Hamiltonian dynamo

Conservation of canonical angular momentum shows that charged particles are typically constrained to stay within a poloidal Larmor radius of a poloidal magnetic flux surface. However, more detailed consideration shows that particles with a critical charge to mass ratio can have zero canonical angular momentum and so be both immune from centrifugal force and not constrained to stay in the vicinity of a specific flux surface. Suitably charged dust grains can have zero canonical angular momentum and in the presence of a gravitational field will spiral inwards across poloidal magnetic surfaces toward the central object and accumulate. This accumulation results in a gravitationally-driven dynamo, i.e., a mechanism for converting gravitational potential energy into a battery-like electric power source.Comment: 14 pages, 1 figur

Experimental Identification of the Kink Instability as a Poloidal Flux Amplification Mechanism for Coaxial Gun Spheromak Formation

The magnetohydrodynamic kink instability is observed and identified experimentally as a poloidal flux amplification mechanism for coaxial gun spheromak formation. Plasmas in this experiment fall into three distinct regimes which depend on the peak gun current to magnetic flux ratio, with (I) low values resulting in a straight plasma column with helical magnetic field, (II) intermediate values leading to kinking of the column axis, and (III) high values leading immediately to a detached plasma. Onset of column kinking agrees quantitatively with the Kruskal-Shafranov limit, and the kink acts as a dynamo which converts toroidal to poloidal flux. Regime II clearly leads to both poloidal flux amplification and the development of a spheromak configuration.Comment: accepted for publication in Physical Review Letter

Modelling of subgrid-scale phenomena in supercritical transitional mixing layers: an a priori study

A database of transitional direct numerical simulation (DNS) realizations of a supercritical mixing layer is analysed for understanding small-scale behaviour and examining subgrid-scale (SGS) models duplicating that behaviour. Initially, the mixing layer contains a single chemical species in each of the two streams, and a perturbation promotes roll-up and a double pairing of the four spanwise vortices initially present. The database encompasses three combinations of chemical species, several perturbation wavelengths and amplitudes, and several initial Reynolds numbers specifically chosen for the sole purpose of achieving transition. The DNS equations are the Navier-Stokes, total energy and species equations coupled to a real-gas equation of state; the fluxes of species and heat include the Soret and Dufour effects. The large-eddy simulation (LES) equations are derived from the DNS ones through filtering. Compared to the DNS equations, two types of additional terms are identified in the LES equations: SGS fluxes and other terms for which either assumptions or models are necessary. The magnitude of all terms in the LES conservation equations is analysed on the DNS database, with special attention to terms that could possibly be neglected. It is shown that in contrast to atmospheric-pressure gaseous flows, there are two new terms that must be modelled: one in each of the momentum and the energy equations. These new terms can be thought to result from the filtering of the nonlinear equation of state, and are associated with regions of high density-gradient magnitude both found in DNS and observed experimentally in fully turbulent high-pressure flows. A model is derived for the momentum-equation additional term that performs well at small filter size but deteriorates as the filter size increases, highlighting the necessity of ensuring appropriate grid resolution in LES. Modelling approaches for the energy-equation additional term are proposed, all of which may be too computationally intensive in LES. Several SGS flux models are tested on an a priori basis. The Smagorinsky (SM) model has a poor correlation with the data, while the gradient (GR) and scale-similarity (SS) models have high correlations. Calibrated model coefficients for the GR and SS models yield good agreement with the SGS fluxes, although statistically, the coefficients are not valid over all realizations. The GR model is also tested for the variances entering the calculation of the new terms in the momentum and energy equations; high correlations are obtained, although the calibrated coefficients are not statistically significant over the entire database at fixed filter size. As a manifestation of the small-scale supercritical mixing peculiarities, both scalar-dissipation visualizations and the scalar-dissipation probability density functions (PDF) are examined. The PDF is shown to exhibit minor peaks, with particular significance for those at larger scalar dissipation values than the mean, thus significantly departing from the Gaussian behaviour

A Model of Reduced Kinetics for Alkane Oxidation Using Constituents and Species for N-Heptane

The reduction of elementary or skeletal oxidation kinetics to a subgroup of tractable reactions for inclusion in turbulent combustion codes has been the subject of numerous studies. The skeletal mechanism is obtained from the elementary mechanism by removing from it reactions that are considered negligible for the intent of the specific study considered. As of now, there are many chemical reduction methodologies. A methodology for deriving a reduced kinetic mechanism for alkane oxidation is described and applied to n-heptane. The model is based on partitioning the species of the skeletal kinetic mechanism into lights, defined as those having a carbon number smaller than 3, and heavies, which are the complement of the species ensemble. For modeling purposes, the heavy species are mathematically decomposed into constituents, which are similar but not identical to groups in the group additivity theory. From analysis of the LLNL (Lawrence Livermore National Laboratory) skeletal mechanism in conjunction with CHEMKIN II, it is shown that a similarity variable can be formed such that the appropriately non-dimensionalized global constituent molar density exhibits a self-similar behavior over a very wide range of equivalence ratios, initial pressures and initial temperatures that is of interest for predicting n-heptane oxidation. Furthermore, the oxygen and water molar densities are shown to display a quasi-linear behavior with respect to the similarity variable. The light species ensemble is partitioned into quasi-steady and unsteady species. The reduced model is based on concepts consistent with those of Large Eddy Simulation (LES) in which functional forms are used to replace the small scales eliminated through filtering of the governing equations; in LES, these small scales are unimportant as far as the overwhelming part of dynamic energy is concerned. Here, the scales thought unimportant for recovering the thermodynamic energy are removed. The concept is tested by using tabular information from the LLNL skeletal mechanism in conjunction with CHEMKIN II utilized as surrogate ideal functions replacing the necessary functional forms. The test reveals that the similarity concept is indeed justified and that the combustion temperature is well predicted, but that the ignition time is over-predicted, a fact traced to neglecting a detailed description of the processes leading to the heavies chemical decomposition. To palliate this deficiency, functional modeling is incorporated into this conceptual reduction in addition to the modeling the evolution of the global constituent molar density, the enthalpy evolution of the heavies, the contribution to the reaction rate of the unsteady lights from other light species and from the heavies, the molar density evolution of oxygen and water, and the mole fractions of the quasisteady light species. The model is compact in that there are only nine species-related progress variables. Results are presented showing the performance of the model for predicting the temperature and species evolution. The model reproduces the ignition time over a wide range of equivalence ratios, initial pressure, and initial temperature

Mixing rules for multicomponent mixture mass diffusion coefficients and thermal diffusion factors

Mixing rules are derived for mass diffusion coefficient and thermal diffusion factor matrices by developing compatibility conditions between the fluid mixture equations obtained from nonequilibrium thermodynamics and Grad’s 13-moment kinetic theory. The mixing rules are shown to be in terms of the species mole fractions and binary processes. In particular, the thermal diffusion factors for binary mixtures obtained by the Chapman–Enskog expansion procedure are suitably generalized for many-component mixtures. Some practical aspects of the results are discussed including the utilization of these mixing rules for high pressure situations

Neutrino magnetohydrodynamics

A new neutrino magnetohydrodynamics (NMHD) model is formulated, where the effects of the charged weak current on the electron-ion magnetohydrodynamic fluid are taken into account. The model incorporates in a systematic way the role of the Fermi neutrino weak force in magnetized plasmas. A fast neutrino-driven short wavelengths instability associated with the magnetosonic wave is derived. Such an instability should play a central role in strongly magnetized plasma as occurs in supernovae, where dense neutrino beams also exist. In addition, in the case of nonlinear or high frequency waves, the neutrino coupling is shown to be responsible for breaking the frozen-in magnetic field lines condition even in infinite conductivity plasmas. Simplified and ideal NMHD assumptions were adopted and analyzed in detail

Laboratory simulations of astrophysical jets and solar coronal loops: new results

An experimental program underway at Caltech has produced plasmas where the shape is neither fixed by the vacuum chamber nor fixed by an external coil set, but instead is determined by self-organization. The plasma dynamics is highly reproducible and so can be studied in considerable detail even though the morphology of the plasma is both complex and time-dependent. A surprising result has been the observation that self-collimating MHD-driven plasma jets are ubiquitous and play a fundamental role in the self-organization. The jets can be considered lab-scale simulations of astrophysical jets and in addition are intimately related to solar coronal loops. The jets are driven by the combination of the axial component of the J×B force and the axial pressure gradient resulting from the non-uniform pinch force associated with the flared axial current density. Behavior is consistent with a model showing that collimation results from axial non-uniformity of the jet velocity. In particular, flow stagnation in the jet frame compresses frozen-in azimuthal magnetic flux, squeezes together toroidal magnetic field lines, thereby amplifying the embedded toroidal magnetic field, enhancing the pinch force, and hence causing collimation of the jet

Statistical Model of Multicomponent-Fuel Drop Evaporation for Many-Drop Flow Simulations

A statistical formulation is developed describing the composition in an evaporating multicomponent-fuel liquid drop and in the gas phase surrounding it. When a complementary discrete-component model is used, it is shown that, when drops are immersed in a carrier gas containing fuel vapor, condensation of species onto the drop results in the development of a minor peak in the liquid composition probability distribution function (PDF). This peak leads to a PDF shape that can be viewed as a combination of two gamma PDFs, which is determined by five parameters. A model is developed for calculating the parameters of the two combined gamma PDFs. Extensive tests of the model for both diesel and gasoline show that the PDF results replicate accurately the discrete model predictions. Most important, the mean and variance of the composition at the drop surface are in excellent agreement with the discrete model. Results from the model show that although the second peak is minor for the liquid PDF, its corresponding peak for the vapor distribution at the drop surface has a comparable magnitude to and sometimes exceeds that corresponding to the first peak. Four-parameter models are also exercised, and it is shown that they are unable to capture the physics of the problem