University of Maryland

Parvez N. Guzdar

Institute for Plasma Research
University of Maryland
College Park, Maryland
20742-3511

A.V. Williams Bldg., Room 3317
Phone: (301) 405-1590
Email: guzdar@glue.umd.edu
FAX: (301) 405-1678

ResearchPublications

I am a Senior Research Scientist in the Institute for Plasma Research at the University of Maryland, College Park. I received my doctorate in physics from Gujarat University, India, in 1976. My current research interests are: (1) two-dimensional/three-dimensional nonlinear fluid simulations of fusion, space and magnetospheric plasmas, (2) two-dimensional simulations of Navier-Stokes fluids, (3) two-dimensional modeling of self-focusing of radio waves in high-latitude ionospheric plasmas, and (4) pattern formation in the two-dimensional complex Ginzburg-Landau equations.


Research

b Limit Disruptions in Tokamaks
Kleva, Guzdar

The most serious impediment to the practical utilization of tokamaks as fusion reactors is the limitation on the plasma thermal pressure imposed by disruptions. The ratio b of the thermal pressure to the pressure of the confining magnetic field provides a measure of the efficiency of a magnetic confinement fusion reactor. Operation at high b is very desirable because it yields a large fusion reaction rate relative to the cost of the confining magnetic field. However, experimental attempts to increase the ratio b beyond a critical limit bc have been thwarted by an abrupt, catastrophic loss of confinement. Not only do these disruptions limit b and, therefore, limit the efficiency of a tokamak, but the disruptions themselves can cause serious damage to the reactor. Displayed below are results from 3D MHD simulations performed on the T3E machine at NERSC.

The figure above shows three-dimensional isosurfaces of the pressure as the instability develops along ridges dominantly aligned along the ambient magnetic field.

The figure on the right is a poloidal projection of the pressure that shows the hot plasma fingers as they reach the wall. Our nonlinear simulations of tokamak stability reproduce the salient features of the disruptive loss of confinement. High b toroidal equilibria are linearly unstable to ballooning modes that grow on the pressure gradient on the large R side of the magnetic axis, where R is the major radius of the torus. The convection cells associated with the unstable ballooning modes transport the thermal energy toward the wall at large R in hot plasma ridges whose two-dimensional projection in the poloidal plane resembles fingers. As the hot central plasma is transported out in R in fingers, the average pressure gradient is reduced so that it no longer completely balances the Lorentz force of the confining magnetic field. As a result, there is a small net inward force in R caused by the unbalanced Lorentz force. At lower bthis force generates an axisymmetric flow that opposes the growth of the fingers to the wall at large R, thereby stabilizing the plasma nonlinearly and maintaining confinement. However, as b increases the growth rate of the plasma fingers towards the wall becomes so rapid that there is insufficient time for the self-consistently generated axisymmetric flow to halt their progress before they strike the wall, and confinement is lost.

This research is supported by the U.S. Department of Energy.


High Latitude Ionospheric Modification by Radio Waves
Gondarenko, Guzdar

The artificial generation of large and small-scale ionospheric irregularities using HF-heaters can impact a variety of space and ground-based communications, navigation, and surveillance systems. The understanding of the generation of these irregularities is thus a key objective. One mechanism that can create the irregularities is the thermal self-focusing of the high-powered heater waves. Over the last few years we have developed a 2D nonlinear code which investigates the spatio-temporal development of such irregularities. In the accompanying figure we show the development of the heater wave (left panels), the electron temperature (middle panels) and the plasma density (right panels) at t = 5.6s near the reflection height of a linearly stratified ionosphere. The heater wave generates the irregularities by the process of thermal self-focusing near the critical height and the irregularities, diffuse along the magnetic field lines (vertical direction) producing the field aligned structures. The goal of our studies is to investigate the spectrum of these irregularities for different ionospheric conditions and different polarizations and power density of the heater wave.

This research is supported by the National Science Foundation and the Office of Naval Research.


