ASHEE: a compressible, equilibrium-Eulerian model for volcanic ash plumes

A new fluid-dynamic model is developed to numerically simulate the non-equilibrium dynamics of polydisperse gas-particle mixtures forming volcanic plumes. Starting from the three-dimensional N-phase Eulerian transport equations for a mixture of gases and solid particles, we adopt an asymptotic expansion strategy to derive a compressible version of the first-order non-equilibrium model, valid for low concentration regimes and small particles Stokes $St<0.2$. When $St < 0.001$ the model reduces to the dusty-gas one. The new model is significantly faster than the Eulerian model while retaining the capability to describe gas-particle non-equilibrium. Direct numerical simulation accurately reproduce the dynamics of isotropic turbulence in subsonic regime. For gas-particle mixtures, it describes the main features of density fluctuations and the preferential concentration of particles by turbulence, verifying the model reliability and suitability for the simulation of high-Reynolds number and high-temperature regimes. On the other hand, Large-Eddy Numerical Simulations of forced plumes are able to reproduce their observed averaged and instantaneous properties. The self-similar radial profile and the development of large-scale structures are reproduced, including the rate of entrainment of atmospheric air. Application to the Large-Eddy Simulation of the injection of the eruptive mixture in a stratified atmosphere describes some of important features of turbulent volcanic plumes, including air entrainment, buoyancy reversal, and maximum plume height. Coarse particles partially decouple from the gas within eddies, modifying the turbulent structure, and preferentially concentrate at the eddy periphery, eventually being lost from the plume margins due to the gravity. By these mechanisms, gas-particle non-equilibrium is able to influence the large-scale behavior of volcanic plumes.Comment: 29 pages, 22 figure

Electron-scale shear instabilities: magnetic field generation and particle acceleration in astrophysical jets

Strong shear flow regions found in astrophysical jets are shown to be important dissipation regions, where the shear flow kinetic energy is converted into electric and magnetic field energy via shear instabilities. The emergence of these self-consistent fields make shear flows significant sites for radiation emission and particle acceleration. We focus on electron-scale instabilities, namely the collisionless, unmagnetized Kelvin-Helmholtz instability (KHI) and a large-scale dc magnetic field generation mechanism on the electron scales. We show that these processes are important candidates to generate magnetic fields in the presence of strong velocity shears, which may naturally originate in energetic matter outburst of active galactic nuclei and gamma-ray bursters. We show that the KHI is robust to density jumps between shearing flows, thus operating in various scenarios with different density contrasts. Multidimensional particle-in-cell (PIC) simulations of the KHI, performed with OSIRIS, reveal the emergence of a strong and large-scale dc magnetic field component, which is not captured by the standard linear fluid theory. This dc component arises from kinetic effects associated with the thermal expansion of electrons of one flow into the other across the shear layer, whilst ions remain unperturbed due to their inertia. The electron expansion forms dc current sheets, which induce a dc magnetic field. Our results indicate that most of the electromagnetic energy developed in the KHI is stored in the dc component, reaching values of equipartition on the order of $10^{-3}$ in the electron time-scale, and persists longer than the proton time-scale. Particle scattering/acceleration in the self generated fields of these shear flow instabilities is also analyzed

New numerical approaches for modeling thermochemical convection in a compositionally stratified fluid

Seismic imaging of the mantle has revealed large and small scale heterogeneities in the lower mantle; specifically structures known as large low shear velocity provinces (LLSVP) below Africa and the South Pacific. Most interpretations propose that the heterogeneities are compositional in nature, differing in composition from the overlying mantle, an interpretation that would be consistent with chemical geodynamic models. Numerical modeling of persistent compositional interfaces presents challenges, even to state-of-the-art numerical methodology. For example, some numerical algorithms for advecting the compositional interface cannot maintain a sharp compositional boundary as the fluid migrates and distorts with time dependent fingering due to the numerical diffusion that has been added in order to maintain the upper and lower bounds on the composition variable and the stability of the advection method. In this work we present two new algorithms for maintaining a sharper computational boundary than the advection methods that are currently openly available to the computational mantle convection community; namely, a Discontinuous Galerkin method with a Bound Preserving limiter and a Volume-of-Fluid interface tracking algorithm. We compare these two new methods with two approaches commonly used for modeling the advection of two distinct, thermally driven, compositional fields in mantle convection problems; namely, an approach based on a high-order accurate finite element method advection algorithm that employs an artificial viscosity technique to maintain the upper and lower bounds on the composition variable as well as the stability of the advection algorithm and the advection of particles that carry a scalar quantity representing the location of each compositional field. All four of these algorithms are implemented in the open source FEM code ASPECT

A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows IV: full Boltzmann and Model Equations

Fluid dynamic equations are valid in their respective modeling scales. With a variation of the modeling scales, theoretically there should have a continuous spectrum of fluid dynamic equations. In order to study multiscale flow evolution efficiently, the dynamics in the computational fluid has to be changed with the scales. A direct modeling of flow physics with a changeable scale may become an appropriate approach. The unified gas-kinetic scheme (UGKS) is a direct modeling method in the mesh size scale, and its underlying flow physics depends on the resolution of the cell size relative to the particle mean free path. The cell size of UGKS is not limited by the particle mean free path. With the variation of the ratio between the numerical cell size and local particle mean free path, the UGKS recovers the flow dynamics from the particle transport and collision in the kinetic scale to the wave propagation in the hydrodynamic scale. The previous UGKS is mostly constructed from the evolution solution of kinetic model equations. This work is about the further development of the UGKS with the implementation of the full Boltzmann collision term in the region where it is needed. The central ingredient of the UGKS is the coupled treatment of particle transport and collision in the flux evaluation across a cell interface, where a continuous flow dynamics from kinetic to hydrodynamic scales is modeled. The newly developed UGKS has the asymptotic preserving (AP) property of recovering the NS solutions in the continuum flow regime, and the full Boltzmann solution in the rarefied regime. In the mostly unexplored transition regime, the UGKS itself provides a valuable tool for the flow study in this regime. The mathematical properties of the scheme, such as stability, accuracy, and the asymptotic preserving, will be analyzed in this paper as well