Search

The test particle approach is a widely used method for studying the dynamics of charged particles in complex electromagnetic fields and has been successful in explaining particle energization in turbulent plasmas. However, this approach is fundamentally not self-consistent, as test particles do not generate their own electromagnetic fields and therefore do not interact with their surroundings realistically. In this work, we compare the energization of a population of test protons in a magnetofluid to that of a plasma composed of self-consistent particles. We use a compressible Hall magnetohydrodynamic (CHMHD) model for the test particle case and a hybrid particle-in-cell (HPIC) approach for the self-consistent case, conducting both 2D and 3D simulations. We calculate the rate of energization and conversion to thermal energy in both models, finding a higher temperature for the test particle case. Additionally, we examine the distribution of suprathermal particles and find that, in the test particle scenario, these particles eventually occupy the entire domain, while in the self-consistent case, suprathermal particles are confined to specific regions. We conclude that while test particles capture some qualitative features of their self-consistent counterparts, they miss finer phenomena and tend to overestimate energization.

Direct numerical simulations of 3D compressible MHD turbulence were performed in order to study the relation between waves modes and coherent structures and the consequent energization of test particles. Moreover, the question of which is the main mechanism of this particle energization is rigorously discussed. In particular, using the same initial conditions, we analyzed the non-linear and linear evolution of a turbulent state along with the case of randomized phases. Then, the behavior of the linear and non-linear simulations were compared through the study of time evolution of particle kinetic energy and preferential concentration. Also, spatio temporal spectra were used to identify the presence of wave modes and quantify the fraction of energy around the MHD modes in linear and non-linear simulations. Finally, the variation of the correlation time of the external forcing is studied in detail along with the effect on the particle energization (and clustering) and the presence of wave modes. More specifically, particle energization tends to decrease when the fraction of linear energy increase, supporting the idea that energization by structures is the dominant mechanism for particle energization instead of resonating with wave modes as suggested by Fermi energization theory.

Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer heliosphere. This shortcoming leaves numerous central questions incompletely answered. Deficiencies remain in at least two important subjects, Space Weather and fundamental plasma turbulence theory, due to a lack of a more complete understanding of the space-time structure of dynamic plasmas. Only with multispacecraft measurements over suitable spans of spatial separation and temporal duration can these ambiguities be resolved. We note that these characterizations apply to turbulence across a wide range of scales, and also equally well to shocks, flux ropes, magnetic clouds, current sheets, stream interactions, etc. In the following, we will describe the basic requirements for resolving space-time structure in general, using turbulence' as both an example and a principal target or study. Several types of missions are suggested to resolve space-time structure throughout the Heliosphere., Comment: White Paper submitted to: Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033. arXiv admin note: substantial text overlap with arXiv:1903.06890

In this work we perform direct numerical simulations of three-dimensional magnetohydrodynamics with a background magnetic field, representing solar wind plasma, and introduce test particles to explore how a turbulent electromagnetic environment affects them. Our focus is on the terms of the electric field present in the generalized Omh's Law that are usually dismissed as unimportant. These are the Hall and the electronic pressure (EP) terms, but we concentrate primarily on the latter. We discover that the EP term generates an acceleration of the particles, which represent protons, in the direction parallel to the background magnetic field, in contrast to the known preferential perpendicular energization. By studying the electric field itself, we are able to detect the type of structures of the EP field that produce such parallel acceleration. These are thin and elongated structures placed on top of a monotonic and near-zero background. A statistical study to understand the real significance of the electronic pressure term is also performed.

