The turbulen flow in an asymmetric diffuser has been en studied by the means of Reynold average Navier-Stokes equations with both differential and explict algebraic expressions to model the Reynolds stress tensor. Modifications to the differential stress model have been derived, using the inverse turbulence timescale to obtain the dissipation of turbuence kinetic energy. The explicit algebraic Reynolds stress model has been used in combination with a two-equation platform to close the system of equations. Modifications made to the transport equation for the inverse turbulence timescale has made it possible to substantially relax the deman on near-wall resolution of this quantity. The rapid growth wth present in the original formulation can be treated as an explicit function of the wall-normal distance. In order to use the new formulation for the transport equation, an equation has as been derived to obtain the shortest distance bettwee a point and the closest wall, regardles of the geometric complexity of the domain. An explicit algebraic expression to model the passive scalar flux vector has been investigated using a comparison with a standard eddy-diffusivity model in the asymmetric diffuser. Results show a substantial improvement of the complexity of the scalar field and scalar flux vector in sepaarated flows. Automated code generation has been used in all the above studies to generate versatile model testing tools for general two-dimensional geometries. Finite element formulations are used for these tools.

Broadband extremely low-frequency wave emissions below the ion plasma frequency have been observed by a number of spacecraft and rockets on auroral field lines. The importance of these broadband emissions for transverse ion heating and electron acceleration in the auroral regions is now reasonably well established. However, the exact mechanism(s) for mediating this energy transfer and the wave mode(s) involved are not well known. In this thesis we focus on the identification of broadband waves by different methods. Two wave analysis methods, involving different approximations and assumptions, give consistent results concerning the wave mode identification. We find that much of the broadband emissions can be identified as a mixture of ion acoustic, electrostatic ion cyclotron and, ion Bernstein waves, which all can be described as different parts of the same dispersion surface in the linear theory of waves in homogeneous plasma. A new result is that ion acoustic waves occur on auroral magnetic field lines. These are found in relatively small regions interpreted as acceleration regions without cold (tens of eV) electrons. From interferometry we also determine the phase velocity and k vector for parallel and oblique ion acoustic waves. The retrieved characteristic phase velocity is of the order of the ion acoustic speed and larger than the thermal velocity of the protons. The typical wavelength is around the proton gyro radius and always larger than the Debye length which is consistent with ion acoustic waves. We have observed quasi-static parallel electric fields associated with the ion acoustic waves in regions with large-scale currents. Waves, in particular ion acoustic waves, can create an anomalous resistivity due to wave-particle interaction when electrons are retarded or trapped by the electric wave-field. To maintain the large-scale current, a parallel electric field is set up, which then can accelerate a second electron population to high velocities.

This thesis deals with energetic oxygen ions (i.e. single-charged atomic oxygen ions, O+) at altitudes higher than 5 Earth radii (RE) and at latitudes above 75 (toward 90) degrees invariant latitude (deg ILAT) in the dayside polar magnetosphere observed by Cluster. The instrument used in this study is CIS (Cluster Ion Spectrometry experiment) / CODIF (a time-of-flight ion COmposition and DIstribution Function analyser), which covers an energy range from »10 eV up to 38 keV. Cluster detected O+ with energies more than 1 keV (hereafter termed “keV O+”), indicating that energization and/or acceleration process(es) take place in the dayside high-altitude (inside magnetopause) and high-latitude region. These O+ are outflowing (precisely, upward-going along the geomagnetic field lines), and these outflowing keV O+ show a heated (or energized) signature in the velocity distribution as well. First, outflowing O+ are observed at the poleward cusp and/or the mantle formed a partial shell-like configuration seen in the velocity distribution. Second, the latitudinal distribution of outflowing O+ (most of them have energies less than 1 keV statistically) observed below 7 RE is consistent with velocity filter effect by the polar convection, while the latitudinal distribution of outflowing keV O+ observed above 7 RE cannot be explained by velocity filter effect only, i.e. this indicates that additional energization and/or acceleration takes place at higher altitudes in the dayside polar region. Thirdly, a tendency to observe outflowing keV O+ for during different geomagnetic conditions is studied. The keV O+ above 9 RE is more often for K p¸5 rather than for K p•3. However the energy of O+ is not dependent on ASY /SYM indices. Finally, the dependence on the solar wind conditions is also studied. The energization and/or acceleration of outflowing O+ is controlled by both solar wind moments (except solar wind electric field) and strong southward interplanetary magnetic field (IMF) at the time scale of tens of minutes at only higher altitudes. Further examination shows that solar wind dependence is different at three regions: one is the poleward cusp, another is the low-altitude polar cap, and finally the high-altitude polar cap, combining all the results. There is (a) new energization and/or acceleration process(es) at the high-altitude polar cap. On the other hand, flux enhancement of O+ observed above 5 RE is also controlled by solar wind moments (e.g. solar wind electric field) and strong southward IMF, however the ionospheric changes play a more important role on the flux enhancement of O+.

