Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plasma blobs”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Low-state disks and low-beta disks

Stellar black hole candidates (BHCs) exhibit bimodal spectral states. We calculate nonthermal disk spectra, demonstrating that a large photon index (alpha (sub x) approximately 2-3) observed in the soft (high) state is due to a copious soft photon supply, whereas soft photon starvation leads to a smaller index (alpha (sub x) approximately 1.5-2) in the hard (low) state. Thus, the absence of the soft component flux in the low state cannot be due to obscuration. A possible disk configuration during the low state is discussed. We proposed that a low-state disk may be a low-beta disk in which magnetic pressure may exceed gas pressure becuase of the suppression of field escape by a strong shear. As a result, disk material will take the form of blobs constricted by mainly toroidal magnetic fields. Fields are dissipated mainly by occasional reconnection events with a huge energy release. This will account for large-amplitude, aperiodic X-ray variations (flickering) and high-energy radiation with small alpha(sub x) from hard state BHCs and possibly from active galactic nuclei. Further, we propose a hysteretic relation between the mass-flow rate and plasma-beta, a ratio of gas pressure to magnetic pressure, for the spectral evolution of transient BHCs. The disk is in the low-beta state in quiescence and early rise. The low-beta disk is optically thin and affected by advection. A hard-to-soft transition occurs before the peak luminosity, since there is no advection-dominated branch at higher luminosities. An optically thick, high-beta disk appears at small radii. In the decay phase of the light curve, the standard-type disk becomes effectively optically thin, when a soft-hard transition is triggered. High-beta plasmas in the main body shrink to form minute blobs, and low-beta coronal plasma fills interblob space.

Mineshige, Shin↗
DOE OSTI2021

Scattering of Radio Frequency Waves by Density Filaments

The edge region and the scrape-off layer of magnetically confined fusion devices, like tokamaks and stellarators, are replete with turbulent plasma that is a mixture of coherent, blob or filament like, structures [1] and incoherent fluctuations [2]. The variation in the density due to turbulence can be comparable to or greater than the ambient density [2]. As part of an overall effort to optimize the efficiency of operation, radio frequency (RF) waves are commonly used for heating fusion plasmas and, in tokamaks, for generating plasma current needed for confinement and controlling instabilities. The RF waves are excited by antenna structures that are placed near the wall of a fusion device. In order to deliver energy and momentum to charged particles in the core of fusion plasmas, RF waves have to propagate through the turbulent plasma. In present fusion devices, the scrap-off layer and the edge plasma region is of the order of a few centimeters. In reactor type devices, like ITER, this region is expected to be of the order of tens of centimeters. Since the efficiency of operation of a fusion reactor is of prime importance, it is imperative that we understand the effect of turbulence on RF waves. The fluctuations in density lead to changes in the plasma permittivity. As in conventional electrodynamics, the propagation of RF waves through different dielectric media is subject to reflection, refraction, and diffraction. In this paper, we summarize our theoretical and computational studies on the propagation of RF waves through filamentary structures present in the scrape-off layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

On how structures convey non-diffusive turbulence spreading

Abstract We report on comprehensive experimental studies of turbulence spreading in edge plasmas. These studies demonstrate the relation of turbulence spreading and entrainment to intermittent convective density fluctuation events or bursts (i.e. blobs and holes). The non-diffusive character of turbulence spreading is thus elucidated. The turbulence spreading velocity (or mean jet velocity) manifests a linear correlation with the skewness of density fluctuations, and increases with the auto-correlation time of density fluctuations. Turbulence spreading by positive density fluctuations is outward, while spreading by negative density fluctuations is inward. The degree of symmetry breaking between outward propagating blobs and inward propagating holes increases with the amplitude of density fluctuations. Thus, blob-hole asymmetry emerges as crucial to turbulence spreading. These results highlight the important role of intermittent convective events in conveying the spreading of turbulence, and constitute a fundamental challenge to existing diffusive models of spreading.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Interaction of radio frequency waves with cylindrical density filaments: scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas , vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci. , vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

Physics↗
DOE OSTI2022

Interaction of radio frequency waves with cylindrical density filaments -- scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas, vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci., vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Krash-ESL-Final report

