A numerical scheme has been developed to solve the quasilinear form of the transonic stream function equation. The method is applied to compute steady two-dimensional axisymmetric solar wind-type problems. A single, perfect, non-dissipative, homentropic and polytropic gas-dynamics is assumed. The four equations governing mass and momentum conservation are reduced to a single nonlinear second order partial differential equation for the stream function. Bernoulli's equation is used to obtain a nonlinear algebraic relation for the density in terms of stream function derivatives. The vorticity includes the effects of azimuthal rotation and Bernoulli's function and is determined from quantities specified on boundaries. The approach is efficient. The number of equations and independent variables has been reduced and a rapid relaxation technique developed for the transonic full potential equation is used. Second order accurate central differences are used in elliptic regions. In hyperbolic regions a dissipation term motivated by the rotated differencing scheme of Jameson is added for stability. A successive-line-overrelaxation technique also introduced by Jameson is used to solve the equations. The nonlinear equationfor the density is a double valued function of the stream function derivatives. The velocities are extrapolated from upwind points to determine the proper branch and Newton's method is used to iteratively compute the density. This allows accurate solutions with few grid points. The applications first illustrate solutins to solar wind models. The equations predict that the effects of vorticity must be confined near the surface and far away the streamlines must resemble the spherically symmetric solution. Irrotational and rotational flows show this behavior. The streamlines bend toward the rotation axis for rapidly rotating models because the coriolis force is much larger than the centrifugal force. Models of galactic winds are computed by considering the flow exterior to a surface which surrounds a uniform density oblate spheroid. Irrotational results with uniform outward mass flux show streamlines bent toward the equator and nearly spherical sonic surfaces. Rotating models for which Bernoulli's function is not constant show the sonic surface is deformed consistent with the one-dimensional theory.

The solar wind interaction with an icy comet is studied through a model problem. A hybrid simulation is done of a box with evenly distributed water ions and protons, where initially the water ions are stationary, and protons move with the speed of the solar wind. The purpose of the thesis is to investigate the interaction between the two species through the convective electric field, and focus is on early acceleration of pick-up ions, and deflection of the solar wind. It is relevant to the cometary case, because it enables study of the physics of this interaction, without involving other mechanisms, such as bow shock, magnetic field pile-up and draping. The species are found to exchange kinetic energy similar to a damped oscillator, where the dampening is caused by kinetic energy being transferred to the magnetic field. At early times, i.e. times smaller than the gyration time for the water ions, the solar wind does not lose much speed when it is deflected. For comparable number densities, the solar wind can be deflected more than 90° at early times, and loses more speed, and water ions are picked up faster. The total kinetic energy of the system decreases when energy builds up in the magnetic field. The nature of the energy exchange is strongly dependent on the number density ratio between water ions and protons. A density instability with behaviour similar to a plasma beam instability forms as energy in the magnetic field increases, and limits the amount of time the simulation preserves total energy, for the particular hybrid solver used. There is a discussion on the structure of the density instability, and it is compared to cometary simulations.

La thèse porte sur le rôle de la rotation planétaire dans la structure globale de l'interaction vent solaire/magnétosphère à partir de simulations magnétohydrodynamiques (MHD). Les magnétosphères planétaires du système solaire présentent une incroyable diversité, et notamment dans leurs configurations respectives de l'inclinaison de leur axe magnétique par rapport à leur axe de rotation. La durée des périodes de rotation par rapport au temps de relaxation de chaque magnétosphère diffère aussi d'une planète à l'autre. On distingue ainsi les rotateurs lents (Mercure et la Terre), pour lesquels le temps de relaxation est plus court que la période de rotation, des rotateurs rapides (Jupiter, Saturne, Uranus et Neptune). Dans le cas du rotateur lent Mercure, on s'intéresse à l'influence des paramètres du vent solaire sur la structure globale du champ magnétique et de l'écoulement. En appui à la mission spatiale BepiColombo, nous présentons des simulations effectuées pour deux modèles différents de champ magnétique herméen. Nous détaillons le rôle des fronts d'onde MHD stationnaires, en particulier les fronts stationnaires de mode lent dans la magnétogaine. Saturne présente la particularité d'avoir un axe magnétique parfaitement aligné avec son axe de rotation. C'est donc un cas de rotateur rapide stationnaire, qui nous permet d'étudier la structure globale du champ magnétique et de l'écoulement pour différentes orientations de l'IMF, mais aussi pour différentes vitesses de rotation de la planète. Enfin, le cas d'une configuration quelconque, avec un grand angle entre l'axe magnétique et l'axe de rotation planétaire, est étudié en présence d'un vent solaire magnétisé en s'inspirant de la configuration d'Uranus au solstice et à l'équinoxe. Dans la configuration « solstice », c'est à dire lorsque l'axe de rotation pointe vers le Soleil, on montre qu'une structure de nature alfvénique se développe en hélice dans la queue de la magnétosphère, et que les zones de reconnexion entre le champ magnétique planétaire et l'IMF, qui forment aussi une double hélice, ralentissent la progression de la structure alfvénique. A l'équinoxe, lorsque l'axe de rotation est toujours dans le plan de l’écliptique mais perpendiculaire à la direction Soleil-Uranus, la structure en hélice disparaît. / The topic of the thesis is the part of planetary rotation in the global structure of the solar wind interaction with planetary magnetospheres using MHD simulations. We discuss the distinction between slow and fast rotators from a MHD point of view. In the case of a non-rotating magnetosphere (as is the one of Mercury), the part of standing MHD modes is studied, along with a method to identify them in simulations. A fast-rotating but stationary magnetosphere (inspired by the case of Saturn) is presented in details and provides a good test to validate the new version of the AMRVAC code allowing for any configuration regarding the respective directions of the planetary spin axis, planetary magnetic axis, solar wind inflow direction, and IMF orientation. Finally, a random configuration, with a large angle between the planetary spin and magnetic axis, is analyzed for the first time in presence of a magnetized solar wind, using configurations inspired from the planet Uranus at solstice and equinox.

