Search NASA⌕ Search

SEARCH · Search NASA

Results for “Large Plasma Device”

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 127 records · Page 7

Oxidation of gallium arsenide in a plasma multipole device. Study of the MOS structures obtained

The oxygen plasma oxidation of GaAs was studied in order to obtain extremely high frequency responses with MOS devices. In the multipole system a homogeneous oxygen plasma of high density can easily be obtained in a large volume. This system is thus convenient for the study of plasma oxidation of GaAs. The electrical properties of the MOS diodes obtained in this way are controlled by interface states, located mostly in the upper half of the band gap where densities in the 10 to the 13th power/(sq cm) (eV) range can be estimated. Despite these interface states the possibility of fabricating MOSFET transistors working mostly in the depletion mode for a higher frequency cut-off still exists.

Gourrier, S.↗

Particle resonances in stellarators

Resonances of high energy particles in magnetic confinement devices due to electromagnetic instabilities can strongly modify the particle distribution, leading to a reduction in fusion power and even discharge termination and particle loss to the device walls through an avalanche. The existence of a mode particle resonance depends on the properties of the equilibrium and particle parameters, and their number, location, and density can vary with device design. Recently, the advent of more powerful computing capabilities and advanced theoretical understanding has led to the design of non-axisymmetric devices or stellarators, which could prove to be more advantageous than tokamaks. Stellarators have the advantage of being immune to major disruptions because of the very low plasma current. One of the problems shared by both types of devices is the existence of resonances in particle orbits, which can lead to large amplitude high frequency instabilities and subsequent induced particle loss. We examine the number of resonances, their location, and dependence on particle energy for some stellarator designs.

43 PARTICLE ACCELERATORS↗

Data-driven picosecond X-ray imaging for quantitative plasma-induced shock characterization

Imaging dynamic events, especially shockwave behavior, is key to advancing high-energy-density (HED) research. Recent advances in fourth- and fifth-generation X-ray light sources allow for high-resolution imaging of fast phenomena, but limited beam time necessitates maximizing data acquisition. We present a benchtop-scale pulsed plasma device submerged in liquid heptane, capable of generating dynamic events at rates exceeding 10 Hz, supporting the field’s data-driven goals by producing large, high-quality imaging datasets. Using X-ray phase contrast imaging (XPCI) at the Advanced Photon Source, we imaged weak shockwaves (Mach ~ 1.2) in heptane interacting with plasma-induced cavitation bubbles, causing deviation from Rankine-Hugoniot behavior; to our knowledge, this represents the first direct imaging of such interaction. Our quantitative analysis offers insight into weak shock phenomena and energy-focusing applications in pulsed plasmas. These results highlight the potential for large datasets to advance dynamic HED research at current light source facilities, and have implications for fields such as inertial confinement fusion, plasma-enhanced chemical processing, and biomedical applications.

36 MATERIALS SCIENCE↗

Line identification of boron and nitrogen emissions in EUV and VUV wavelength ranges in the impurity powder dropping experiments of LHD and its application to spectroscopic diagnostics

An impurity powder dropper was installed in the 21st campaign of the Large Helical Device experiment (Oct. 2019–Feb. 2020) under a collaboration between the National Institute for Fusion Science and the Princeton Plasma Physics Laboratory for the purposes of real-time wall conditioning and edge plasma control. In order to assess the effective injection of the impurity powders, spectroscopic diagnostics were applied to observe line emission from the injected impurity. Therefore, extreme-ultraviolet (EUV) and vacuum-ultraviolet (VUV) emission spectra were analyzed to summarize observable impurity lines with B and BN powder injection. Emission lines released from B and N ions were identified in the EUV wavelength range of 5–300 Å measured using two grazing incidence flat-field EUV spectrometers and in the VUV wavelength range of 300–2400 Å measured using three normal incidence 20 cm VUV spectrometers. BI–BV and NIII–NVII emission lines were identified in the discharges with the B and BN powder injection, respectively. Useful B and N emission lines which have large intensities and are isolated from other lines were successfully identified as follows: BI (1825.89, 1826.40) Å (blended), BII 1362.46 Å, BIII (677.00, 677.14, 677.16) Å (blended), BIV 60.31 Å, BV 48.59 Å, NIII (989.79, 991.51, 991.58) Å (blended), NIV 765.15 Å, NV (209.27, 209.31) Å (blended), NVI 1896.80 Å, and NVII 24.78 Å. Applications of the line identifications to the advanced spectroscopic diagnostics were demonstrated, such as the vertical profile measurements for the BV and NVII lines using a space-resolved EUV spectrometer and the ion temperature measurement for the BII line using a normal incidence 3 m VUV spectrometer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gasdynamic Mirror Fusion Propulsion Experiment