We investigated incompressible and compressible magnetohydrodynamic (MHD) energy cascade rates in the solar wind at different heliocentric distances. We used in situ magnetic field and plasma observations provided by the Parker Solar Probe (PSP) mission and exact relations in fully developed turbulence. To estimate the compressible cascade rate, we applied two recent exact relations for compressible isothermal and polytropic MHD turbulence, respectively. Our observational results show a clear increase of the absolute value of the compressible and incompressible cascade rates as we get closer to the Sun. Moreover, we obtained an increase in both isothermal and polytropic cascade rates with respect to the incompressible case as compressibility increases in the plasma. Further discussion about the relation between the compressibility and the heliocentric distance is carried out. Furthermore, we compared both exact relations as compressibility increases in the solar wind and although we note a slightly trend to observe larger cascades using a polytropic closure, we obtained essentially the same cascade rate in the range of compressibility observed. Finally, we investigated the signed incompressible and compressible energy cascade rates and its connection with the real cascade rate., Comment: 23 pages. arXiv admin note: text overlap with arXiv:2102.11781

We present a high-order numerical study of the dependence of the dynamo onset on different electromagnetic boundary conditions, in convecting Boussinesq flows forced by a temperature gradient. Perfectly conducting boundaries, vacuum, and mixed electromagnetic boundary conditions are considered, using a method that treats fields and boundary conditions with close to spectral accuracy. Having one or two conducting boundaries greatly facilitates dynamo action. For the mixed case it is shown that the critical magnetic Reynolds number becomes independent of the Rayleigh number, Ra, for sufficiently large Ra.

The effect of precession in a rotating sphere filled with fluid was studied with direct numerical simulations, both in the incompressible hydrodynamics (HD) and magnetohydrodynamics (MHD) scenarios. In both cases the asymptotic state and its dependence with both rotating and precession frequency was analyzed. For the MHD case no self-sustaining dynamos were found for the prograde precession case, whereas on the other hand a critical retrograde precession frequency was found above which dynamo action is self-sustained. It was also found that these correspond to small-scale dynamos with a developed turbulent regime. Furthermore, it is observed the presence of reversals of the magnetic dipole moment with greater waiting times between reversals for smaller precession frequencies.

We consider the simplest optimal control problem with one nonregular mixed inequality constraint, i.e. when its gradient in the control can vanish on the zero surface. Using the Dubovitskii--Milyutin theorem on the approximate separation of convex cones, we prove a first or der necessary condition for a weak minimum in the form of the so-called local minimum principle, which is formulated in terms of functions of bounded variation, integrable functions, and Lebesgue--Stieltjes measures, and does not use functionals on the space of measurable bounded functions. Two illustrative examples are given. The work is based on results by Milyutin., Comment: Text in 26 pages, no figures

Various exact laws governing compressible magnetohydrodynamic (MHD) and Hall-MHD (CHMHD) turbulence have been derived in recent years. Other than their fundamental theoretical interest, these laws are generally used to estimate the energy dissipation rate from spacecraft observations in order to address diverse problems related, e.g., to heating of the solar wind (SW) and magnetospheric plasmas. Here we use various $1024^3$ direct numerical simulation (DNS) data of free-decay isothermal CHMHD turbulence obtained with the GHOST code (Geophysical High-Order Suite for Turbulence) to analyze two of the recently derived exact laws. The simulations reflect different intensities of the initial Mach number and the background magnetic field. The analysis demonstrates the equivalence of the two laws in the inertial range and relates the strength of the Hall effect to the amplitude of the cascade rate at sub-ion scales. When taken in their general form (i.e., not limited to the inertial range) some subtleties regarding the validity of the stationarity assumption or the absence of the forcing in the simulations are discussed. We show that the free-decay nature of the turbulence induces a shift from a large scale forcing towards the presence of a scale-dependent reservoir of energy fueling the cascade or dissipation. The reduced form of the exact laws (valid in the inertial range) ultimately holds even if the stationarity assumption is not fully verified., Comment: 14 pages, 10 figures

