Title page for ETD etd-03142009-040454

Type of Document Master's Thesis
Author Ismail, Atikah
URN etd-03142009-040454
Title Fourier spectral methods for numerical modeling of ionospheric processes
Degree Master of Science
Department Electrical Engineering
Advisory Committee
Advisor Name Title
Scales, Wayne A. Committee Chair
Besieris, Ioannis M. Committee Member
Claus, Richard O. Committee Member
  • Plasma dynamics
Date of Defense 1994-05-05
Availability restricted
Fourier spectral and pseudospectral methods are used in numerical modeling of ionospheric processes, namely macroscopic evolution of naturally and artificially created ionospheric density irregularities. The simulation model consists of two-dimensional electrostatic nonlinear fluid plasma equations that describe the plasma evolution. The spectral and pseudospectral methods are used to solve the spatial dependence of these self-consistent equations. They are chosen over the widely used finite difference and finite element techniques since spectral methods are straightforward to implement on nonlinear equations. They are at least as accurate as finite difference simulations. A potential equation solver is developed to solve the nonlinear potential equation iteratively. Time integration is accomplished using a combination of leapfrog and leapfrog-trapezoidal methods. A FORTRAN program is developed to implement the simulation model. All calculations are performed in the Fourier domain.

The simulation model is tested by considering three types of problems. This is accomplished by specifying an initial density (Pedersen conductivity) profile that represents slab model density, density enhancement (due to releases such as barium), or density depletion (due to late times effects of electron attachment material releases) in the presence of a neutral wind. The evolution of the irregularities is monitored and discussed. The simulation results agree with similar results obtained using finite difference methods. A comparison is made between the ionospheric depletion and enhancement problems. Our results show that, given the same parameters and perturbation level, the depletion profiles bifurcate much faster than that of the enhancement. We argue that this is due to the larger growth rate in the E X B interchange instability of the density depletion case.

  Filename       Size       Approximate Download Time (Hours:Minutes:Seconds) 
 28.8 Modem   56K Modem   ISDN (64 Kb)   ISDN (128 Kb)   Higher-speed Access 
[VT] LD5655.V855_1994.I863.pdf 5.32 Mb 00:24:38 00:12:40 00:11:05 00:05:32 00:00:28
[BTD] next to an author's name indicates that all files or directories associated with their ETD are accessible from the Virginia Tech campus network only.

Browse All Available ETDs by ( Author | Department )

dla home
etds imagebase journals news ereserve special collections
virgnia tech home contact dla university libraries

If you have questions or technical problems, please Contact DLA.