Nuclear fusion appears to be the most promising concept for producing extremely high specific impulse rocket engines. One particular fusion concept which seems to be particularly well suited for fusion propulsion applications is the gasdynamic mirror (GDM). This device would operate at much higher plasma densities and with much larger LD ratios than previous mirror machines. Several advantages accrue from such a design. First, the high LA:) ratio minimizes to a large extent certain magnetic curvature effects which lead to plasma instabilities causing a loss of plasma confinement. Second, the high plasma density will result in the plasma behaving much more Re a conventional fluid with a mean free path shorter than the length of the device. This characteristic helps reduce problems associated with "loss cone" microinstabilities. An experimental GDM device is currently being constructed at the NASA Marshall Space Flight Center to provide an initial assessment of the feasibility of this type of propulsion system. Initial experiments are expected to commence in the late fall of 2000.

Emrich, William J., Jr.↗

Tomographic Reconstruction of the Neutron Time-Energy Spectrum from a Dense Plasma Focus

A dense plasma focus (DPF) nuclear fusion device is an attractive pulsed neutron source in many applications, due to its relatively large neutron yield, produced in a short time duration. In order to design a DPF that generates neutrons within a specified time profile, or generates a neutron energy spectrum with specific properties, it is necessary to be able to characterize and model the results of the fusion process releasing the neutrons. The time-energy spectrum of fusion neutrons is an ideal quantity to use to validate multiphysics codes that simulate the pinch and fusion processes, because it is a quantity that requires the physics of each of the stages leading up to and ending in the fusion reaction to be simulated correctly. In particular, since DPF fusion neutrons are not monoenergetic -- and there can often be several fusion pinches creating neutrons -- a computer simulation matching high quality neutron spectrum measurements provides great confidence in the fidelity of the simulation. In order to make such a comparison, it is first necessary to have quality measurements from which to infer the spectrum. In this work we pose neutron spectroscopy as the classical tomographic inverse problem from neutron time of flight data at multiple distances, but enhanced by using an additional measurement of the time profile of the fusion pinch near its source and by using detector pairs set up in a geometry that allows for scatter background subtraction. The detector pairs enhance the quality of the time of flight measurements, and the additional constraint posed by the measured time profile allows for reconstructions discretized as finely as the time measurements and in energy as finely as 100 keV, without the problem being underdetermined. We present results from a Deuterium-fueled DPF at the U.S. Department of Energy's Nevada National Security Site and show that we can infer the time-energy spectrum from our measurements for both single and multi-pinch fusion reactions with equal fidelity.

Catenacci, Jared↗

Software For Advanced Large-scale Analysis Of Magnetic Confinement For Numerical Design, Engineering & Research (salamander)