In the present work, we study the energization and displacement of heavy ions through the use of test particles interacting with the electromagnetic fields of magnetohydrodynamic (MHD) turbulence. These fields are obtained from pseudospectral direct numerical solutions (DNSs) of the compressible three-dimensional MHD equations with a strong background magnetic field. We find particle energization to be predominantly perpendicular as the ions become heavier (lower charge-to-mass ratio) and that high displacement is detrimental for perpendicular energization. On the other hand, perpendicular displacement is unaffected by the charge-to-mass ratio, which we explain with a simple guide center model. Using Voronoi tessellation along with this model, we analyze preferential concentration and find that particles behave as tracers in the perpendicular plane, clustering in regions with $\nabla_\perp\cdot\mathbf{u}_\perp < 0$. These regions also have $(\nabla\times\mathbf{E})_z < 0$, which is optimal for perpendicular energization, thus providing a mechanism to understand precedent results.

Using direct numerical simulations (DNSs), the interaction between linear waves and turbulence under the compressible magnetohydrodynamic (CMHD) approach was studied. A set of DNSs in three dimensions for a spatial resolution of $128^3$ and $256^3$ were performed. A parametric study was carried out varying the sonic Mach number, the mean magnetic field and the compressibility amplitude of the forcing. Spatio-temporal spectra of the magnetic energy were built and analyzed, allowing for direct identification of all wave modes in a CMHD turbulent system and quantification of the amount of energy in each mode as a function of the wave number. Thus, linear waves were detected, that is Alfv\'en waves and fast and slow magnetosonic waves. Furthermore, different responses of the plasma were found according to whether the Mach number or the mean magnetic field was varied. On the other hand, making use of spatio-temporal spectra and two different integration methods, we accurately quantified the amount of energy present in each of the normal modes. Finally, although the presence of linear waves was observed, in all the cases studied the system was mainly dominated by the non-linear dynamics of the plasma., Comment: arXiv admin note: substantial text overlap with arXiv:1707.09190

A high-order method to evolve in time electromagnetic and velocity fields in conducting fluids with non-periodic boundaries is presented. The method has a small overhead compared with fast FFT-based pseudospectral methods in periodic domains. It uses the magnetic vector potential formulation for accurately enforcing the null divergence of the magnetic field, and allowing for different boundary conditions including perfectly conducting walls or vacuum surroundings, two cases relevant for many astrophysical, geophysical, and industrial flows. A spectral Fourier continuation method is used to accurately represent all fields and their spatial derivatives, allowing also for efficient solution of Poisson equations with different boundaries. A study of conducting flows at different Reynolds and Hartmann numbers, and with different boundary conditions, is presented to study convergence of the method and the accuracy of the solenoidal and boundary conditions.

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence., Comment: White paper submitted to the PLASMA 2020 Decadal Survey Committee

We consider the time-optimal problem for a classical system of "double integrator" under the presence of a linear state constraint. By using the Maximum Principle of Dubovitskii and Milyutin, we determine a complete synthesis of optimal trajectories and qualitatively analyze the corresponding Lagrange multipliers., Comment: 27 pages, 9 figures

We present a Fourier Continuation-based parallel pseudospectral method for incompressible fluids in cuboid non-periodic domains. The method produces dispersionless and dissipationless derivatives with fast spectral convergence inside the domain, and with very high order convergence at the boundaries. Incompressibility is imposed by solving a Poisson equation for the pressure. Being Fourier-based, the method allows for fast computation of spectral transforms. It is compatible with uniform grids (although refined or nested meshes can also be implemented), which in turn allows for explicit time integration at sufficiently high Reynolds numbers. Using a new parallel code named SPECTER we illustrate the method with two problems: channel flow, and plane Rayleigh-B\'enard convection under the Boussinesq approximation. In both cases the method yields results compatible with previous studies using other high-order numerical methods, with mild requirements on the time step for stability.

