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.
XGC-6D is a geometric full-f Particle-In-Cell code on unstructured meshes including XGC-like unstructured mesh, tetrahedral mesh, etc. The code is particularly suitable to study boundary plasma with steep gradient and frequency higher than ion gyro-frequency. References: Wang Z, Qin H, Sturdevant B, Chang CS. Geometric electrostatic particle-in-cell algorithm on unstructured meshes. Journal of Plasma Physics. 2021;87(4):905870406.
The quasi-steady evolution of solar magnetic fields in response to gradual photospheric changes is considered, with particular attention given to the threshold of a sudden eruption in the solar atmosphere. The formal model of an evolving, force-free field dependent on two Cartesian coordinates is extended to a field which is not force free but in static equilibrium with plasma pressure and gravity. The basic physics is illustrated through the evolution of a loop-shaped electric current sheet enclosing a potential bipolar field with footpoints rooted in the photosphere. A free-boundary problem is posed and then solved for the equilibrium configuration of the current sheet in a hydrostatically supported isothermal atmosphere. As the footpoints move apart to spread a constant photospheric magnetic flux over a larger region, the equilibria available extend the field to increasing heights.
XGC-6D is a geometric full-f Particle-In-Cell code on unstructured meshes including XGC-like unstructured mesh, tetrahedral mesh, etc. The code is particularly suitable to study boundary plasma with steep gradient and frequency higher than ion gyro-frequency. References: Wang Z, Qin H, Sturdevant B, Chang CS. Geometric electrostatic particle-in-cell algorithm on unstructured meshes. Journal of Plasma Physics. 2021;87(4):905870406.
A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0,2π] 2 ×[0,1] to a toroidal domain in $\mathbb{R}$ 3 . This approach employs an action built on the squared Jacobian and radial length. Extensive testing is conducted on general toroidal boundaries using a global Fourier–Zernike basis via action minimisation. The results demonstrate successful coordinate construction capable of accurately describing strongly shaped toroidal domains. The coordinate construction is successfully applied to the computation of three-dimensional magnetohydrodynamic equilibria in the GVEC code where the use of traditional coordinate construction by interpolation from the boundary failed.
Here, an extended thermal particle boundary condition is devised to more efficiently and accurately model laser–plasma interactions in overdense plasmas. Particle-in-cell simulations of such interactions require many particles per cell, and a large region of background plasma is often necessary to correctly mimic a semi-infinite plasma and avoid electron refluxing from a truncated plasma. For long-pulse lasers of many picoseconds, such constraints can become prohibitively expensive. Here, an extended particle boundary condition (absorber) is designed that instantaneously stops and re-emits energetic particles streaming toward the simulation boundary over a defined region, allowing sufficient time and space for a suitably cool return current to develop in the background plasma. Tunable parameters of the absorber are explained, and simulations using the absorber with a 3-ps laser are shown to accurately reproduce those of a causally separated boundary while requiring only 20% the number of particles.
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.
Here, the physics of astrophysical jets can be divided into three regimes: (i) engine and launch (ii) propagation and collimation, (iii) dissipation and particle acceleration. Since astrophysical jets comprise a huge range of scales and phenomena, practicality dictates that most studies of jets intentionally or inadvertently focus on one of these regimes, and even therein, one body of work may be simply boundary condition for another. We first discuss long standing persistent mysteries that pertain the physics of each of these regimes, independent of the method used to study them. This discussion makes contact with frontiers of plasma astrophysics more generally. While observations theory, and simulations, and have long been the main tools of the trade, what about laboratory experiments? Jet related experiments have offered controlled studies of specific principles, physical processes, and benchmarks for numerical and theoretical calculations. We discuss what has been accomplished on these fronts. Although experiments have indeed helped us to understand certain processes, proof of principle concepts, and benchmarked codes, they have yet to solved an astrophysical jet mystery on their own. A challenge is that experimental tools used for jet-related experiments so far, are typically not machines originally designed for that purpose, or designed with specific astrophysical mysteries in mind. This presents an opportunity for a different way of thinking about the development of future platforms: start with the astrophysical mystery and build an experiment to address it.
Selected problems in magnetospheric plasma physics are critically reviewed. The discussion is restricted to questions that are 'global' in nature (i.e., involve the magnetosphere as a whole) and that are beyond the stage of systematic survey or isolated study requirements. Only low-energy particle aspects are discussed. The article focuses on the following subjects: (1) the effect of the interplanetary magnetic field on the topography, topology, and stability of the magnetospheric boundary; (2) solar-wind plasma entry into the magnetosphere; (3) plasma storage and release mechanisms in the magnetospheric tail; and (4) magnetic-field-aligned currents and magnetosphere-ionosphere interactions. A brief discussion of the prospects for the solution of these problems during and after the International Magnetospheric Study is given.
Numerical simulations were used to investigate the electron bunching induced by a laser pulse in its penetration into a plasma that has a nonuniform density profile at the vacuum – plasma boundary, which reaches a plateau. A mechanism was found to exist for this plasma which fosters a denser electron grouping in the generated bunch in comparison with the plasma with a sharp boundary. This mechanism is related to the nonstationarity of the wake wave produced by the laser pulse at the beginning of the plasma density plateau. (paper)
Many applications in magnetic confinement fusion require the efficient calculation of surface integrals with singular integrands. The singularity subtraction approaches typically used to handle such singularities are complicated to implement and low-order accurate. In contrast, we demonstrate that the Kapur–Rokhlin quadrature scheme is well-suited for the logarithmically singular integrals encountered for a toroidally axisymmetric confinement system, is easy to implement and is high-order accurate. As an illustration, we show how to apply this quadrature scheme for the efficient and accurate calculation of the normal component of the magnetic field due to the plasma current on the plasma boundary, via the virtual-casing principle.
The plasma potential measured by cylindrical and planar Langmuir probes has been shown to differ from the plasma potential measured by emissive probes in the neighborhood of the presheath near a negatively biased electrode immersed in a weakly collisional low temperature argon plasma. There are two principal results demonstrated in this paper. First, while it is well known that Langmuir probes cannot reliably measure plasma potentials inside of sheaths, results presented here demonstrate that the problem persists in presheaths, the quasineutral plasma bordering sheaths. It is known that emissive probes analyzed in the limit of zero emission accurately measure the plasma potential in the sheath. It is now clear that they are the only known electrostatic probe technique able to measure the plasma potential accurately throughout the presheath. Second, it is shown that the difference between potential measurements made by Langmuir probes and emissive probes in the body of the plasma, farther than a presheath distance from the boundary, is not proportional to T e , as has been previously claimed.
Based on He(2+) and H(-) ion composition data from the Plasma Composition Experiment on ISEE 1, a scenario is proposed for the solar wind penetration of the earth's magnetic tail, which does not require that the solar wind plasma be magnetized. While this study does not take issue with the notion that earth's magnetic field merges with the solar wind magnetic field on a regular basis, it focuses on certain aspects of interaction between the solar wind particles and the earth's field, e.g, the fact that the geomagnetic tail always has a plasma sheet, even during times when the physical signs of magnetic merging are weak or absent. It is argued that the solar plasma enters along slots between the tail lobes and the plasma sheet, even quite close to earth, convected inward along the plasma sheet boundary layer or adjacent to it, by the electric fringe field of the ever present low-latitude magnetopause boundary layer (LLBL). The required E x B drifts are produced by closing LLBL equipotential surfaces through the plasma sheet.
During this reporting period effort continued in the areas: (1) understanding the mechanisms responsible for substorm onset, and (2) application of a fundamental description of field-aligned currents and parallel electric fields to the plasma-sheet boundary layer.
Resistive wall implementation for tokamaks has been achieved by both a 'thin wall' approximation in which time-varying fields are matched across a boundary for all modes; and, by meshing the region from the first wall to the outer boundary and setting certain fields to mimic a resistive wall. Both implementations were tested on existing tokamaks.
This review deals with the propagation of waves, especially radio waves in the ionosphere. In the macroscopic electromagnetic theory, the mathematical structure of wave propagation problems depends entirely on the properties of the dielectric operator in a magnetically nonpermeable medium. These properties can be deduced from general discussions of symmetry and considerations of physical principles. When the medium is specifically the ionosphere, various physical phenomena may occur. Because of a large number of parameters, it is desirable to define a parameter space. A point in the parameter space corresponds to a specific plasma. The parameter space is subdivided into regions whose boundaries correspond to conditions of resonance and cutoff. As the point crosses these boundaries, the refractive index surface transforms continuously.
A physical model for a thruster discharge was developed, consisting of a spatially diverging plasma sustained electrically between a small ring cathode and a larger ring anode in a cylindrical chamber with an axial magnetic field. The associated boundary-value problem for the coupled partial differential equations with mixed boundary conditions, which describe the electric potential and the plasma velocity fields, was solved in closed form. By means of quantum-mechanical perturbation theory, a formula for the number S(E) of atoms sputtered on the average by an ion of energy E was derived from first principles. The boundary-value problem describing the diffusion of the sputtered atoms through the surrounding rarefied electron-ion plasma to the system surfaces of ion propulsion systems was formulated and treated analytically. It is shown that outer boundary-value problems of this type lead to a complex integral equation, which requires numerical resolution.
The paper is concerned with magnetosphere boundary phenomena observed by the Imp 7 magnetic field, plasma, and plasma wave instruments in 1972 and 1973. Boundary locations for a 15-month period are surveyed, and the different types of crossings are described. The spacecraft crosses the dawn and dusk boundaries near 25 earth radii downstream, and the physical processes at the Imp 7 magnetopause appear to be intermediate between those observed over the poles and those observed at the lunar orbit. The Imp 7 orbit also traverses a downstream region near where 'fireball' phenomena occur. Electromagnetic wave modes detected in the broad low-frequency channel of the wave instrument are analyzed, and the interpretation of data of this type suggests that the broad low-frequency channel is sensitive to oscillations in the lower hybrid resonance region of the spectrum.
The unusually high solar wind pressure and strongly southward IMF on May 4, 1998, pushed the magnetopause well into the geosynchronous orbit which exposed the POLAR satellite to the magnetosheath and solar wind. We use a gasdynamic convected field model to predict the magnetosheath quantities and then compare them with the in situ observations. The model prediction helps to reduce the uncertainty in the timing of the solar wind arrival time and provides a reference value for each physical parameter. It also helps to resolve the location of the satellite during strong magnetic fluctuations near the magnetopause. The plasma measurements from the TIDE instrument, in conjunction with the magnetometer measurements, indicate that there is a magnetospheric boundary layer during the event. There are also transient signatures near the magnetopause which may be caused by magnetospheric flux transfer events.