As magnetic confinement fusion energy gains traction internationally to enable abundant energy production, designing components for fusion systems is a pressing challenge. During the planned lifetime of a fusion device, components evolve in extreme environments and must withstand large, repeated thermal loads and bombardment by 14 MeV neutrons, plasma ions, and neutral particles (deuterium, tritium, and helium), corrosive conditions, etc. All these physical processes take place simultaneously, interact in intricate ways, and impose important constraints that can affect performance. Experimental data is rare and costly to obtain, making design particularly challenging. Predictive computational frameworks must be an integral part of an accelerated and cost-effective design process by modeling fusion system performance in simulated environments. To better understand component degradation and operational impacts on their performance, the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER) is designed as an open-source, fully integrated, multiphysics, multiscale, NQA-1 compliant framework facilitating 3D, high-fidelity fusion system modeling. To that end, SALAMANDER is a MOOSE-based framework, and therefore leverages MOOSE upstream libraries such as PETSc and libMesh to deliver sophisticated finite element, finite volume, and nonlinear solver technology for fusion energy simulations. SALAMANDER couples MOOSE physics module capabilities—such as thermal hydraulics, heat conduction, Navier-Stokes, and thermomechanics—with tritium transport via TMAP8, neutronics via Cardinal, and nascent particle-in-cell capabilities. Direct simulation Monte Carlo methods will be used to address neutral transport near the walls. By coupling all these physics in an integrated application, SALAMANDER will enable high-fidelity modeling of irradiation levels and plasma exposure conditions of plasma facing components and their impact on heat and tritium distributions, as well as the resulting mechanical constraints experienced by the plasma facing components and performance of blanket systems. Furthermore, SALAMANDER will be particularly suited for engineering studies thanks to the stochastic tool module readily available in MOOSE, allowing for extended uncertainty quantification and risk analysis studies. It is also able to use computer-aided design (CAD) meshes to model complex geometries, which is indispensable for fusion systems. SALAMANDER therefore supports design, safety, engineering, and research projects for magnetic confinement fusion systems