We study the spatio-temporal behavior of the Els\"asser variables describing magnetic and velocity field fluctuations, using direct numerical simulations of three-dimensional magnetohydrodynamic turbulence. We consider cases with relatively small, intermediate, and large values of a mean background magnetic field, and with null, small, and high cross-helicity (correlations between the velocity and the magnetic field). Wavenumber-dependent time correlation functions are computed for the different simulations. From these correlation functions, the decorrelation time is computed and compared with different theoretical characteristic times: the local non-linear time, the random-sweeping time, and the Alfv\'enic time. It is found that decorrelation times are dominated by sweeping effects for low values of the mean magnetic field and for low values of the cross-helicity, while for large values of the background field or of the cross-helicity and for wave vectors sufficiently aligned with the guide field, decorrelation times are controlled by Alfv\'enic effects. Finally, we observe counter-propagation of Alfv\'enic fluctuations due to reflections produced by inhomogeneities in the total magnetic field. This effect becomes more prominent in flows with large cross-helicity, strongly modifying the propagation of waves in turbulent magnetohydrodynamic flows., Comment: 16 pages, 40 figures

We consider some critical claims concerning our above paper, and reply to these claims.

In this work, direct numerical simulations of the compressible fluid equations in turbulent regimes are performed. The behavior of the flow is either dominated by purely turbulent phenomena or by the generation of sound waves in it. Previous studies suggest that three different types of turbulence may happen at the low Mach number limit in polytropic flows: Nearly incompressible (NI), modally equipartitioned compressible (MEC) and compressible wave (CW). The distinction between these types of turbulence is investigated here applying different kinds of forcing. Scaling of density fluctuations with Mach number, comparison of the ratio of transverse velocity fluctuations to longitudinal fluctuations and spectral decomposition of fluctuations are used to distinguish the nature of these solutions. From the study of the spatio-temporal spectra and correlation times we quantify the contribution of the waves to the total energy of the system. Also, in the dynamics of a compressible flow, three associated correlation times are considered: the non-linear time of local interaction between scales, the sweeping time or non-local time of large scales on small scales, and the time associated with acoustic waves (sound). We observed that different correlation times dominate depending on the wave number (k), the Mach number and the type of forcing.

Plasma dynamics is a multi-scale problem that involves many spatial and temporal scales. Turbulence connects the disparate scales in this system through a cascade that is established by nonlinear interactions. Most astrophysical plasma systems are weakly collisional, making a fully kinetic Vlasov description of the system essential. Use of reduced models to study such systems is computationally desirable but careful benchmarking of physics in different models is needed. We perform one such comparison here between fully kinetic Particle-In-Cell (PIC) model and a two-fluid model that includes Hall physics and electron inertia, with a particular focus on the sub-proton scale electric field. We show that in general the two fluid model captures large scale dynamics reasonably well. At smaller scales the Hall physics is also captured reasonably well by the fluid code but electron features show departures from the fully kinetic model. Implications for use of such fluid models are discussed.

A parametric study of the magnetic dipole behavior in resistive incompressible MHD inside a rotating sphere is performed, using direct numerical simulations and considering Reynolds and Ekman numbers as controlling parameters. The tendency is to obtain geodynamo-like magnetic dipole reversal regimes for sufficiently small Ekman and large Reynolds numbers. The typical dipole latitude obtained in the reversal regime is around 40 degrees (with respect to the rotation axis of the sphere). A statistical analysis of waiting times between dipole reversals is also performed, obtaining a non-Poissonian distribution of waiting times, indicating long-term memory effects. We also report the presence of a $1/f$ frequency power spectrum in the magnetic dipole time-series, which also shows a tendency to grow toward lower frequencies as the Ekman number is decreased., Comment: 14 pages, 9 figures