During the period of the current Grant, the effort of the UCSD researchers was focused on: i) The application of the gyrokinetic code COGENT code to study the kinetic drift wave instability and its impact on blob dynamics, and ii) The modeling of edge plasma transport in collisional and semi-collisional regimes, and iii) Impurity transport in edge plasma and the incorporation of plasma-impurity collision operator (hereafter unlike particle collision operator) into COGENT. Over the years four UCSD graduate students (J. Angus, W. Lee, Y. Zhang, and A. Knyazev), funded under this Grant, contributed to the COGENT development/utilization and defended their PhD theses, which are based, partly, on their research related to the COGENT. Three of them (J. Angus, W. Lee, and Y. Zhang) have received “The Best MAE Graduate Student of the Year” awards. Whereas W. Lee after defending his PhD thesis returned back to his home country (South Korea), J. Angus, and Y. Zhang are continuing their research in plasma physics in the LLNL and LANL, whereas A. Knyazev currently is a postdoctoral fellow at Columbia University working on the theoretical aspects of advanced stellarators. The results obtained in our studies were disseminated in many papers published in first-class refereed scientific journals and presented on prestigious scientific meetings including the Annual APS DPP Meeting, Sherwood Fusion Theory Conference, IAEA Fusion Energy Conferences, International Conferences on Plasma-Surface Interactions in Controlled Fusion Devices, International Workshops on Edge Plasma Theory in Fusion Devices, etc. In what follows we highlight the major accomplishments of the UCSD Team in the main areas of our research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of neutral interactions on parallel transport and blob dynamics in gyrokinetic scrape-off layer simulations

The effect of neutral interactions on scrape-off layer (SOL) turbulence is investigated in a continuum gyrokinetic code that has been coupled to a continuum kinetic model of neutral transport. This extends the work of a previous paper, which compared two NSTX SOL simulations in simple helical geometry, one with neutrals and one without. The former included electron-impact ionization, charge exchange, and wall recycling. Here, the case with neutrals is compared to a gyrokinetic-only simulation that includes an effective ionization source to separate the effect of sourcing from charge exchange collisions. It is observed that sourcing accounts for many features of the simulated SOL with neutrals, including density and temperature magnitudes and reduced normalized density fluctuations, but differences persist. In particular, a flatter density profile results due to changes in parallel transport when neutral collisions are included, illustrating the importance of neutral drag on global plasma properties. An analysis of coherent turbulent structures, or blobs, in these simulations demonstrates the case with neutrals has slower and larger blobs. Here, a series of seeded blob simulations corroborates the blob velocity observation. In general, the blob motion does not contribute significantly to radial transport in these simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Physics of Edge-Core Coupling by Inward Turbulence Propagation

The dynamics of edge-core coupling is critically important to the optimization of magnetically confined fusion plasmas. Since early proposals, there has been persistent speculation that inward propagation of turbulence from the boundary is a possible means to energize the edge-core coupling region. However, the detailed mechanism of this process has remained a mystery until recent experiments observed that regular, intense gradient relaxation events generated blob-void pairs very close to the last closed flux surface. Blobs ($\tilde{n}$ >0) propagate outward and detach from the bulk plasma, while voids ($\tilde{n}$ <0) propagate inward, and so stir the core plasma. Here, in this work, we demonstrate that this heretofore ignored process of void emission can drive a broad turbulent layer of width ∼100 𝜌 𝑠 , for typical parameters. The mechanism is the Cherenkov emission of drift waves from inward-propagating voids. The model shows promise to resolve several questions surrounding the shortfall problem and the strong turbulence in the edge-core coupling region.

drift waves↗
DOE OSTI2026

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A three-point velocity estimation method for two-dimensional coarse-grained imaging data

Time delay and velocity estimation methods have been widely studied subjects in the context of signal processing, with applications in many different fields of physics. The velocity of waves or coherent fluctuation structures is commonly estimated as the distance between two measurement points divided by the time lag that maximizes the cross correlation function between the measured signals, but this is demonstrated to result in erroneous estimates for two spatial dimensions. We present an improved method to accurately estimate both components of the velocity vector, relying on three non-aligned measurement points. We introduce a stochastic process describing the fluctuations as a superposition of uncorrelated pulses moving in two dimensions. Using this model, we show that the three-point velocity estimation method, using time delays calculated through cross correlations, yields the exact velocity components when all pulses have the same velocity. The two- and three-point methods are tested on synthetic data generated from realizations of such processes for which the underlying velocity components are known. The results reveal the superiority of the three-point technique. Finally, we demonstrate the applicability of the velocity estimation on gas puff imaging data of strongly intermittent plasma fluctuations due to the radial motion of coherent, blob-like structures at the boundary of the Alcator C-Mod tokamak.

Materials Science↗
NASA NTRS2004

Magnetic Fields and Flows in Open Magnetic Structures

