Solar wind simulation with the Pioneer spacecraft.
Plasma wind tunnel simulation of solar wind effect on Pioneer spacecraft
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.
Plasma wind tunnel simulation of solar wind effect on Pioneer spacecraft
How the addition of ions at rest (mass loading) affects the temperature of a flowing plasma in a MHD approximation is investigated, using analytic theory and time dependent, three-dimensional MHD simulations of plasma flow past Io. The MHD equations show that the temperature can increase or decrease relative to the background, depending on the local sonic Mach number M(S), of the flow. For flows with M(S) of greater than sq rt 9/5 (when gamma = 5/3), mass loading increases the plasma temperature. However, the simulations show a nonlinear response to the addition of mass. If the mass loading rate is large enough, the temperature increase may be smaller than expected, or the temperature may actually decrease, because a large mass loading rate slows the flow and decreases the thermal energy of the newly created plasma.
A computer program is being developed to simulate the interaction of a plasma with a conducting disk. Two configurations are examined: (1) the conductor is a 'button' in the center of a larger dielectric disk, and (2) the conducting disk is covered by a dielectric disk to the same size with a circular hole in the center of the dielectric, exposing a region of conductor. Results of the electrostatic plasma simulation are presented both with and without secondary electron emission from the dielectric; characteristic curves and voltage profiles are included.
Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.
Abstract not provided
Photon sources produced by laser beams with moderate laser intensities, up to 10{sup 14 }W/cm{sup 2}, are being developed for many industrial applications. The performance requirements for high volume manufacture devices necessitate extensive experimental research supported by theoretical plasma analysis and modeling predictions. We simulated laser produced plasma sources currently being developed for several applications such as extreme ultraviolet lithography using 13.5% ± 1% nm bandwidth, possibly beyond extreme ultraviolet lithography using 6.× nm wavelengths, and water-window microscopy utilizing 2.48 nm (La-α) and 2.88 nm (He-α) emission. We comprehensively modeled plasma evolution from solid/liquid tin, gadolinium, and nitrogen targets as three promising materials for the above described sources, respectively. Results of our analysis for plasma characteristics during the entire course of plasma evolution showed the dependence of source conversion efficiency (CE), i.e., laser energy to photons at the desired wavelength, on plasma electron density gradient. Our results showed that utilizing laser intensities which produce hotter plasma than the optimum emission temperatures allows increasing CE for all considered sources that, however, restricted by the reabsorption processes around the main emission region and this restriction is especially actual for the 6.× nm sources.
Transport modeling based on inclusion of turbulence spreading and residual stresses shows internal rotation reversals and polarity reversal of cold pulses, with a clear indication of nonlocal transport effects due to fast spreading in the turbulence intensity field. The effects of turbulence spreading and residual stress are calculated from the gradient of the turbulence intensity. In the model presented in this paper, the flux is carried by the turbulence intensity field, which in itself is subject to radial transport effects. The pulse polarity inversion and the rotation profile reversal positions are close to the radial location of the stable/unstable transition. Both effects have no direct explanation within the framework of classical transport modeling, where the fluxes are related directly to the linear growth rates, the turbulence intensity profile is not considered and the corresponding residual stress is absent. Our simulations are in qualitative agreement with measurements from ohmically heated plasmas. Rotation reversal at a finite radius is found in situations not displaying saturated confinement, which we identify as situations where the plasma is nearly everywhere unstable. As an additional and new effect, the model predicts a perturbation of the velocity profile following a cold pulse from the edge. This allows direct experimental confirmation of both the existence of residual stress caused by turbulence intensity profiles and fundamental ideas of transport modeling presented here.
We report temporally resolved simultaneous measurements of the turbulent Reynolds stresses in both the parallel and perpendicular directions and the corresponding particle fluxes in the fusion relevant cylindrical magnetized plasma device Controlled Shear Decorrelation eXperiment (CSDX). CSDX simulates the plasma conditions of multiple plasma instabilities that can arise in the scrape-off layer of fusion devices. In this study, we designed and used a 6-tip Langmuir probe in a novel yet simple design to simultaneously measure all the three dimensional components (radial, azimuthal, and axial) of fluctuations in velocity from the floating potentials and plasma densities with high temporal resolution. From these, we calculated the parallel and perpendicular Reynolds stress and the particle fluxes in addition to the density and potential spectra and the cross phase between different quantities. We can obtain radial profiles of all the aforementioned plasma quantities, which are extremely useful for studying plasma turbulence due to multiple instabilities. We have also cross-checked the time averaged velocity profiles from the probe with laser induced fluorescence measurements of the mean plasma velocity for some common plasma source parameters.
The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (NCSX). Zonal flow oscillations and long-time damping affect the zonal flow evolution, and the zonal flow residual goes to zero for small radial wavenumber. The oscillation frequency and damping rate depend on the bounce-averaged radial particle drift in accordance with theory. While each flux tube on a flux surface is unique, several different flux tubes in HSX or NCSX can reproduce the zonal flow damping from a flux-surface calculation given an adequate parallel extent. The flux-surface or flux-tube calculations can accurately reproduce the full-volume long-time residual for moderate kx, but the oscillation and damping time scales are longer in local representations, particularly for small kx approaching the system size.
Global neoclassical simulations of a model equilibrium of the W7-X stellarator find an ambipolar electric field with either an ion root or an electron root solution depending on the temperature ratio between electrons and ions. The ambipolar electric field is then used as an equilibrium field in the turbulence simulations of ion temperature gradient (ITG) instability. The shear of the ambipolar electric field has modest effects on the ITG linear instability, nonlinear saturation amplitude, and turbulent transport in the ion root case. However, in the electron root case, the ambipolar electric field significantly reduces the linear ITG growth rate, turbulence intensity, and radial correlation length, resulting in an ion heat conductivity comparable to the neoclassical transport level in the strong shear region.
The incorporation of lower hybrid (LH) wave spectrum broadening in the poloidal wavenumber (k θ ) space at the last close field surface (LCFS) is reported to lead to better agreement of the modeled LH wave current profile with the experimental results [Baek et al., Nucl. Fusion 61, 106034 (2021)]. To further understand its underlying mechanism and find the possible influence factors, effects of wave scattering caused by drift-wave type density fluctuation on the probability distribution of the LH wave polar refractive index (N θ ) at the LCFS are studied under density parameters in the scrape-off-layer. According to a scattering model [P. T. Bonoli and E. Ott, Phys Fluids 25(2), 359–375 (1982)], scattering probability and scattering angle distribution are two main factors that determine the degree of spectral broadening. Studies presented here show that the total scattering probability increases first and then decreases as the wave propagates toward a smaller normalized radius of poloidal magnetic flux (ρ). The degree of spectral broadening is found to depend on the density and density fluctuation together by changing the intensity and a proportion of the geometrical optics approximation term and the E×B drift term in the scattering model. Additionally, the fluctuation correlation length can significantly modify the probability distribution of N θ at the LCFS, which is found to significantly change the LH wave current profile.
Not Available
Abstract not provided.
Abstract not provided.
Explore the source record for details and available documents.
Explore the source record for details and available documents.