3D Simulation of the Dynamics of High Latitude Plasma Patches
Gondarenko, Guzdar

The ionospheric plasma at high latitude is known to display varied characteristics. During the phase when the interplanetary magnetic field is southward, large scale plasma patches/blobs, typically hundreds of kilometers, are formed in the cusp region and convect to the polar cap over periods of several hours. These patches are seen to have small scale irregularities, which are typically ten to a hundred times smaller than the size of the patch. We have developed a three-dimensional code for the plasma patch to gain an understanding of this meso-scale structuring. Earlier two-dimensional simulations showed that the structuring would lead to rapid fragmentation of the patch, contrary to observations. Introduction of the third dimension along the earth's field line causes the long wavelength instabilities to be stabilized. Thus this work will hopefully lead to a more complete understanding of the selection of the scalelengths of the structuring in the nonlinear phase and better agreement with observations.

The top figure shows five isosurfaces for the initial patch, with the peak density in red. After about an hour the patch develops irregularities (shown in the lower figure) due to the gradient drift instability and secondary Kelvin-Helmholtz instabilities.

This work is supported by the National Science Foundation.


Novel Algorithms
Guzdar

The study of the nonlinear dynamics of a higher dimensional system represented by one-dimensional and two-dimensional partial differential equations requires fast, efficient algorithms for generating long time-series data sets (for computing suitable thermodynamic averages). The algorithms remove severe time-step restrictions encountered by earlier workers and therefore can be solved on low-end computational platforms, which are readily accessible. One equation that has been solved using such techniques is the Kuramoto-Shivashinsky equation. The accompanying diagram shows the space-time evolution of a scalar field represented by the KS equation. Another equation is the 3D Ginzburg-Landau equation. The stability of 3D scrolls, rings, and spirals has been studied in collaboration with Prof. Edward Ott, Michael Gabbay, and Keeyeol Nam. More recently we have applied these techniques to develop an efficient algorithm for coupled nonlinear Schrodinger equations representing modulation and Raman scattering of laser beams propagating in an optical fiber. This work is in collaboration with Prof. Raj Roy and graduate student Bhaskar Khubchandani.