Observations of modifications of the electron temperature in the F-region produced by powerful high-frequency waves transmitted in X-mode are presented. The experiments were performed during quiet nighttime conditions with low ionospheric densities so no reflections occurred. Nevertheless temperature enhancements of the order of 300-400K were obtained. The modifications found can be well described by the theory of Ohmic heating by the pump wave and both temporal and spatial changes are reproduced.  A brief overview of several different experimental campaigns at EISCAT facilities in the period from October 2006 to February 2008 are also given pointing out some interesting features from the different experiments. The main focus is then on the campaign during October 2006 and modifications of the electron temperature in the F-region.

This thesis is focused on the role of Alfvén waves in the energy transformation and transport in the magnetosphere. Different aspects of Alfvén wave generation, propagation and dissipation are considered. The study involves analysis of experimental data from the Freja, Polar and Cluster spacecraft, as well as theoretical development. An overview of the linear theory of Alfvén waves is presented, including the effects of fnite parallel electron inertia and fnite ion gyroradius, and nonlinear theory is developed for large amplitude Alfvén solitons and structures. The methodology is presented for experimental identification of dispersive Alfvén waves in a frame moving with respect to the plasma, which facilitates the resolution of the space-time ambiguity in such measurements. Dispersive Alfvén waves are identified on field lines from the topside ionosphere up to the magnetopause and it is suggested they play an important role in magnetospheric physics. One of the processes where Alfvén waves are important is the establishment of the field aligned current system, which transports the energy from the reconnection regions at the magnetopause to the ionosphere, where a part of the energy is dissipated. The main mechanism for the dissipation in the top-side ionosphere is related to wave-particle interactions leading to particle energization/heating. An observed signature of such a process is the presence of parallel energetic electron bursts associated with dispersive Alfvén waves. The accelerated electrons (electron beams) are unstable with respect to the generation of high frequency plasma wave modes. Therefore this thesis also demonstrates an indirect coupling between low frequency Alfvén wave and high frequency oscillations.

The space around Earth, confined in the terrestrial magnetosphere, is to some extent shielded from the Sun's solar wind plasma and magnetic field. During certain conditions, however, strong interaction can occur between the solar wind and the magnetosphere, resulting in magnetospheric activity of several forms, among which substorms and storms are the most prominent. A general framework for how these processes work have been outlayed through the history of research, however, there still remain questions to be answered. The most striking example regards the onset of substorms, where both the onset cause and location in the magnetosphere/ionosphere are still debated. These are clearly not easily solved problems, since a substorm is a global process, ideally requiring simultaneous measurements in the magnetotail and ionosphere. Investigated in this work are temporal and spatial scales for substorm and convection processes in the Earth's magnetosphere and ionosphere. This is performed by combining observations from a number of both ground-based and spacecraft-borne instruments. The observations indicate that the magnetotail's cross-section is involved to a larger spatial extent than previously considered in the substorm process. Furthermore, convection changes result in topological changes of the magnetosphere on a fast time scale. The results show that the magnetosphere is, on a global magnetospheric scale, highly dynamic during convection changes and ensuing substorms.

Plasma physics governs the area of interactions between charged particles. As 99% of the visible universe is in a plasma state, it is an important topic in astronomy and space physics, where we already at an altitude of 60 km reach the plasma environment surrounding our planet in the form of the ionosphere. The search for fusion, the source of power for the sun, as well as industrial use have been the main topics for earth bound plasma reasurch. A plasma is composed of charged particles which interact by the electromagnetic force. In the kinetic description, via the Vlasov-Maxwell equations, the system is described in terms of probability distribution functions for each particle species, expressed in terms of particles position and velocity. The particles interact via self-consistent fields as determined by Maxwell's equations. For understanding the complex behaviour of the system, we need numerical solvers. These come in two flavours, Lagrangian methods, dealing with the moving around of synthetic particles, and Eulerian methods, which solve the set of partial differential, Vlasov and Maxwell equations. To perform the computations within reasonable time, we need to distribute our calculations on multiple machines, i.e. parallel programming, with the best possible matching between our computational needs and the need of splitting algorithms to adapt to our processing environment. Paper I studies electron and ion beams within a Lagrangian and fluid model and compare the results with experimental observations. This is continued with studies of a full kinetic system, using an Eulerian solver, for a closer look at electron-ion interactions in relation to ionospheric observations, (Papers II and IV). To improve the performance of the Eulerian solver it was parallelised (Paper III). The thesis is ending with the possibility to observe ultrahigh energy neutrinos from an orbiting satellite by using the Moon's surface as a detector Paper V.