Three-dimensional direct numerical simulations are used to study the energy cascade rate in isothermal compressible magnetohydrodynamic turbulence. Our analysis is guided by a two-point exact law derived recently for this problem in which flux, source, hybrid, and mixed terms are present. The relative importance of each term is studied for different initial subsonic Mach numbers $M_S$ and different magnetic guide fields ${\bf B}_0$. The dominant contribution to the energy cascade rate comes from the compressible flux, which depends weakly on the magnetic guide field ${\bf B}_0$, unlike the other terms whose modulus increase significantly with $M_S$ and ${\bf B}_0$. In particular, for strong ${\bf B}_0$ the source and hybrid terms are dominant at small scales with almost the same amplitude but with a different sign. A statistical analysis made with an isotropic decomposition based on the SO(3) rotation group is shown to generate spurious results in presence of ${\bf B}_0$, when compared with an axisymmetric decomposition better suited to the geometry of the problem. Our numerical results are compared with previous analyses made with in-situ measurements in the solar wind and the terrestrial magnetosheath.

Using spatio-temporal spectra we show direct evidence of excitation of magnetosonic and Alfv\'en waves in three-dimensional compressible magnetohydrodynamic turbulence at small Mach numbers. For the plasma pressure dominated regime, or high $\beta$ regime (with $\beta$ the ratio between fluid and magnetic pressure), and for the magnetic pressure dominated regime, or low $\beta$ regime, we study magnetic field fluctuations parallel and perpendicular to a guide magnetic field $\textbf{B}_0$. In the low $\beta$ case we find excitation of compressible and incompressible fluctuations, with a transfer of energy towards Alfv\'enic modes and to a lesser extent towards magnetosonic modes. In particular, we find signatures of the presence of fast magnetosonic waves in a scenario compatible with that of weak turbulence. In the high $\beta$ case, fast and slow magnetosonic waves are present, with no clear trace of Alfv\'en waves, and a significant part of the energy is carried by two-dimensional turbulent eddies.

We consider a class of optimal control problems with a state constraint and investigate a trajectory with a single boundary interval (subarc). Following R.V. Gamkrelidze, we differentiate the state constraint along the boundary subarc, thus reducing the original problem to a problem with mixed control-state constraints, and show that this way allows one to obtain the full system of stationarity conditions in the form of A.Ya. Dubovitskii and A.A. Milyutin, including the sign definiteness of the measure (state constraint multiplier), i.e., the nonnegativity of its density and atoms at junction points. The stationarity conditions are obtained by a two-stage variation approach, proposed in this paper. At the first stage, we consider only those variations, which do not affect the boundary interval, and obtain optimality conditions in the form of Gamkrelidze. At the second stage, the variations are concentrated on the boundary interval, thus making possible to specify the stationarity conditions and obtain the sign of density and atoms of the measure., Comment: 30 pages

In this paper we analyze the effect of dynamical three-dimensional MHD turbulence on test particle acceleration, and compare how this evolving system affects particle energization by current sheets interaction, against frozen-in-time fields. To do this we analize the ensamble particle acceleration for static electromagnetic fields extracted from direct numerical simulations of the MHD equations, and compare with the dynamical fields. We show that a reduction in particle acceleration in the dynamical model results from the particle trapping in the field lines, which forces the particles to remain in a moving current sheet that suppress the longer exposure at the strong electric field gradients located between structures, which is an efficient particle acceleration mechanism. In addition, we analize the effect of anisotropy caused by the mean magnetic field. It is well known that for sufficiently strong external fields, the system suffers a transition towards a two-dimensional flow. This causes an increment in the size of the coherent structures, resulting in a magnetized state of the particles and the reduction of the particle energization., Comment: 13 pages,7 figures

The investigation of the diffusive transport of charged particles in a turbulent magnetic field remains a subject of considerable interest. Research has most frequently concentrated on determining the diffusion coefficient in the presence of a mean magnetic field. Here we consider diffusion of charged particles in fully three-dimensional isotropic turbulent magnetic fields with no mean field, which may be pertinent to many astrophysical situations. We identify different ranges of particle energy depending upon the ratio of the Larmor radius of the charged particle to the characteristic outer length scale of the turbulence. Two different theoretical models are proposed to calculate the diffusion coefficient, each applicable to a distinct range of particle energies. The theoretical results are compared with those from computer simulations, showing good agreement.