Publications since 1997

  1. "Nonlinear Stability Limit in High b Tokamaks," Robert G. Kleva and Parvez N. Guzdar, Phys. Plasmas 7, 1163 (2000).PDF

  2. "Diffraction Model of Ionospheric Irregularity-Induced Heater-Wave Pattern Detected on the WIND satellite," P. N. Guzdar, N. A. Gondarenko, K. Papadopoulos, G. M. Milikh, A. S. Sharma, P. Rodriguez, Yu. V. Tokarev, Yu. I. Belov, and S. L. Ossakow, Geophys. Res. Lett. 27, 317-320 (2000).Abstract

  3. "Lagrangian Chaos and the Effect of Drag on the Enstrophy Cascade in Two-Dimensional Turbulence," Keeyeol Nam, Edward Ott, Thomas M. Antonsen, Jr., and Parvez N. Guzdar, Phys. Rev. Lett. 84, 5134 (2000).PDF

  4. "Collisionless Nonideal Ballooning Modes," Robert G. Kleva and Parvez N. Guzdar, Phys. Plasmas 6, 116-121 (1999). PDF

  5. "Reply to Comment on b Limit Disruptions in Tokamaks," R. G. Kleva and P. N. Guzdar, Phys. Rev. Lett. 82, 5414 (1999).PDF

  6. "k Spectrum of Finite Lifetime Passive Scalars in Lagrangian Chaotic Fluid Flows," Keeyeol Nam, Thomas M. Antonsen, Jr., and Parvez N. Guzdar, Phys. Rev. Lett. 83, 3426 (1999).PDF

  7. "Cross-Field Transport Due to Low-Frequency Oscillations in the Auroral Region: A Three-Dimensional Simulations," Supriya B. Ganguli, Parvez N. Guzdar, Valeriy V. Gavrishchaka, Warren A. Krueger, and Paul E. Blanchard, J. Geophys. Res. 104, 4297 (1999).Abstract

  8. "Spatio-Temporal Development of Filaments Due to Thermal Self-Focusing Instability near the Critical Surface in the Ionosphere," N. A. Gondarenko, P. N. Guzdar, G. M. Milikh, A. S. Sharma, K. Papadopoulos, and S. L. Ossakow, Izv. VUZ Radiof. 42 (7), 670-681 (1999).Abstract

  9. "Excitation of Short-Scale Density Structures by Drift Waves during Ionospheric Heating," S. N. Antani and P. N. Guzdar, Geophys. Res. Lett. 26, 3285-3288 (1999).Abstract

  10. "Gradient Drift Instability in High Latitude Plasma Patches: Ion Inertial Effects," N. A. Gondarenko and P. N. Guzdar, Geophys. Res. Lett. 26, 3345-3348 (1999).Abstract

  11. "Three-Dimensional Simulations of the Ionospheric Plasma Transport in the Presence of the Structured Field-Aligned Flows," Valeriy V. Gavrishchaka, Supriya B. Ganguli, and Parvez N. Guzdar, J. Geophys. Res. 104, 22,511-22,524 (1999).Abstract

  12. "Three-Dimensional Nonlinear Simulations of the Gradient Drift Instability in a High-Latitude Ionosphere," P. N. Guzdar, N. A. Gondarenko, and P. K. Chaturvedi, Radio Sci. 33, 1901-1913 (1998). Abstract

  13. "The Thermal Self-Focusing Instability near the Critical Surface in the High-Latitude Ionosphere," P. N. Guzdar, P. K. Chaturvedi, K. Papadopoulos, and S. L. Ossakow, J. Geophys. Res. 103, 2231-2237 (1998). Abstract

  14. "b Limit Disruptions in Tokamaks," R. G. Kleva and P. N. Guzdar, Phys. Rev. Lett. 80, 3081 (1998). PDF

  15. "The Derivation of Equations for Fluctuations and Transport in Flux-Tube Geometries," J. J. Martinell, P. N. Guzdar, and A. B. Hassam, Phys. Plasmas 5, 1273 (1998). PDF

  16. "Route to Chaos for a Two-Dimensional Externally Driven Flow," R. Braun, F. Feudel, and P. Guzdar, Phys. Rev. E 58, 1927 (1998).PDF

  17. "The Dynamics of Scroll Wave Filaments in the Complex Ginzburg-Landau Equation," Michael Gabbay, Edward Ott, and Parvez N. Guzdar, Physica D 118, 371 (1998).Abstract

  18. "Reconnection of Vortex Filaments in the Complex Ginzburg-Landau Equation," Michael Gabbay, Edward Ott, and Parvez N. Guzdar, Phys. Rev. E 58, 2576 (1998).PDF

  19. "Stability of Spiral Wave Vortex Filaments with Phase Twists," Keeyeol Nam, Edward Ott, Parvez N. Guzdar, and Michael Gabbay, Phys. Rev. E 58, 2580 (1998).PDF

  20. "Turbulence and the Formation of Transport Barriers in Finite Plasma," B. N. Rogers, J. F. Drake, Y. T. Lau, P. N. Guzdar, A. B. Hassam, S. V. Novakovskii, and A. Zeiler, Plasma Phys. Controlled Nucl. Fusion Res., IAEA, Vienna, 1997.

  21. "Charging of Substrates Irradiated by Particle Beams," P. N. Guzdar, A. S. Sharma, and S. K. Guharay, Appl. Phys. Lett. 71, 3302 (1997).PDF

  22. "Motion of Scroll Wave Filaments in the Complex Ginzburg-Landau Equation," Michael Gabbay, Edward Ott, and Parvez N. Guzdar, Phys. Rev. Lett. 78, 2012 (1997).PDF


Last updated August 2000