In space and astrophysical plasmas, turbulence is responsible for transferring energy from large scales driven by violent events or instabilities, to smaller scales where turbulent energy is ultimately converted into plasma heat by dissipative mechanisms. In the inertial range, the self-similar turbulent energy cascade to smaller spatial scales is driven by the nonlinear interaction between counterpropagating Alfvén waves, denoted Alfvén wave collisions. For the more realistic case of the collision between two initially separated Alfvén wavepackets (rather than previous idealized, periodic cases), we use a nonlinear gyrokinetic simulation code, AstroGK, to demonstrate three key properties of strong Alfvén wave collisions: they (i) facilitate the perpendicular cascade of energy and (ii) generate current sheets self-consistently, and (iii) the modes mediating the nonlinear interaction are simply Alfvén waves. Once the turbulent cascade reaches the ion gyroradius scale, the Alfvén waves become dispersive and the turbulent energy starts to dissipate, energizing the particles via wave-particle interactions with eventual dissipation into plasma heat. The novel Field-Particle Correlation technique determines how turbulent energy dissipates into plasma heat by identifying which particles in velocity-space experience a net gain of energy. By utilizing knowledge of discrete particle arrival times, we devise a new algorithm called PATCH (Particle Arrival Time Correlation for Heliophysics) for implementing a field-particle correlator onboard spacecraft. Using AstroGK, we create synthetic spacecraft data mapped to realistic phase-space resolutions of modern spacecraft instruments. We then utilize Poisson statistics to determine the threshold number of particle counts needed to resolve the velocity-space signature of ion Landau damping using the PATCH algorithm.

<p>The Neutral Particle Detector (NPD) is a new type of instrumentation for energetic neutral atom (ENA) diagnostics. This thesis deals with development of the NPD sensor designed as a part of the plasma and neutral particle packages ASPERA-3 and ASPERA-4 on board Mars Express and Venus Express, the European Space Agency (ESA) satellites to Mars and Venus, respectively. It describes how the NPD sensors were designed, developed, tested and calibrated. </p><p>It also presents the first scientific results obtained with NPD during its operation at Mars. </p><p>The NPD package consists of two identical detectors, NPD1 and NPD2. Each detector has a 9^o x 90^o intrinsic field-of-view divided into three sectors. The ENA detection principle is based on the surface interaction technique. NPD detects ENA differential fluxes within the energy range of 100 eV to 10 keV and is capable of resolving hydrogen and oxygen atoms by time-of-flight (TOF) measurements or pulse height analysis.</p><p>During the calibration process the detailed response of the sensor was defined, including properties such as an angular response function and energy dependent efficiency of each of the sensor sectors for different ENA species. </p><p>Based on the NPD measurements at Mars the main scientific results reported so far are:</p><p>- observation of the Martian H-ENA jet / cone and its dynamics, </p><p>- observations of ENA emissions from the Martian upper atmosphere, </p><p>- measurements of the hydrogen exosphere density profile at Mars, </p><p>- observations of the response of the Martian plasma environment to an interplanetary shock, </p><p>- observations of the H-ENA fluxes in the interplanetary medium.</p>

Thesis (M.Sc. (Physics))--North-West University, Potchefstroom Campus, 2008.

