A numerical model was developed to understand the time evolution of a wake structure around a CubeSat moving in a plasma with transonic speed. A cubeSat operates in the F2 layer of ionosphere with an altitude of 300 − 600 Km. The average plasma density varies between 10−6cm−3 − 10−9cm−3 and the temperature of ions and electrons is found between 0.1−0.2 eV. The study of a wake structure can provide insights for its effects on the measurements obtained from space instruments. The CubeSat is modeled to have a metal surface, which is a realistic assumption, with a negative electric potential. To solve the equations of plasma, the numerical difference equations were obtained by discretizing the fluid equations of the plasma along with nonlinear Poisson’s equation. The electrons were assumed to follow the Boltzmann’s relation and the dynamics of ions was followed using the fluid equations. The initial and boundary conditions for the evolution of the structure are discussed. The computation was compared to the analytical solution for a 1D problem before being applied to the 2D model. There was a good agreement between the numerical and analytical solution. In the 2D simulation, we observe the formation of plasma wake structure around the CubeSat. The plasma wake structure consists of rarefaction region where ion density and ion velocity decreases compared to the initial density and velocity.

When plasmas come in contact with the boundaries that confine them, various complex processes occur between the plasma and the materials in the boundary. These processes, called plasma-material interactions (PMI) lead to physical and chemical modifications in the materials and in the plasma. In the case of a tokamak, a magnetic confinement fusion reactor, the interactions between the plasma and the material in the bounding walls can negatively impact the performance and service life of the reactor. Furthermore, PMI are also found in other areas of significant engineering interest, such as plasma-based spacecraft propulsion engines, where interactions affect the transport properties of the plasma and consequently the performance of the engine. Therefore, gaining a fundamental understanding of the various plasma-material interactions is necessary for the development and improvement of these devices. PMI are dictated by the plasma sheath, a layer of net positive charge that forms at the plasma-boundary interface. The sheath regulates the energy and particle fluxes to the boundary, mediating the interactions. Sheaths, however, are only stable and well-developed when the ions enter the sheath with a speed equal to or greater than the `Bohm speed'. The Bohm speed is a landmark result in sheath theory and various mathematical expressions for it have been derived from fluid and kinetic treatment of plasmas. Although these models are widely used, they are only accurate in cases where the thickness of the sheath is negligible when compared to the scale length of the plasma in consideration. These cases are said to satisfy the `asymptotic limit'. To resolve this, a new Bohm speed model that considers the effects of transport terms such as the electron heat flux, thermal force, and temperature isotropization has been recently proposed [Y. Li et al., Physical Review Letters (2022)]. The model is verified using particle-in-cell (PIC) kinetic simulations and is shown to accurately predict the Bohm speed in cases away from the asymptotic limit. This thesis investigates the new model using the continuum kinetic approach on the Gkeyll software framework. The continuum kinetic approach numerically solves the Vlasov-Maxwell equations using the discontinuous Galerkin method and captures the kinetic phenomena of the plasma without needing to track individual particles. Multiple collisional cases ranging from a Knudsen number of 20 to 5000 are considered in a 1X3V simulation domain using the Lenard-Bernstein collisional operator. The results of the continuum kinetic simulations are benchmarked to the PIC simulation results. It is concluded that across a wide range of collisionalities, the continuum kinetic method captures much of the same physics as the PIC method while offering noise-free results. However, there is a discrepancy between the Bohm speed prediction and the simulation results in the continuum kinetic case. This discrepancy is explored and significant error in the collisional integral derived transport terms between the continuum kinetic method and PIC method is found, suggesting that the difference in collisional operator may be the source of the discrepancy. Nevertheless, the sheath profiles developed in the PIC simulations and the continuum kinetic simulations are in reasonable agreement. / Master of Science / Nuclear fusion is a process in which two light atomic nuclei (like hydrogen) fuse to form a heavier nucleus (like helium) and release tremendous amounts of energy. The resultant energy from these reactions powers the sun and has the potential to provide clean energy for our terrestrial needs. Harnessing fusion energy has been one of the biggest scientific and engineering challenges of our time due to various reasons. One of these reasons is the interaction of plasma, which is the fuel for the fusion reaction, and the materials of the walls that bound the plasma. These interactions are called plasma-material interactions (PMI) and can affect the longevity and performance of fusion reactors. The main governing phenomenon behind these interactions is the plasma sheath, a layer of plasma that is formed when the plasma encounters a boundary. For a sheath to form it is also necessary that ions in the plasma possess a speed greater than the so-called `Bohm speed'. While many expressions have been derived for the Bohm speed, these expressions are only valid when there is a clear sheath entrance that divides the bulk plasma and the sheath. This condition is not satisfied in many cases of interest. Instead, a sheath-transition region is found to exist between the bulk plasma and the sheath. A recently proposed Bohm speed model [Y. Li et al., Physical Review Letters (2022)] resolves this. This model is accurate in cases where the sheath-transition region exists and is derived by considering previously overlooked transport physics. In this work, this new model is studied using a different computational approach known as `continuum kinetics' using an open-source solver called Gkeyll. The results of the continuum kinetic simulations are compared to the results used to verify the model.

No description available.

<p>A numerical simulation case involving space plasma and the evolution of instabilities that generates very fast electrons, i.e. approximately at half of the speed of light, is used as a test bed for scientific visualisation techniques. A visualisation system was developed to provide interactive real-time animation and visualisation of the simulation results. The work focuses on two themes and the integration of them. The first theme is the storage and management of the large data sets produced. The second theme deals with how the Visualisation System and Visual Objects are tailored to efficiently visualise the data at hand. </p><p>The integration of the themes has resulted in an interactive real-time animation and visualisation system which constitutes a very powerful tool for analysis and understanding of the plasma physics processes. The visualisations contained in this work have spawned many new possible research projects and provided insight into previously not fully understood plasma physics phenomena.</p>

A numerical simulation case involving space plasma and the evolution of instabilities that generates very fast electrons, i.e. approximately at half of the speed of light, is used as a test bed for scientific visualisation techniques. A visualisation system was developed to provide interactive real-time animation and visualisation of the simulation results. The work focuses on two themes and the integration of them. The first theme is the storage and management of the large data sets produced. The second theme deals with how the Visualisation System and Visual Objects are tailored to efficiently visualise the data at hand. The integration of the themes has resulted in an interactive real-time animation and visualisation system which constitutes a very powerful tool for analysis and understanding of the plasma physics processes. The visualisations contained in this work have spawned many new possible research projects and provided insight into previously not fully understood plasma physics phenomena.