Open magnetic structures connect the solar surface to the heliosphere and are thus of great interest in solar-terrestrial physics. This talk is primarily an observational review of what is known about magnetic fields and particularly flows in such regions with special focus on coronal holes and origins of the fast solar wind. First evidence of the connection between these two features was seen in correlations of Skylab data with in situ measurements of the solar wind soon after the discovery of coronal holes, which are now known to emanate from unipolar magnetic regions at the photosphere. Subsequently many observations of have been made, ranging from oscillations in the underlying photosphere and chromosphere, to possible beginnings of the solar wind as observed by Doppler shifts in high chromospheric and transition-region lines, to coronagraphic time-lapse studies of outward-moving blobs of material which perhaps trace elements of solar-wind plasma. Some of the many unresolved and controversial issues regarding details of these observations and their association with the solar wind will be discussed.

Jones, Harrlson P.↗
NASA NTRS2012

Formation and Evolution of a Multi-Threaded Prominence

We investigate the process of formation and subsequent evolution of prominence plasma in a filament channel and its overlying arcade. We construct a three-dimensional time-dependent model of a filament-channel prominence suitable to be compared with observations. We combine this magnetic field structure with one-dimensional independent simulations of many flux tubes. The magnetic structure is a three-dimensional sheared double arcade, and the thermal non-equilibrium process governs the plasma evolution. We have found that the condensations in the corona can be divided into two populations: threads and blobs. Threads are massive condensations that linger in the field line dips. Blobs are ubiquitous small condensations that are produced throughout the filament and overlying arcade magnetic structure, and rapidly fall to the chromosphere. The total prominence mass is in agreement with observations. The threads are the principal contributors to the total mass, whereas the blob contribution is small. The motion of the threads is basically horizontal, while blobs move in all directions along the field. The peak velocities for both populations are comparable, but there is a weak tendency for the velocity to increase with the inclination, and the blobs with motion near vertical have the largest values of the velocity. We have generated synthetic images of the whole structure in an H proxy and in two EUV channels of the AIA instrument aboard SDO. These images show the plasma at cool, warm and hot temperatures. The theoretical differential emission measure of our system agrees very well with observations in the temperature range log T = 4.6-5.7. We conclude that the sheared-arcade magnetic structure and plasma dynamics fit well the abundant observational evidence.

Luna, M.↗
DOE OSTI2021

Millimeter-wave beam scattering and induced broadening by plasma turbulence in the TCV tokamak

The scattering of millimeter-wave beams from electron density fluctuations and the associated beam broadening are experimentally demonstrated. Using a dedicated setup, instantaneous deflection and (de-)focusing of the beam due to density blobs on the beam path are shown to agree with full-wave simulations. The detected time-averaged wave power transmitted through the turbulent plasma is reproduced by the radiative-transfer model implemented in the WKBeam code, which predicts a ~50% turbulence-induced broadening of the beam cross-section. Finally, the role of core turbulence for the considered geometry is highlighted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

On the thermodynamics of the plasma sheet

The present study reinvestigates the evidence for nonadiabatic transport in the quiet central plasma sheet, using AMPTE IRM data from the plasma sheet boundary layer and active times selected on the basis of large AE values. It is found that as the plasma is transported from the plasma sheet boundary layer into the central plasma sheet, both its temperature and its density (n) increase. The plasma obeys the relation p varies as n exp 4/3 for quiet times (AE is less than 100 nT) and p varies as n exp 5/3 for AE greater than 300 nT. The temperature in the quiet plasma sheet is usually less than 6 keV, and high-temperature values are more likely to be observed in what is defined as the active plasma sheet. It is suggested that the plasma sheet contains a mixture of high-entropy 'bubbles' and low-entropy 'blobs.' It is argued that these either merge or are lost from the tail before they are convected into the near-earth tail.

Baumjohann, W.↗
DOE OSTI2023

Strongly intermittent far scrape-off layer fluctuations in Alcator C-Mod plasmas close to the empirical discharge density limit

Abstract Intermittent plasma fluctuations in the boundary region of the Alcator C-Mod device were comprehensively investigated using data time-series from gas puff imaging and mirror Langmuir probe diagnostics. Fluctuations were sampled during stationary plasma conditions in ohmically heated, lower single null diverted configurations with scans in both line-averaged density and plasma current, with Greenwald density fractions up to 0.85. Utilizing a stochastic model, we describe the plasma fluctuations as a super-position of uncorrelated pulses, with large-amplitude events corresponding to blob-like filaments moving through the scrape-off layer (SOL). A deconvolution method is used to estimate the pulse arrival times and amplitudes. The analysis reveals a significant increase of pulse amplitudes and waiting times as the line-averaged density approaches the empirical discharge density limit. Broadened and flattened average radial profiles are thus accompanied by strongly intermittent and large-amplitude fluctuations. Although these filaments are arriving less frequently at high line-averaged densities, we show that there are significant increases in radial far-SOL particle and heat fluxes which will further enhance plasma–wall interactions. The stochastic model has been used as a framework for study of the scalings in the intermittency parameter, flux and mean amplitude and waiting times, and is being used to inform predictive capability for the effects of filamentary transport as a function of Greenwald fraction.