The Cluster spacecraft have been used to investigate auroral electric fields and field-aligned currents (FACs) at geocentric distances between 4 and 7 Re. The electric fields have been measured by the EFW instrument, consisting of two pairs of spherical probes, and the FACs have been calculated from measurements of the magnetic field by the FGM fluxgate magnetometer. CIS ion and PEACE electron measurements have also been used. Event studies as well as statistical studies have been used to determine the characteristics of the auroral electric fields. In two events where regions of both spatial and temporal electric field variations could be identified, the quasi-static electric fields were, compared to the Alfvén waves, found to be more intense and contribute more to the downward Poynting flux. With the use of the four Cluster spacecraft, the quasi-static electric field structures were found to be relatively stable on the time scale of at least half a minute. Quasi-static electric fields were found throughout the altitude range covered by Cluster in the auroral region. The electric field structures were found both in the upward and downward current regions. Bipolar and monopolar electric fields, corresponding to U- and S-shaped potential structures, have been found at different plasma boundaries, consistent with the view that the plasma conditions and the geometry of the current system are related to the shape of the electric field. The type of the bipolar electric field structures (convergent or divergent) was further found to be consistent with the FAC direction. The typical scale sizes of the electric field structures have been determined to be between 4 and 5 km, when mapped to ionospheric altitude. The most intense FACs associated with intense electric fields were found for small FAC widths. The widths of upward and downward FACs were similar. / QC 20100730

A magnetospheric substorm is manifested in a variety of phenomena observed both in space and on the ground. Two electromagnetic signatures are the Substorm Chorus Event (SCE) and Pi2 pulsations. The SCE is a Very Low Frequency (VLF) radio phenomenon observed on the ground after the onset of the substorm expansion phase. It consists of a band of VLF chorus with rising upper and lower cutoff frequencies. These emissions are thought to result from Doppler-shifted cyclotron resonance between whistler mode waves and energetic electrons which drift into an observer’s field of view from an injection site around midnight. The ascending frequency of the emission envelope has been attributed to the combined effects of energy dispersion due to gradient and curvature drifts and the modification of the resonance conditions resulting from the radial component of the E × B drift. Two numerical models have been developed which simulate the production of a SCE. One accounts for both radial and azimuthal electron drifts but treats the wave-particle interaction in an approximate fashion, while the other retains only the azimuthal drift but rigorously calculates both the electron anisotropy and the wave growth rate. Results from the latter model indicate that the injected electron population should have an enhanced high-energy tail in order to produce a realistic SCE. Pi2 are damped Ultra Low Frequency (ULF) pulsations with periods between 40 and 150 s. The impulsive metamorphosis of the nightside inner magnetosphere during the onset of the substorm expansion phase is accompanied by a broad spectrum of magnetohydrodynamic (MHD) waves. Over a limited range of local times around midnight these waves excite field line resonances (FLRs) on field lines connected with the auroral zone. Compressional waves propagate into the inner magnetosphere, where they generate cavity mode resonances. The uniform frequency of Pi2 pulsations at middle and low latitudes is a consequence of these cavity modes. A number of Pi2 events were identified at times when the Cluster constellation was located in the nightside inner magnetosphere. Electric and magnetic field data from Cluster were used to establish the existence of both cavity and field line resonances during these events. The associated Poynting flux indicated negligible radial or field-aligned energy flow but an appreciable azimuthal flux directed away from midnight. / QC 20100920

<p>The turbulen flow in an asymmetric diffuser has been en studied by the means of Reynold average Navier-Stokes equations with both differential and explict algebraic expressions to model the Reynolds stress tensor. Modifications to the differential stress model have been derived, using the inverse turbulence timescale to obtain the dissipation of turbuence kinetic energy. The explicit algebraic Reynolds stress model has been used in combination with a two-equation platform to close the system of equations. Modifications made to the transport equation for the inverse turbulence timescale has made it possible to substantially relax the deman on near-wall resolution of this quantity. The rapid growth wth present in the original formulation can be treated as an explicit function of the wall-normal distance. In order to use the new formulation for the transport equation, an equation has as been derived to obtain the shortest distance bettwee a point and the closest wall, regardles of the geometric complexity of the domain. An explicit algebraic expression to model the passive scalar flux vector has been investigated using a comparison with a standard eddy-diffusivity model in the asymmetric diffuser. Results show a substantial improvement of the complexity of the scalar field and scalar flux vector in sepaarated flows. Automated code generation has been used in all the above studies to generate versatile model testing tools for general two-dimensional geometries. Finite element formulations are used for these tools.</p>