We present a three-dimensional study of the plasma dynamics at the flank magnetopause of the Earth's magnetosphere during mainly northward interplanetary magnetic field (IMF) periods. Two-fluid simulations show that the initial magnetic shear at the magnetopause and the field line bending caused by the dynamics itself (in a configuration taken as representative of the properties of the flank magnetopause) influence both the location where the Kelvin-Helmholtz (KH) instability and the induced magnetic reconnection take place and their nonlinear development. The KH vortices develop asymmetrically with respect to the Earth's equatorial plane where the local KH linear growth rate is maximal. Vortex driven reconnection processes take place at different latitudes, ranging from the equatorial plane to mid-latitude regions, but only in the hemisphere that turns out to be the less KH unstable. These results suggest that KH-induced reconnection is not limited to specific regions around the vortices (inside, below or above), but may be triggered over a broad and continuous range of locations in the vicinity of the vortices.

The evolution of the filamentation instability produced by two counter-streaming pair plasmas is studied with particle-in-cell (PIC) simulations in both one (1D) and two (2D) spatial dimensions. Radiation friction effects on particles are taken into account. After an exponential growth of both the magnetic field and the current density, a nonlinear quasi-stationary phase sets up characterized by filaments of opposite currents. During the nonlinear stage, a strong broadening of the particle energy spectrum occurs accompanied by the formation of a peak at twice their initial energy. A simple theory of the peak formation is presented. The presence of radiative losses does not change the dynamics of the instability but affects the structure of the particle spectra.

A covariant action principle for ideal relativistic magnetohydrodynamics (MHD) in terms of natural Eulerian field variables is given. This is done by generalizing the covariant Poisson bracket theory of Marsden et al., which uses a noncanonical bracket to effect constrained variations of an action functional. Various implications and extensions of this action principle are also discussed. Two significant by-products of this formalism are the introduction of a new divergence-free 4-vector variable for the magnetic field, and a new Lie-dragged form for the theory.

The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime.

Many astrophysical problems, ranging from structure formation in cosmology to dynamics of elliptical galaxies, refer to slow processes of evolution of essentially collisionless self-gravitating systems. In order to determine the relevant quasi-equilibrium configuration at time t from given initial conditions, it is often argued that such slow evolution may be approximated in terms of adiabatic evolution, for the calculation of which efficient semi--analytical techniques are available. Here we focus on the slow process of evolution, induced by dynamical friction of a host stellar system on a minority component of "satellites", that we have investigated in a previous paper, to determine to what extent an adiabatic description might be applied. The study is realized by comparing directly N--body simulations of the stellar system evolution (in two significantly different models) from initial to final conditions in a controlled numerical environment. We demonstrate that for the examined process the adiabatic description is going to provide incorrect answers, not only quantitatively, but also qualitatively. The two classes of models considered exhibit generally similar trends in evolution, with one exception noted in relation to the evolution of the total density profile. This simple conclusion should be taken as a warning against the indiscriminate use of adiabatic growth prescriptions in studies of structure of galaxies.

The main goal of this paper is to set up a numerical laboratory for the study of the slow evolution of the density and of the pressure tensor profiles of an otherwise collisionless stellar system, as a result of the interactions with a minority component of heavier "particles".The effects that we would like to study are those attributed to slow collisional relaxation and generically called "dynamical friction". We start by revisiting the problem of the sinking of a satellite inside an initially isotropic, non-rotating, spherical galaxy, which we follow by means of N-body simulations using about one million particles. We then consider a quasi-spherical problem, in which the satellite is fragmented into a set of many smaller masses with a spherically symmetric initial density distribution. In a wide set of experiments, designed in order to bring out effects genuinely associated with dynamical friction, we are able to demonstrate the slow evolution of the density profile and the development of a tangentially biased pressure in the underlying stellar system, while we briefly address the issue of the circularization of orbits induced by dynamical friction on the population of fragments. The results of the simulations presented here and others planned for future investigations allow us to study the basic mechanisms of slow relaxation in stellar systems and thus may be of general interest for a variety of problems, especially in the cosmological context. Here our experiments are conceived with the specific goal to clarify some mechanisms that may play a role in the evolution of an elliptical galaxy as a result of the interaction between the stars and a significant population of globular clusters or of the merging of a large number of small satellites.