No description available.

The Neutral Particle Detector (NPD) is a new type of instrumentation for energetic neutral atom (ENA) diagnostics. This thesis deals with development of the NPD sensor designed as a part of the plasma and neutral particle packages ASPERA-3 and ASPERA-4 on board Mars Express and Venus Express, the European Space Agency (ESA) satellites to Mars and Venus, respectively. It describes how the NPD sensors were designed, developed, tested and calibrated. It also presents the first scientific results obtained with NPD during its operation at Mars. The NPD package consists of two identical detectors, NPD1 and NPD2. Each detector has a 9o x 90o intrinsic field-of-view divided into three sectors. The ENA detection principle is based on the surface interaction technique. NPD detects ENA differential fluxes within the energy range of 100 eV to 10 keV and is capable of resolving hydrogen and oxygen atoms by time-of-flight (TOF) measurements or pulse height analysis. During the calibration process the detailed response of the sensor was defined, including properties such as an angular response function and energy dependent efficiency of each of the sensor sectors for different ENA species. Based on the NPD measurements at Mars the main scientific results reported so far are: - observation of the Martian H-ENA jet / cone and its dynamics, - observations of ENA emissions from the Martian upper atmosphere, - measurements of the hydrogen exosphere density profile at Mars, - observations of the response of the Martian plasma environment to an interplanetary shock, - observations of the H-ENA fluxes in the interplanetary medium.

The Sun’s corona expands outward, populating the solar system with plasma. This plasma is known as the solar wind. The solar wind carries with it the Sun’s magnetic field, which is also known as the interplanetary magnetic field (IMF). The resulting configuration of the IMF creates a current sheet at solar equatorial latitudes, which the Earth crosses as it orbits the Sun. When the Earth is on one side of the current sheet it is in a sector where the IMF is directed largely away from or toward the Sun. On the other side of the current sheet the IMF is in opposite direction. The crossing of the current sheet is known as a sector boundary crossing (SBC). The solar wind and IMF properties change significantly near the current sheet, and this affects the Earth’s ionosphere. The Super Dual Auroral Radar Network (SuperDARN) high frequency (HF) radar data rates from 2001-2011 were examined using several techniques: a superposed epoch analysis, a fast fourier transform (FFT) analysis, and a cross–correlation analysis. Data from multiple instruments were analyzed in this study. These include the solar wind and IMF data from spacecraft, observations of charged particles precipitating into the Earth’s ionosphere, echoes from ground–based SuperDARN radars, and data from gound–based neutron monitors that detect galactic cosmic rays. Solar wind and IMF properties change significantly across a sector boundary. An increase in the IMF magnitude of about 30% occurs on the day of the SBC, and the IMF returns to pre–crossing values over the next two days. There is a decrease in the solar wind speed of about 15% the day before and the day of the SBC, and the solar wind density doubles at the time of the SBC. The polarity of the SBC does not appear to affect the solar wind and IMF. A peak in the data rate of SuperDARN echoes from both the ionosphere and ground occurs within one day of the SBC, though the variability of these data is quite large. The hemispherical power, which is an estimation of the electron energy flux precipitating into the ionosphere derived from satellite observations, increases following a SBC. Satellite particle data also revealed that the equatorward auroral oval boundary moves equatorward following a SBC. The cosmic ray counts at the Earth’s surface appear to be unaffected by the SBC. The solar wind and ionosphere data sets exhibited strong periodicities, and these were harmonics of the synodic rotational period of the Sun (approximately 27 days). Common periodicities observed were 27 days, 13.5 days, 9 days, 6.75 days and 5.4 days. There was a dominant 9–day periodicity observed in the solar wind and ionospheric data from 2005–2008, but was not observed in the solar 10.7 cm wavelength electromagnetic flux. The 9-day periodicity in the solar wind during this period has been linked to three persistent features on the Sun that produced corotating high–speed streams, or areas of fast solar wind. The parameters whose change did not depend on the polarity of the SBC had periodicities that were half that of the SBCs. From the cross–correlation analysis some relationships between the data sets became evident. For periods of high solar wind speed there were low SuperDARN data rates, and vice versa. The solar wind speed and hemispherical power were found to be well correlated, while the hemispherical power and the SuperDARN scatter occurrence were found to be anticorrelated. The solar wind changes appear to be affecting the state of the ionosphere, likely through particle precipitation. The SuperDARN scatter occurrence has been shown in past studies to be most greatly affected by changes in the electron density profile of the ionosphere, which can be influenced by changes in particle precipitation. These results demonstrate a link between the solar wind and the state of the ionosphere.

No description available.