Turbulence in the ionosphere is important to understand because it can negatively affect communication signals. This work examines different scenarios in the ionosphere in which turbulence may develop. The two main causes of turbulence considered in this work are the gradient drift instability (GDI) and the Kelvin-Helmholtz instability (KHI). The likelihood of the development of the GDI during the August 17, 2017 total solar eclipse is studied numerically. This analysis uses the ``Sami3 is Also a Model of the Ionosphere" (SAMI3) model to study the effect of the eclipse on the plasma density. The calculated GDI growth rates are small compared to how quickly the eclipse moves over the Earth. Therefore, the GDI is not expected to occur during the solar eclipse. A novel 2D electrostatic pseudo-spectral fluid model is developed to study the growth of these two instabilities and the problem of ionospheric turbulence in general. To focus on the ionospheric turbulence, a set of perturbed governing equations are derived. The model accurately captures the GDI growth rate in different limits; it is also benchmarked to the evolution of instability development in different collisional regimes of a plasma cloud. The newly developed model is used to study if the GDI is the cause of density irregularities observed in subauroral polarization streams (SAPS). Data from Global Positioning System (GPS) scintillations and the Super Dual Auroral Radar Network (SuperDARN) are used to examine the latitudinal density and velocity profiles of SAPS. It is found that the GDI is stabilized by velocity shear and therefore will only generate density irregularities in regions of low velocity shear. Furthermore, the density irregularities cannot extend through regions of large velocity shear. In certain cases, the turbulence cascade power laws match observation and theory. The transition between the KHI and the GDI is studied by understanding the effect of collisions. In low collisionality regimes, the KHI is the dominant instability. In high collisionality regimes, the GDI is the dominant instability. Using nominal ionospheric parameters, a prediction is provided that suggests that there exists an altitude in the upper textit{F} region ionosphere above which the turbulence is dominated by the KHI. / Doctor of Philosophy / In the modern day, all wireless communication signals use electromagnetic waves that propagate through the atmosphere. In the upper atmosphere, there exists a region called the ionosphere, which consists of plasma (a mixture of ions, electrons, and neutral particles). Because ions and electrons are charged particles, they interact with the electromagnetic communication signals. A better understanding of ionospheric turbulence will allow for aid in forecasting space weather as well as improve future communication equipment. Communication signals become distorted as they pass through turbulent regions of the ionosphere, which negatively affects the signal quality at the receiving end. For a tangible example, when Global Positioning System (GPS) signals pass through turbulent regions of the ionosphere, the resulting position estimate becomes worse. This work looks at two specific causes of ionospheric turbulence: the gradient drift instability (GDI) and the Kelvin-Helmholtz instability (KHI). Under the correct background conditions, these instabilities have the ability to generate ionospheric turbulence. To learn more about the GDI and the KHI, a novel simulation model is developed. The model uses a method of splitting the equations such that the focus is on just the development of the turbulence while considering spatially constant realistic background conditions. The model is shown to accurately represent results from previously studied problems in the ionosphere. This model is applied to an ionospheric phenomenon known as subauroral polarization streams (SAPS) to study the development of the GDI and the KHI. SAPS are regions of the ionosphere with large westward velocity that changes with latitude. The shape of the latitudinal velocity profile depends on many other factors in the ionosphere such as the geomagnetic conditions. It is found that for certain profiles, the GDI will form in SAPS with some of these examples matching observational data. At higher altitudes, the model predicts that the KHI will form instead. While the model is applied to just the development of the GDI and the KHI in this work, it is written in a general manner such that other causes of ionospheric turbulence can be easily studied in the future.

The complexity of numerical ocean circulation models requires careful checking with a variety of test problems. The purpose of this paper is to develop a test problem involving Rossby and gravity waves in a two-layer fluid in a channel. The goal is to compute very accurate solutions to this test problem. These solutions can then be used as a part of the checking process for numerical ocean circulation models. Here, Chebychev pseudo-spectral methods are used to solve the governing equations with a high degree of accuracy. Chebychev pseudo-spectral methods can be described in the following way: For a given function, find the polynomial interpolant at a particular non-uniform grid. The derivative of this polynomial serves as an approximation to the derivative of the original function. This approximation can then be inserted to differential equations to solve for approximate solutions. Here, the governing equations reduce to an eigenvalue problem with eigenvectors and eigenvalues corresponding to the spatial dependences of modal solutions and the frequencies of those solutions, respectively. The results of this method are checked in two ways. First, the solutions using the Chebychev pseudo-spectral methods are analyzed and are found to exhibit the properties known to belong to physical Rossby and gravity waves. Second, in the special case where the two-layer model degenerates to a one-layer system, some analytic solutions are known. When the numerical solutions are compared to the analytic solutions, they show an exponential rate of convergence. The conclusion is that the solutions computed using the Chebychev pseudo-spectral methods are highly accurate and could be used as a test problem to partially check numerical ocean circulation models. / Graduation date: 2012

The Poiseuille Flow is the flow of a viscous incompressible fluid in a channel between two infinite parallel plates. The behaviour of flow is properly described by the well-known Navier-Stokes Equations. The fact that Navier-Stokes equations are partial differential equations makes their solution difficult. They can rarely be solved in closed form. On the other hand, numerical techniques can be applied successfully to the well-posed partial differential equations. In the present study pseudo-spectral method is implemented to analyze the Poiseuille Flow. The pseudo-spectral method is a high-accuracy numerical modelling technique. It is an optimum choice for the Poiseuille flow analysis due to the flows simple geometry. The method makes use of Fourier Transform and by handling operations in the Fourier space reduces the difficulty in the solution. Fewer terms are required in a pseudo-spectral orthogonal expansion to achieve the same accuracy as a lower order method. Karhunen-Lo&egrave / ve (KL) decomposition is widely used in computational fluid dynamics to achieve reduced storage requirements or construction of relatively low-dimensional models. In this study the KL basis is extracted from the flow field obtained from the direct numerical simulation of the Poiseuille flow.

A prestudy indicated superior accuracy and convergence properties of apseudo-spectral method compared to a spline-based method implemented byCòrdoba et al. in 2005 when solving the α-patches problem. In this thesis wefurther investigate the numerical properties of the pseudo-spectral method and makeit more robust by implementing the Nonequispaced Fast Fourier Transform. Wepresent a more detailed overview and analysis of the pseudo-spectral method and theα-patches problem in general and conclude that the pseudo-spectral method issuperior in regards to accuracy in periodic settings.