Abstract.One possibility for propellantless propulsion in space is to use the momentum flux of the solar wind. A way to set up a solar wind sail is to have a set of thin long wires which are kept at high positive potential by an onboard electron gun so that the wires repel and deflect incident solar wind protons. The efficiency of this so-called electric sail depends on how large force a given solar wind exerts on a wire segment and how large electron current the wire segment draws from the solar wind plasma when kept at a given potential. We use 1-D and 2-D electrostatic plasma simulations to calculate the force and present a semitheoretical formula which captures the simulation results. We find that under average solar wind conditions at 1 AU the force per unit length is (5±1)×10 −8 N/m for 15 kV potential and that the electron current is accurately given by the well-known orbital motion limited (OML) theory cylindrical Langmuir probe formula. Although the force may appear small, an analysis shows that because of the very low weight of a thin wire per unit length, quite high final speeds (over 50 km/s) could be achieved by an electric sailing spacecraft using today's flight-proved components. It is possible that artificial electron heating of the plasma in the interaction region could increase the propulsive effect even further.Keywords. General or miscellaneous (Instruments useful in three or more fields; New fields (not classifiable under other headings); Techniques applicable in three or more fields)

This paper reviews Vlasov-based numerical methods used to model plasma in space physics and astrophysics. Plasma consists of collectively behaving charged particles that form the major part of baryonic matter in the Universe. Many concepts ranging from our own planetary environment to the Solar system and beyond can be understood in terms of kinetic plasma physics, represented by the Vlasov equation. We introduce the physical basis for the Vlasov system, and then outline the associated numerical methods that are typically used. A particular application of the Vlasov system is Vlasiator, the world’s first global hybrid-Vlasov simulation for the Earth’s magnetic domain, the magnetosphere. We introduce the design strategies for Vlasiator and outline its numerical concepts ranging from solvers to coupling schemes. We review Vlasiator’s parallelisation methods and introduce the used high-performance computing (HPC) techniques. A short review of verification, validation and physical results is included. The purpose of the paper is to present the Vlasov system and introduce an example implementation, and to illustrate that even with massive computational challenges, an accurate description of physics can be rewarding in itself and significantly advance our understanding. Upcoming supercomputing resources are making similar efforts feasible in other fields as well, making our design options relevant for others facing similar challenges.

For decades, monochromatic large‐scale ultralow frequency (ULF) waves with a period of about 30 s have been observed upstream of the quasi‐parallel bow shock. These waves typically propagate obliquely with respect to the interplanetary magnetic field (IMF), while the growth rate for the instability causing the waves is maximized parallel to the magnetic field. It has been suggested that the mechanism for the oblique propagation concerns wave refraction due to the spatial variability of the suprathermal ions, originating from the E × B drift component. We investigate the ULF foreshock under a quasi‐radial IMF with Vlasiator, which is a newly developed global hybrid‐Vlasov simulation solving the Vlasov equation for protons, while electrons are treated as a charge‐neutralizing fluid. We observe the generation of the 30 s ULF waves and compare their properties to previous literature and multipoint Time History of Events and Macroscale Interactions during Substorms (THEMIS) spacecraft observations. We find that Vlasiator reproduces the foreshock ULF waves in all reported observational aspects. We conclude that the variability of the density and velocity of the reflected back streaming ions determines the large‐scale structure of the foreshock, which affects the wave frequency, wavelength, and oblique propagation. We conclude that the wave refraction may also be at work for radial IMF conditions, which has earlier been thought of as an exception to the refraction mechanism due to the small E × B drift component. We suggest that additional refraction may be caused by the large‐scale spatial variability of the density and velocity of the back streaming ions.

No abstract

No abstract