Physics↗
DOE OSTI2022

Cross-field and parallel dynamics of SOL filaments in TCV

Using recently installed scrape-off layer diagnostics on the tokamak à configuration variable, we characterise the poloidal and parallel properties of turbulent filaments. We access both attached and detached divertor conditions across a wide range of core densities (f G ϵ [0.09, 0.66]) in diverted L-mode plasma configurations. With a gas puff imaging (GPI) diagnostic at the outer midplane we observed filaments with a monotonic increase in radial velocity (from 390 m s –1 to 800 m s –1 ) and cross-field radii (from 8.5 mm to 13.4 mm) with increasing core density. Interpreting the filament behaviour in the context of the two-region model, we find that they populate the ideal-interchange regime (C i ) in discharges at very low densities, and the resistive X (RX)-point regime for all other discharges. The scaling of filament velocity versus size shows good agreement with this interpretation. These results are discussed and compared with previous probe-based measurements for similar conditions, which mostly placed filaments in TCV in the resistive ballooning (RB) regime. In addition, for the first time in TCV, the parallel filament extension is studied by magnetically aligning the GPI measurements at the outboard midplane with a reciprocating probe in the divertor. In agreement with the filaments being in the ideal-interchange and the RX-point regimes, they are found to extend beyond the X-point into the outer divertor leg.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A reduced-model (nSOLT) simulation of neutral recycling effects on plasma turbulence in the divertor region of MAST-U

The 2D scrape-off-layer turbulence code (nSOLT), which includes 1D kinetic neutral–plasma interactions, is applied to study effects of neutral recycling on plasma turbulence for parameters illustrative of the MAST-U divertor region. Neutral recycling is modeled by injecting a fraction of the parallel plasma flux to the divertor back into the simulation domain as a source of Franck–Condon-distributed neutrals. Stationary sources, concentrated at the magnetic separatrix, model plasma streaming into the divertor region from the upstream scrape-off-layer and sustain plasma turbulence absent neutral recycling. Starting from one such no-neutrals equilibrium, we initiate recycling in a numerical experiment designed to diagnose and identify the effects of various neutral–plasma interactions on the divertor plasma, divertor turbulence, and plasma exhaust. The onset of recycling triggers an initial burst of enhanced cross field plasma transport that is quelled by ionization cooling and charge–exchange (CX) friction, with growing neutral pressure, leading to a quiescent, turbulence-free state. Diagnosis of this transient burst reveals that (1) the sudden increase in plasma density due to ionization dominates the onset of the burst, (2) electron cooling due to ionization increases collisionality and disconnects blob filaments from the sheath, and (3) CX friction drives tripole polarization of a blob that can dominate the curvature-driven dipole polarization, leading to the stagnation of blob propagation and reduced radial turbulent transport. It is shown that CX friction is negligible compared to sheath physics in determining equilibrium mean flow shearing rates, for parameters considered herein (specifically a short connection length to the divertor target), while it can significantly reduce interchange-instability growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1997

Numerical Simulation of Rotation-Driven Plasma Transport In the Jovian Magnetosphere

A Jupiter version of the Rice Convection Model (RCM-J) was developed with support of an earlier NASA SR&T grant. The conversion from Earth to Jupiter included adding currents driven by centrifugal force, reversing the planetary magnetic field, and rescaling various parameters. A series of informative runs was carried out, all of them solving initial value problems. The simulations followed an initial plasma torus configuration as it fell apart by interchange instability. Some conclusions from the simulations were the following: 1. We confirmed that, for conventional values of the torus density and ionospheric conductance, the torus disintegrates by interchange instability on a time scale of approx. one day, which is 1-2 orders of magnitude shorter than the best estimates of the average residence time of plasma in the torus. 2. In the model, the instability could be slowed to an arbitrary degree by the addition of sufficient impounding energetic particles, as suggested earlier by Siscoe et al (1981). However, the observed energetic particles do not seem sufficient to guarantee impoundment (e.g., Mauk et al., 1996). 3. Whether inhibited by impoundment or not, the interchange was found to proceed by the formation of long fingers, which get thinner as they get longer. This picture differed dramatically from the conventional radial-diffusion picture (e.g., Siscoe and Summers (1981)), more superficially with the outward-moving-blob picture (Pontius and Hill, 1989). The obvious limitation of the original RCM-J was that it could not represent a plasma source. We could represent the decay of a pre-existing torus, but we could not represent the way ionization of material from Io continually replenishes the plasma. We consequently were precluded from studying a whole set of fundamental issues of torus theory, including whether the system can come to a steady state.

Wolf, Richard A.↗