Researcher profile

M. S. Rosin

M. S. Rosin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

physics.plasm-ph astro-ph.CO astro-ph.GA astro-ph.HE physics.class-ph physics.comp-ph physics.flu-dyn physics.space-ph math-ph math.MP math.NA nlin.SI

Trust snapshot

Quick read

Trust 21 - Emerging

6works

0followers

12topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2013arXiv

An entropy based thermalization scheme for hybrid simulations of Coulomb collisions

We formulate and test a hybrid fluid-Monte Carlo scheme for the treatment of elastic collisions in gases and plasmas. While our primary focus and demonstrations of applicability are for moderately collisional plasmas, as described by the Landau-Fokker-Planck equation, the method is expected to be applicable also to collision processes described by the Boltzmann equation. This scheme is similar to the previously discussed velocity-based scheme [R. Caflisch, et. al, Multiscale Modeling & Simulation 7, 865, (2008)] and the scattering-angle-based scheme [A.M. Dimits, et. al, Bull. APS 55, no. 15 (2010, Abstract: XP9.00006)], but with a firmer theoretical basis and without the inherent limitation to the Landau-Fokker-Planck case. It gives a significant performance improvement (e.g., error for a given computational effort) over the velocity-based scheme. These features are achieved by assigning passive scalars to each simulated particle and tracking their evolution through collisions. The method permits a detailed error analysis that is confirmed by numerical results. The tests performed are for the evolution from anisotropic Maxwellian and a bump-on-tail distribution.

preprint2013arXiv

Multilevel Monte Carlo simulation of Coulomb collisions

We present a new, for plasma physics, highly efficient multilevel Monte Carlo numerical method for simulating Coulomb collisions. The method separates and optimally minimizes the finite-timestep and finite-sampling errors inherent in the Langevin representation of the Landau-Fokker-Planck equation. It does so by combining multiple solutions to the underlying equations with varying numbers of timesteps. For a desired level of accuracy epsilon, the computational cost of the method is order(epsilon^{-2}) or order(epsilon^{-2} (\ln epsilon)^2), depending on the underlying discretization, Milstein or Euler-Maruyama respectively. This is to be contrasted with a cost of order(epsilon^{-3}) for direct simulation Monte Carlo or binary collision methods. We successfully demonstrate the method with a classic beam diffusion test case in 2D, making use of the Levy area approximation for the correlated Milstein cross terms, and generating a computational saving of a factor of 100 for epsilon = 10^{-5}. We discuss the importance of the method for problems in which collisions constitute the computational rate limiting step, and its limitations.

preprint2012arXiv

Stability and energetics of Bursian diodes

We present an analysis of the stability, energy and torque properties of a model Bursian diode in a one dimensional Eulerian framework using the cold Euler-Poisson fluid equations. In regions of parameter space where there are two sets of equilibrium solutions for the same boundary conditions, one solution is found to be stable and the other unstable to linear perturbations. Following the linearly unstable solutions into the non-linear regime, we find they relax to the stable equilibrium. A description of this process in terms of kinetic, potential and boundary-flux energies is given, and the relation to a Hamiltonian formulation is commented upon. A non-local torque integral theorem, relating the prescribed boundary data to the average current in the domain, is also provided. These results should prove useful for understanding Bursian diodes in general, as well as for control applications and benchmarking numerical codes.

preprint2011arXiv

A nonlinear theory of the parallel firehose and gyrothermal instabilities in a weakly collisional plasma

Weakly collisional plasmas dynamically develop pressure anisotropies with respect to the magnetic field. These anisotropies trigger plasma instabilities at scales just above the ion Larmor radius ρ_i and much below the mean free path λ_{mfp}. They have growth rates of a fraction of the ion cyclotron frequency - much faster than either the global dynamics or local turbulence. The instabilities dramatically modify the transport properties and, therefore, the macroscopic dynamics of the plasma. Their nonlinear evolution drives pressure anisotropies towards marginal stability, controlled by the plasma beta β_i. Here this nonlinear evolution is worked out for the simplest analytically tractable example - the parallel firehose instability. In the nonlinear regime, both analytical theory and the numerical solution predict secular growth of magnetic fluctuations. They develop a k^{-3} spectrum, extending from scales somewhat larger than ρ_i to the maximum scale that grows secularly with time (~t^{1/2}); the relative pressure anisotropy (\pperp-\ppar)/\ppar tends to the marginal value -2/β_i. The marginal state is achieved via changes in the magnetic field, not particle scattering. When a parallel ion heat flux is present, the firehose mutates into the new gyrothermal instability (GTI), which continues to exist up to firehose-stable values of pressure anisotropy, which can be positive and are limited by the heat flux. The nonlinear evolution of the GTI also features secular growth of magnetic fluctuations, but the spectrum is eventually dominated by modes around the scale ~ρ_i l_T/λ_{mfp}, where l_T is the scale of the parallel temperature variation. Implications for momentum and heat transport are speculated about. This study is motivated by the dynamics of galaxy cluster plasmas.

preprint2011arXiv

Beyond the Child-Langmuir Limit

This paper describes a new solution formulation for fully nonlinear and unsteady planar flow of an electron beam in a diode. Using characteristic variables - i.e., variables that follow particle paths - the solution is expressed through an exact analytic, but implicit, formula for any choice of incoming velocity $v_0$, electric field $E_0$ and current $J_0$. For steady solutions, this approach clarifies the origin of the maximal current $J_max$, derived by Child and Langmuir for $v_0=0$ and by Jaffe for $v_0>0$. The implicit formulation is used to find (1) unsteady solutions having constant incoming flux $J_0>J_max$, which leads formation of a virtual cathode, and (2) time-periodic solutions whose average flux exceeds the adiabatic average of $J_max$.

preprint2010arXiv

Magnetofluid dynamics of magnetized cosmic plasma: firehose and gyrothermal instabilities

Both global dynamics and turbulence in magnetized weakly collisional cosmic plasmas are described by general magnetofluid equations that contain pressure anisotropies and heat fluxes that must be calculated from microscopic plasma kinetic theory. It is shown that even without a detailed calculation of the pressure anisotropy or the heat fluxes, one finds the macroscale dynamics to be generically unstable to microscale Alfvenically polarized fluctuations. Two instabilities are considered in detail: the parallel firehose instability (including the finite-Larmor-radius effects that determine the fastest growing mode) and the gyrothermal instability (GTI). The latter is a new result - it is shown that a parallel ion heat flux destabilizes Alfvenically polarized fluctuations even in the absence of the negative pressure anisotropy required for the firehose. The main conclusion is that both pressure anisotropies and heat fluxes trigger plasma microinstabilities and, therefore, their values will likely be set by the nonlinear evolution of these instabilities. Ideas for understanding this nonlinear evolution are discussed. It is argued that cosmic plasmas will generically be "three-scale systems," comprising global dynamics, mesoscale turbulence and microscale plasma fluctuations. The astrophysical example of cool cores of galaxy clusters is considered and it is noted that observations point to turbulence in clusters being in a marginal state with respect to plasma microinstabilities and so it is the plasma microphysics that is likely to set the heating and conduction properties of the intracluster medium. In particular, a lower bound on the scale of temperature fluctuations implied by the GTI is derived.