TIME-DEPENDENT SOLPS-ITER SIMULATIONS OF THE TOKAMAK PLASMA BOUNDARY FOR MODEL PREDICTIVE CONTROL
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Self-consistent electric potential around axisymmetric body in stationary flow of rarefied plasma
The basic model of Lindhard and Scharff, known as the local plasma model, is used to study the effects on stopping power of the chemical and physical state of the medium. Unlike previous work with the local plasma model, in which individual electron shifts in the plasma frequency were estimated empirically, he Pines correction derived for a degenerate Fermi gas is shown herein to provide a reasonable estimate, even on the atomic scale. Thus, the model is moved to a complete theoretical base requiring no empirical adjustments, as characteristic of past applications. The principal remaining error is in the overestimation of the low-energy absorption properties that are characteristic of the plasma model in the region of the atomic discrete spectrum, although higher-energy phenomena are accurately represented, and even excitation-to-ionization ratios are given to fair accuracy. Mean excitation energies for covalent-bonded gases and solids, for ionic gases and crystals, and for metals are calculated using first-order models of the bonded states.
A time-dependent, zero-dimensional collisional-radiative model has been developed to study capillary discharge plasmas in view of possibilities for obtaining population inversion leading to laser action in the extreme ultraviolet and soft x-ray spectral region, by amplifying the Balmer-alpha (H-alpha) line of different hydrogenlike ions with nuclear charges from Z = 3 to Z = 6. The model is described in detail, and results for the case of a carbon plasma are presented. Limitations of the model are discussed and comparison is made between our calculations and recent experimental results.
In this work, we explore surrogate models to optimize plasma enhanced atomic layer deposition (PEALD) in high aspect ratio features. In plasma-based processes such as PEALD and atomic layer etching (ALE), surface recombination can dominate the reactivity of plasma species with the surface, which can lead to unfeasibly long exposure times to achieve full conformality inside nanostructures like high aspect ratio vias. Using a synthetic dataset based on simulations of PEALD, we train artificial neural networks to predict saturation times based on cross section thickness data obtained for partially coated conditions. The results obtained show that just two experiments in undersaturated conditions contain enough information to predict saturation times within 10% of the ground truth. A surrogate model trained to determine whether surface recombination dominates the plasma–surface interactions in a PEALD process achieves 99% accuracy. This demonstrates that machine learning can provide a new pathway to accelerate the optimization of PEALD processes in areas such as microelectronics. Our approach can be easily extended to ALE and more complex structures.
Detachment of the divertor plasma during application of resonant magnetic perturbation (RMP) fields is evaluated for hydrogen H-mode plasma during the first pre-fusion power operation (PFPO-1) phase in ITER by 3D plasma boundary modelling with EMC3–EIRENE. Plasma response effects from a linearized, resistive, single fluid MHD model are discussed, which includes partial screening of the externally applied field—but also field amplification near the separatrix. Furthermore, this field amplification is found to play a pivotal role for the magnetic footprint on the divertor targets, but is sensitive to model parameters. Extensions of the footprint beyond the straight portions of the ITER vertical divertor targets, optimized for high stationary heat flux handling, may be possible depending on the level of toroidal rotation in the plasma. Exhaust from the bulk plasma is guided by the helical corrugations (lobes) of the perturbed separatrix, and this results in an upstream heat flux that is distributed over these lobes with lower peak values than in the typical radial heat flux profiles seen in the absence of magnetic perturbations. As a consequence, an earlier onset (with respect to the upstream density) of detachment is found in the traditional strike zone when RMPs are applied, but secondary, non-axisymmetric strike locations appear—and those remain attached at temperatures above 10 eV. Neon seeding can mitigate these non-axisymmetric heat loads, but this becomes less efficient for large magnetic footprints.
A model of a plasma focus is examined wherein large axial electric fields are produced by an imploding current sheet during the final nanoseconds of the collapse phase and where the fields provide a mechanism for creating a beam of electrons of highly suprathermal energies. The expected bremsstrahlung radiation above 100 keV is calculated for such a beam, which has a power-law spectrum, both from electron-deuteron collisions in the focused plasma and when the beam reaches the wall of the device. It is concluded that, since the experimental results indicate little or no radiation above 100 keV originating in the walls, that the electrons in the beam must be decelerated after leaving the plasma and before reaching the wall. Comparisons with the results and the total energy of the device yield qualitative agreement with the expected angular distribution of hard X-rays and reasonable agreement with the total energy in accelerated electrons required to produce the observed total energy in hard X-rays by this mechanism.
Effects of dynamic wall outgassing on divertor plasma during an ELM-like heat pulse are modelled using newly coupled edge plasma code UEDGE and all reaction-diffusion code FACE. Different divertor regimes are simulated depending on the degree of the plasma detachment. It is shown that hydrogen outgassing from the divertor target plate induced by ELM pulse in the semi-detached regime can trigger deep plasma detachment. The induced detachment does not spontaneously reverse after the pulse due to ineffective implantation of low energy plasma hydrogen in the target material. Furthermore, we demonstrate that plasma neutrals play essential role in heat transport in a detached divertor. The modelled hydrogen dynamics in the wall material demonstrates that the heat pulse leads to enhanced hydrogen transport toward the wall bulk.
Interplanetary magnetic field and plasma structure observed with Mariners magnetometer, noting 27-day deviations from Parker spiral model
The reported investigation presents the magnetic observations averaged in a format that facilitates comparison with models of plasma diamagnetic effects. These models demonstrate that a decreased plasma energy density in the lunar vicinity always causes an induced dipole that increases the external field. Hence, an increased plasma energy density in the lunar environs is required to account for the observed magnetic decrease. There are several different plasma populations in the high-latitude geomagnetic tail lobes, including the lunar ions resulting from photoionization of lunar atmosphere neutrals, possible fluxes of polar-wind protons, low-energy particle events which consist of antisolar directed fluxes of 100-eV protons, photoelectrons from the dayside lunar surface, and a high-energy thermal population.
An analytical model is developed for the diffusive equilibrium plasma density distribution in the Io plasma torus. The model has been employed successfully to follow the ray path of plasma waves in the multi-ion Jovian magnetosphere; it would also be valuable for other studies of the Io torus that require a smooth and continuous description of the plasma density and its gradients. Validity of the analytical treatment requires that the temperature of thermal electrons be much lower than the ion temperature and that superthermal electrons be much less abundant than the thermal electrons; these two conditions are satisfied in the warm outer region of the Io torus from L = 6 to L = 10. The analytical solutions agree well with exact numerical calculations for the most dense portion of the Io torus within 30 deg of the equator.
A model in which thermal effects are simulated through use of a multibeam plasma distribution function is developed and investigated to see if solutions which take an initially uniform magnetized plasma to a new uniform state with different field orientation are possible. The momentum conservation integrals are found to admit two classes of such solutions, but only one class exhibits appropriate asymptotic behavior. Extensive numerical integrations have failed to demonstrate the existence of the desired solutions.
Plasmas in contact with liquids can degrade organic molecules in solution, as reactive oxygen and nitrogen species produced in the plasma solvate into the liquid. Immersing small droplets (tens of microns in diameter) in the plasma can more rapidly activate the liquid compared to treating a large volume of liquid with a smaller surface-to-volume ratio. The interactions between a radio frequency glow discharge sustained in He/H 2 O and a water droplet containing formate (HCOO - aq ) immersed in and flowing through the plasma were modeled using a 0-dimensional global plasma chemistry model to investigate these activation processes. HCOO - aq interacts with OH aq , produced from solvation of OH from the gas phase. HCOO - aq concentrations were benchmarked with previously reported experimental measurements. The diameter of the droplet, initial HCOO - aq concentration, and gas flow rate affect only the HCOO - aq concentration and OH aq density, leaving the OH density in the gas phase unaffected. Power deposition and gas mixture (e.g., percentage of H 2 O) change both the gas and liquid phase chemistry. A general trend was observed – during the first portion of droplet exposure to the plasma, OH aq primarily consumes HCOO - aq . However, O2 - aq , a byproduct of HCOO - aq consumption, consumes OH aq once O2 - aq reaches a critically large density. Using HCOO - aq as a surrogate for OH aq sensitive contaminants, combinations of residence time, droplet diameter, water vapor density, and power will determine the optimum remediation strategy.
A model of current collection in Langmuir probes is used to investigate the effect representative non- equilibrium plasmas under various conditions have on the electron temperature and number density that would be calculated through analysis of the probe collection characteristics. The model uses the distribution function to calculate the charged particle flux to a probe and then, for fixed applied voltages, the current values that satisfy continuity in the probes are determined. The triple probe is not scanned in voltage, so there is no practical way using experimental triple probe data to determine if the plasma is in equilibrium. As a consequence, a triple probe analysis typically relies on the assumption that the plasma is in equilibrium. Proceeding from this point, the numerically-generated non-equilibrium triple probe data are analyzed assuming that the plasma is in equilibrium, with the data compared to the initial distribution function inputs of plasma temperature and number density to determine the effect the non-equilibrium distribution has on plasma measurements. The temperature and number density are both significantly affected when a fraction of the particles in the distribution are shifted from the equilibrium configuration into the non-equilibrium part of the distribution function. For all instances studied, the computed electron temperature and number density are extremely sensitive to small deviations from equilibrium ( ≤ 5% of the plasma shifted into the non-equilibrium function). Shifting more of the plasma into the non-equilibrium distribution beyond this initial level does not produce a significant additional shift in the computed plasma properties.
Dusty plasmas are present in comets, in the ring systems of the outer planets, and in the interstellar medium. A two-dimensional particle-in-cell (PIC) model of a dusty plasma is presented in this paper. The PIC code is best suited for modeling the plasma-dust interaction for large grains, with diameters of the order of a centimeter. We have modeled the charging process for an individual dust grain and the associated potential pattern in the surrounding plasma. We have also considered the case of a large number of grains in a plasma, with intergrain separations of the order of the Debye length, and have shown that the plasma becomes depleted and the charge on a dust grain is reduced, as other workers in this field have predicted (cf. C. K. Goertz, 1989). We examine the electron and ion distribution functions in the vicinity of a charged grain and demonstrate that the ions near a grain have clearly been accelerated by the electrostatic potential.
The effect or hot magnetospheric electrons on the polar-plasma outflow was investigated, using a semikinetic model developed by Wilson et al. (1990) and Ho et al. (1991) to simulate the effect. The model is based on a hybrid particle-in-cell approach, in which the H(+) and O(+) ions are treated as adiabatic parallel-drifting gyrocenters injected as the upgoing portions of drifting bi-Maxwellian distributions at 1.6 R(E), while the electrons are treated as a massless neutralizing fluid. The results show that, in order to simulate the polar outflow under the influence of hot magnetospheric electrons, it is necessary to consider the effect of the electron temperature gradient.
A model is investigated for the decay of flare heated coronal loops in which rapid radiative cooling at the loop base creates strong pressure gradients which, in turn, generate large (supersonic) downward flows. The important features of this model which distinguish it from previous models of flare cooling are: (1) Most of the thermal energy of the coronal plasma may be lost by mass motion rather than by conduction or coronal radiation. (2) Flare loops are not isobaric during their decay phase, and large downward velocities are present near the footpoints. (3) The differential emission measure has a strong temperature dependence. These results can account for recent observations of compact flare loops that are not consistent with the previous cooling models.
Nonaxisymmetric magnetic fields such as the intrinsic error field and the externally applied resonant magnetic perturbation (RMP) in a tokamak are known to influence the plasma momentum transport and flow evolution through plasma response, which itself strongly depends on the plasma flow as well. The nonlinear interaction between plasma response and flow has been previously modeled in the conventional error field theory with the “no-slip” condition, which has been recently extended to allow the “free-slip” condition. In this work, we further target this specific process and numerically simulate the nonlinear plasma response and flow evolution in the presence of a single-helicity RMP in a circular-shaped model tokamak configuration, based on the full resistive MHD model in the initial-value code NIMROD. Time evolution of the parallel (to k) flow or “slip frequency” profile and its asymptotic steady state obtained from the NIMROD simulations are compared with both conventional and extended nonlinear response theories. Here, k is the wave vector of the propagating island. Good agreement with the extended theory with free-slip condition has been achieved for the parallel flow profile evolution in response to RMP in all resistive regimes, whereas the difference from the conventional theory with the no-slip condition tends to diminish as the plasma resistivity approaches zero.