Source author record

Andrei N. Simakov

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

physics.plasm-ph physics.space-ph astro-ph.SR physics.comp-ph

Catalog footprint

Research overview

4works

4topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2020arXiv

An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system

We consider the solution of the fully kinetic (including electrons) Vlasov-Ampère system in a one-dimensional physical space and two-dimensional velocity space (1D-2V) for an arbitrary number of species with a time-implicit Eulerian algorithm. The problem of velocity-space meshing for disparate thermal and bulk velocities is dealt with by an adaptive coordinate transformation of the Vlasov equation for each species, which is then discretized, including the resulting inertial terms. Mass, momentum, and energy are conserved, and Gauss's law is enforced to within the nonlinear convergence tolerance of the iterative solver through a set of nonlinear constraint functions while permitting significant flexibility in choosing discretizations in time, configuration, and velocity space. We mitigate the temporal stiffness introduced by, e.g., the plasma frequency through the use of high-order/low-order (HOLO) acceleration of the iterative implicit solver. We present several numerical results for canonical problems of varying degrees of complexity, including the multiscale ion-acoustic shock wave problem, which demonstrate the efficacy, accuracy, and efficiency of the scheme.

preprint2015arXiv

Fluid vs. kinetic magnetic reconnection with strong guide-fields

The fast rates of magnetic reconnection found in both nature and experiments are important to understand theoretically. Recently, it was demonstrated that two-fluid magnetic reconnection remains fast in the strong guide field regime, regardless of the presence of fast-dispersive waves. This conclusion is in agreement with recent results from kinetic simulations, and is in contradiction to the findings in an earlier two-fluid study, where it was suggested that fast-dispersive waves are necessary for fast reconnection. In this paper, we give a more detailed derivation of the analytic model presented in a recent letter, and present additional simulation results to support the conclusions that the magnetic reconnection rate in this regime is independent of both collisional dissipation and system-size. In particular, we present a detailed comparison between fluid and kinetic simulations, finding good agreement in both the reconnection rate and overall length of the current layer. Finally, we revisit the earlier two-fluid study, which arrived at different conclusions, and suggest an alternative interpretation for the numerical results presented therein.

preprint2010arXiv

Analytical Model of Fast Magnetic Reconnection with a Large Guide Field

Analytical theory of fast magnetic reconnection with a large guide field is presented for the first time. We confirm that two distinct reconnection regimes are possible depending on whether the diffusion region thickness $δ$ is larger or smaller than the sound gyroradius $ρ_s$ . The reconnection is slow or Sweet-Parker-like for $δ\gtrsim ρ_s$ , and fast otherwise. In the fast regime, however, we find that ion viscosity $μ$ plays a critical role. In particular, for $δ< ρ_s$ the diffusion region thickness is proportional to $Ha^{-1}$ with $Ha \propto 1/\sqrt{ημ}$ the Hartmann number, and the reconnection rate is proportional to $Pr^{-1/2}$ with $Pr = μ/η$ the Prandtl number and $η$ the resistivity. If the perpendicular ion viscosity is employed for $μ$ the reconnection rate becomes independent of both plasma $β$ and collision frequencies and therefore potentially fast.

preprint2009arXiv

Current sheet bifurcation and collapse in electron magnetohydrodynamics

Inertial effects in nonlinear magnetic reconnection are studied within the context of 2D electron magnetohydrodynamics (EMHD) with resistive and viscous dissipation. Families of nonlinear solutions for relevant current sheet parameters are predicted and confirmed numerically in all regimes of interest. Electron inertia becomes important for current sheet thicknesses $δ$ below the inertial length $d_{e}$ . In this case, in the absence of electron viscosity, the sheet thickness experiences a nonlinear collapse. Viscosity regularizes solutions at small scales. Transition from resistive to viscous regimes shows a nontrivial dependence on resistivity and viscosity, featuring a hysteresis bifurcation. In all accessible regimes, the nonlinear reconnection rate is found to be explicitly independent of the electron inertia and dissipation coefficients.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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