Simon, Pierre-Clement [Idaho National Laboratory (↗

Investigation of radiated-power for low aspect ratio fusion plasmas

Investigating radiated power in fusion plasmas is of utmost importance to understand the effect of undesired impurities, such as metals, in present devices, or also desired impurities, such as noble gasses, to purposefully radiate a large fraction of power in future devices. These studies are especially important for high Z impurities which will play a crucial role in future generation fusion pilot plants. In this work, we have developed a power radiation analysis module, which is used to investigate 2D distributions of impurity densities and radiated power asymmetries caused due to both plasma rotation and the cooling rate dependence on temperature profile for the cases of experimental NSTX plasmas and designed scenarios for the spherical tokamak advanced reactor (STAR), both being low aspect ratio tokamaks. Two different atomic databases have been tested during this work to study their impact on the radiated power distribution, especially for high Z impurities. Also in PRAM, self-consistent calculations of two-dimensional electron, main ion, and impurity ion densities are carried out using one-dimensional input density, temperature, and rotation profiles. In the case of NSTX, discharges with high rotation of ∼ 170 km s, measured with charge exchange recombination spectroscopy, have been investigated. Rotation-induced charge separation, leading to an electrostatic potential, is calculated iteratively to a self-consistent solution while testing high Z impurities to observe any 2D asymmetry in the core radiated power due to centrifugal forces. The STAR design, being much larger (R = 4 m), is projected to have a much lower rotation, and is shown to have low rotation-induced asymmetries, on the order of ten percent or less, between the low field and high field sides. However, another effect not due to rotation but to the dependence of impurity cooling rates on temperature can lead to radiation peaking off-axis, near the plasma edge. This effect is noticeable for argon in NSTX, for example, but can also be enhanced for certain impurities at much higher temperatures projected for STAR (T e0 ~ 32 keV), for example for undesired tungsten or possibly desired xenon.

NSTX↗

Investigation of Recombination Processes In A Magnetized Plasma

Interplanetary travel requires propulsion systems that can provide high specific impulse (Isp), while also having sufficient thrust to rapidly accelerate large payloads. One such propulsion system is the Variable Specific Impulse Magneto-plasma Rocket (VASIMR), which creates, heats, and exhausts plasma to provide variable thrust and Isp, optimally meeting the mission requirements. A large fraction of the energy to create the plasma is frozen in the exhaust in the form of ionization energy. This loss mechanism is common to all electromagnetic plasma thrusters and has an impact on their efficiency. When the device operates at high Isp, where the exhaust kinetic energy is high compared to the ionization energy, the frozen flow component is of little consequence; however, at low Isp, the effect of the frozen flow may be important. If some of this energy could be recovered through recombination processes, and re-injected as neutral kinetic energy, the efficiency of VASIMR, in its low Isp/high thrust mode may be improved. In this operating regime, the ionization energy is a large portion of the total plasma energy. An experiment is being conducted to investigate the possibility of recovering some of the energy used to create the plasma. This presentation will cover the progress and status of the experiment involving surface recombination of the plasma.

Chavers, Greg↗

Generation of Large-Volume High-Pressure Plasma by Spatio-Temporal Control of Space Charge

Due to the fundamental limitations of scaling, very little progress has been made towards achieving a large-volume dense non-equilibrium plasma near atmospheric pressures. Commercially available state-of-the art glow-like atmospheric plasma devices for industrial applications are either narrow tubes or wide slits with narrow openings. Often multiple sources are put together in various configurations to process larger surfaces. The traditional methods of exciting electrodes create spatially fixed electric fields. As a result, the space charge at atmospheric pressure tends to be spatially confined resulting in non-uniformity which eventually leads to instability as the discharge is scaled. Theoretical work done under this project showed that it is possible to generate a spatially rotating electric field by exciting a set of electrodes with phase staggered sinusoidal waveforms. The modeling and simulations were done using plasma fluid models. It was shown that such a field can produce a uniform plasma. At the conclusion of the project, experimental proof of concept with an eight-electrode system in various gases (Air, Helium and Argon) was demonstrated. Power measurements and spectral investigation show that the concept can be used to generate a uniform stable plasma. This plasma source has the potential of opening new applications of nonthermal plasma including combustion of carbon-free fuel. The current limitation of the proposed method in scaling to higher volume and pressure is the need for multiple high voltage amplifiers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dissociation phenomena in electron-beam sustained carbon dioxide lasers

A number of applications are emerging requiring efficient, long pulse, long-life sealed CO2 lasers. Examples include the proposed NASA and ESA wind lidars. Electron-beam sustained discharge devices are strong contenders. Unlike self-sustained discharges, e-beam sustenance readily provides efficient performance from large volume discharges and offers pulse lengths well in excess of the microsecond or so generally associated with self-sustained devices. In the case of the e-beam sustained laser, since the plasma is externally maintained and operated at electric field strengths less than that associated with the glow to arc transition, the discharges can be run even in the presence of strongly attacking species such as O2. Build up of large levels of attacking contaminants is nevertheless undesirable as their presence reduces the current drawn by the plasma and thus the pumping rate to the upper laser level. The impedance rise leads to a mismatch of the pulse forming network with a consequent loss of control over energy deposition, operating E/N, and gain. Clearly CO2 dissociation rates, the influence of dissociation products on the discharge and gain, and tolerance of the discharge to these products need to be determined. This information can then be used to assess co-oxidation catalyst requirements for sealed operation.

Harris, Michael R.↗

Machine learning surrogates for ion energy–angle distributions in thermal and RF plasma sheaths

Ion energy–angle distributions (IEADs) at material surfaces are a critical input for plasma–material interaction (PMI) studies in fusion devices, yet they are computationally expensive to obtain using particle-in-cell (PIC) simulations. In this work, we develop a machine learning surrogate based on a deep deconvolutional neural network (DDeCNN) trained on large databases generated with the hPIC2 code. The surrogate is capable of reconstructing IEADs from sheath parameters for both thermal and radio-frequency (RF) plasmas, including cases with multiple ion species. Across thousands of test cases, the model achieves high accuracy, with over 97 % of predictions classified as good or average based on standard error metrics (MAE, MSE, L2). Even in the more challenging RF and multi-species regimes, the surrogate reliably captures the multi-peak structure of PIC results. Once trained, the surrogate produces IEADs in milliseconds on a common workstation, yielding speedups of six to seven orders of magnitude compared with running a full PIC simulation. This computational gain enables dense parameter scans and direct coupling of IEAD predictions with PMI and erosion models on whole-device scales in fusion-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗

A scoping study of far-SOL main-wall protection limiters for steady-state operation of compact pilot plant tokamaks

We present a novel method for handling steady-state heat fluxes incident on the main wall of pilot plant-scale magnetic fusion devices, based on the utilization of protection limiters in the far scrape-off layer (SOL). This method helps avoid large plasma-wall gaps, without excessively compromising blanket performance. We present an optimization algorithm for determining the appropriate size and scale of these protection limiters given (1) probability distributions of SOL plasma parameters and (2) assumed risk tolerance. As part of this optimization, we have developed an analytic description of parallel heat fluxes across limiter shadows, and an objective cost function (the ‘Far-SOL Marginal Cost’) to quantify the impact that different main-wall thermal management design choices have on reactor capital cost. Applying the model to a midscale fusion pilot plant concept shows that making use of far-SOL protection limiters can reduce capital costs on the order of $500 M, relative to naively increasing the plasma-wall gap. Our analysis demonstrates that the far-SOL power decay length is the highest-leverage plasma assumption for thermal loading of the first wall, and the primary cost driver for main wall thermal management. The relative cost efficiency of protection limiters increases as assumptions on the far-SOL heat flux become more pessimistic. The concepts described in this paper motivate the further development of far-SOL protection limiters as part of larger efforts to design economical core-edge-wall compatible solutions for a fusion pilot plant.

Design under uncertainty↗

Synthesis of Carbon Nanotubes Array by CVD

Since the properties of multi-wall carbon nanotubes (MWCNT) are superior in many devices such as electronics and sensors, many efforts have been involved in synthesizing particular structural or dimensional MWCNT. Uniform aligned MWCNT array is one of the prototype structures for devices such as filed emission device and microelectromechanical systems in which a large length to diameter ratio may also be required. Most aligned MWCNT recently synthesized by plasma enhanced chemical vapor deposition (CVD) have cone shaped structures. This presentation will illustrate aligned MWCNT array synthesized on silicon substrates using thermal CVD that could produce MWCNT with uniform diameter. An array of nickel particles was used as catalyst for MWCNT growth. A thin Ti or Au buffer layer was coated on the substrate prior to depositing nickel particles. Because the MWCNT size depends on the catalyst particle size, the nickel particle size annealed at various temperatures was investigated. MWCNT were grown on the substrate in the temperature range of 700 C - 1000 C and the pressure range of 1 to 300 torr. Methane and hydrogen gases with methane content of 1 - 10 % were used for the MWCNT synthesis. Morphology, length and diameter of MWCNT were determined by scanning electron microscopy and Raman spectroscopy. The detailed results of synthesis and characterizations will be discussed in the presentation.

Zhu, S.↗

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗

MHD Turbulence and Magnetic Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V↗

Temporal study of the current collected by a Langmuir probe biased to high voltages in the wake behind a conducting body

The temporal evolution of the voltage-current characteristic of a spherical Langmuir probe in the near-wake region of a conducting body is investigated in a pulsed plasma device. Measurements of the current collected by a sphere placed on axis behind a conducting disk are presented as a function of the bias voltage and time. The bias voltage on the collecting sphere is large enough, ranging from +500 to -3000 V, to attract ions or electrons across the wake region. In addition, plasma potential profiles as a function of sphere bias voltage and time throughout the nearwake and midwake region are presented. It is found that for a negative sphere bias a threshold voltage exists for current collection. This threshold voltage decreases during the formation of the wake behind the conducting disk. Potential measurements indicate that the sheath region around the negatively biased sphere contracts as the wake region establishes itself. This contraction of the sheath region increases the electric field attracting ions to the sphere and decreases the threshold voltage for current collection. For positive sphere bias voltage the current collected by the sphere does not exhibit a threshold voltage but is monotonically increasing.

Meassick, S.↗