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.
The study will develop mathematical and numerical models of the turbulent ignition process that comprises plasma species interaction with turbulence and principles of ignition kernel growth or extinction based on a flamelet modeling approach.
It is shown here that of the three plasma models for the generation of broadband noise in the magnetotail proposed by Dusenbery and Lyons (1985), two result in the same types of instabilities that are also excited in the Grabbe and Eastman (GE, 1984) model, while the third model introduces other modes not present in the GE model. For the plasma parameters given in the GE model, the ion/ion acoustic instability is a nonresonant one, even though the phase velocity of the excited waves falls within the distribution function of the core ions. The nonlinear saturation mechanism of the ion/ion acoustic instability is the trapping of both the core and the beam ions. In a reply, Dusenbery addresses several ongoing controversies which have resulted from studies of wave particle interactions in the plasma sheet boundary layer.
Numerical modeling, combining the toroidal ideal MHD code GPEC and the nonlinear two-fluid MHD code TM1, was used for comprehensive studies of the plasma response to resonant magnetic perturbations (RMPs) with toroidal mode number n = 1–5 for controlling edge-localized modes (ELMs) in ITER for the standard operation scenario (15 MA Q = 10). Several issues related to RMP ELM control are investigated, including the optimization of the RMP coils configuration, the evaluation of the magnitude of density pump-out and the q95 windows of ELM suppression. GPEC calculates the magnetic response, which consistently includes the very important edge kink/peeling response to static magnetic perturbations. Furthermore, GPEC two-dimensional scans of the relative coil current phasing among the three rows of internal coils, at fixed coil current amplitude, reveal the optimal phasing for the RMP coil configuration with n = 1–5, respectively. The poloidal half wavelength of resonant mode at the edge of plasma calculated by GPEC indicates that the midplane row coils have the best resonant coupling with the plasma for n = 2, while the upper and lower row coils have the best resonant coupling with the plasma for n = 3. Based on the plasma kinetic equilibrium and the GPEC calculations of the magnetic response, TM1 was used to simulate the conditions for RMP field penetration in the ITER pedestal. TM1 shows magnetic island formation at the foot of ITER pedestal with RMP coil current threshold ranging from 4 kAt to 8 kAt with n = 2 to 4. These magnetic islands at the pedestal-foot lead to density pump-out, the magnitude of which scales as ${I}_{\text{RMP}}^{0.5}$ and ranges from 5% to 20% at the pedestal-top when scanning the coil current from 4 to 60 kAt. The density pump-out is found to be weaker for higher n RMP. The nonlinear TM1 simulations also show field penetration at the pedestal-top, where the threshold of RMP coil current depends on the q 95 . The alignment of the magnetic island and the location of the pedestal-top decreases the height and width of the pedestal to suppress ELMs. Simulations by two-dimensional scans of RMP coil current and q 95 reveal the accessible q 95 windows of ELM suppression for both n = 3 and 4 RMPs. The predicted q 95 windows of ELM suppression are very similar to the ones in currently operating tokamaks and the required RMP coil current for ELM suppression is less than 40–50 kAt, which is well within the designed capability for ITER. In addition, the simulations indicate that wide q 95 windows of ELM suppression may be accessible in ITER by operating with dominant n = 4 (or n = 5) RMPs.
The angular distributions of ions in the wake of a body moving through a space plasma computed from three different models are compared in order to investigate wake current depletion ratios under conditions simulating the topside ionosphere and plasmasphere. Results demonstrate the importance of taking into account the thermal flux at low Mach numbers and the angular acceptance of ion detectors in making theory-experiment comparisons. For all models considered, gradients in the angular variations of the fluxes are shown to be steeper near the wake-ambient interface than closer to the maximum rarefaction region.
Recent simulation and experimental results suggest that the magnetic island and flow on a resonant surface often do not satisfy the “no-slip” condition in the steady state. A new theoretical model on a nonlinear plasma response to external magnetic perturbation in the absence of a no-slip condition is proposed. The model is composed of the equations for the evolution of both the width and phase of the magnetic island due to forced reconnection driven by the external magnetic perturbation, and the force-balance equation for the plasma flow. When the island width is much less than the resistive layer width, the island growth is governed by the linear Hahm–Kulsrud–Taylor solution in the presence of time-dependent plasma flow. In the other regime when the island width is much larger than the resistive layer width, the evolution of both the island width and phase can be described using the Rutherford theory. The island solution is used to construct the quasilinear electromagnetic force, which, together with the viscous one, contributes to the nonlinear variation in plasma flow. The no-slip condition assumed in the conventional error field theory is not imposed here, where the island oscillation frequency depends on but does not necessarily equal to the plasma flow frequency at the rational surface.
A thermal equilibrium plasma model is used to process data from an impact ionization time-of-flight mass spectrometer in order to convert the raw ion data to relative abundances of the elemental constituents of cosmic dust particles.
AE Aquarii is a cataclysmic variable with the fastest known rotating magnetized white dwarf (P(sub spin) = 33.08 s). Compared to many intermediate polars, AE Aquarii shows a soft X-ray spectrum with a very low luminosity (LX (is) approximately 10(exp 31) erg per second). We have analyzed overlapping observations of this system with the NuSTAR and the Swift X-ray observatories in 2012 September. We find the 0.5-30 keV spectra to be well fitted by either an optically thin thermal plasma model with three temperatures of 0.75(+0.18 / -0.45), 2.29(+0.96 / -0.82), and 9.33 (+6.07 / -2.18) keV, or an optically thin thermal plasma model with two temperatures of 1.00 (+0.34 / -0.23) and 4.64 (+1.58 / -0.84) keV plus a power-law component with photon index of 2.50 (+0.17 / -0.23). The pulse profile in the 3-20 keV band is broad and approximately sinusoidal, with a pulsed fraction of 16.6% +/- 2.3%. We do not find any evidence for a previously reported sharp feature in the pulse profile.
The RF-Transpond code couples a fluid plasma transport solver with a frequency domain cold plasma RF wave solver in a 1D domain parallel to a strong background magnetic field. A ponderomotive force term proportional to parallel gradients in the electric field strength is included in the transport model in order to describe ponderomotive effects in the scrape-off layer (SOL) of fusion plasmas. The transport and wave codes are verified independently and a coupled case corresponding to experimental parameters from the LArge Plasma Device (LAPD) is presented. The density perturbation ratio R n , calculated to describe ponderomotive force driven modifications, is up to 20% for the simulation inputs used.
In this fourth paper in a series, we present a model of the remarkable temporal and azimuthal variability of the Io plasma torus observed during the Cassini encounter with Jupiter. Over a period of three months, the Cassini Ultraviolet Imaging Spectrograph (UVIS) observed a dramatic variaton in the average torus composition. Superimposed on this long-term variation, is a 10.07-hour periodicity caused by azimuthal variation in plasma composition subcorotating relative to System III longitude. Quite surprisingly, the amplitude of the azimuthal variation appears to be modulated at the beat frequency between the System III period and the observed 10.07-hour period. Previously, we have successfully modeled the months-long compositional change by supposing a factor of three increase in the amount of material supplied to Io's extended neutral clouds. Here, we extend our torus chemistry model to include an azimuthal dimension. We postulate the existence of two azimuthal variations in the number of superthermal electrons in the torus: a primary variation that subcorotates with a period of 10.07 hours and a secondary variation that remains fixed in System III longitude. Using these two hot electron variations, our model can reproduce the observed temporal and azimuthal variations observed by Cassini UVIS.
The influence of plasma beta effects on the edge plasma transport in the Wendelstein 7-X standard configuration is studied systematically by using EMC3-EIRENE combined with a 3D equilibrium code named HINT. The magnetic topology changes induced by plasma beta effects are significantly reflected in plasma transport behaviors and heat flux patterns on divertor targets. After validating the modeling strategy by comparisons with experimental data, the extended simulations for high performance plasmas show that the threshold separatrix density for accessing the power detachment is reduced in higher beta plasmas. Compared with the vacuum field case, the impurity radiation distributions with finite beta effects are modified in the magnetic island region. The divertor heat flux is distributed more evenly along the toroidal direction on the strike line at the vertical target. The strike line on the horizontal target moves towards the pumping gap with an increase in the plasma beta. In addition, the different pressure profiles with the same central beta also result in a modified heat flux pattern on the divertor targets.
We’ve achieved our technical goals: (1) Perform large-scale kinetic simulations of SRS and CBET a range of plasma and laser conditions; (2) Develop a new $\delta f$-Gaussian-mixture algorithm to represent non-Maxwellian distribution functions from particle trapping and time-dependent plasma response; (3) Develop physics-based preliminary nonlinear LPI models for coupling to LRT/ Mazinisin and AMP; and (4) Begin implementation of LPI sub-grid models into Mazinisin.
Reduction of both the plasma density and toroidal flow speed, due to application of the predominantly n = 1 (n is the toroidal mode number) resonant magnetic perturbation (RMP) for controlling the edge localized mode in the HL-2A tokamak, is numerically investigated utilizing the quasi-linear initial-value code MARS-Q (Liu et al 2013 Phys. Plasmas 20 042503). Simulation results reveal that the neoclassical toroidal viscosity (NTV) due to three dimensional fields plays the key role in modifying the plasma momentum and particle transport in the HL-2A discharge. By comparing the modeling results with the measured density pump-out in the experiment, the electron NTV particle flux model, in combination with the free-boundary condition for the axisymmetric change of the density at the plasma edge, is found to yield the best agreement in terms of both the pump-out level and the overall time scale. Further sensitivity studies show that the simulated density pump-out level is reasonably robust against variations in the model assumptions, including the particle diffusion model and the non-ambipolar versus ambipolar NTV particle flux. The latter however affects the time scale for reaching the steady state solution. Finally, it is found that the plasma edge-peeling response, the NTV torque, as well as the plasma momentum and particle transport, all are sensitive to the toroidal phase difference between the upper and lower rows of the RMP coil currents in HL-2A, with the 30 ∘ coil phasing producing the minimal side effects on the plasma.
We review the status of theoretical work on magnetospheric convection in the lower auroral zone and at midlatitudes, and compare with various observed features, such as ionospheric electric fields and plasma flows, the form of the plasmasphere, and the distribution of plasma-sheet particles. We present preliminary results from a new series of computerized convection models, which follow the time evolution of the inner magnetosphere (L less than around 10), self-consistently including ionospheric currents and Birkeland currents, as well as the currents generated in a model plasma sheet with a realistic energy spectrum. We find that the model plasma sheet's inner edge quickly becomes rather sharp. Computed electric field distributions resemble those obtained earlier for a simple single-energy plasma sheet.
The Dynamic Global Core Plasma Model (DGCPM) is a parameterized model for core or thermal plasma in the magnetosphere. The model accounts for dayside ionospheric outflow and nightside inflow. It accounts for the global pattern of convection and corotation. The model is capable of being coupled to ring current and superthermal electron models for the purpose of providing thermal plasma spatial distributions and for the purpose of accepting the dynamic influences of these plasma populations back upon the thermal plasma. The DGCPM is designed to operate alone or to operate as part of a larger integrated package. The convection electric field and magnetic field used within the DGCPM can be shared with models of other plasma populations, in addition to the exchange of parameters important to the collective modeling of whole plasma systems in the inner magnetosphere. This talk will present the features of the DGCPM model code and the various forms of information that can be exchanged with other cooperating codes.
A new model for electron temperature gradient (ETG) modes is developed as a component of the multi-mode anomalous transport module [Rafiq et al., Phys Plasmas 20, 032506 (2013)] to predict a time-dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, β and gradient of β (β'), collisionality, safety factor, and toroidicity, where β is the ratio of the plasma pressure to the magnetic pressure. Furthermore, the electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.
A new model for electron temperature gradient (ETG) modes is developed as a component of the Multi-Mode anomalous transport module [T. Rafiq \textit{et al.,} Phys Plasmas \textbf{20}, 032506 (2013)] to predict a time dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, $\beta$ and gradient of $\beta$ $(\betap)$, collisionality, safety factor, and toroidicity, where $\beta$ is the ratio of plasma pressure to the magnetic pressure. The electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.
Research concerning the transport and distribution of ionospheric plasma in the magnetosphere are reviewed, stressing the dichotomy in explanations given for the low plasma densities outside the plasmasphere. The convection/hot solar plasma model and the convection/loss model are considered. Observations of global ionospheric outflows are compared with theoretical studies. It is suggested that there is a need for a hybrid model of magnetospheric plasma in which terrestrial plasma is both lost into the solar wind and energized and trapped within the magnetosphere, inflating the geomagnetic field and excluding cold plasma from conjugate regions.
The development of lithium magnetoplasmadynamic (MPD) thrusters at JPL is discussed. The following topics are presented in vugraph form: mercury vapor mass flow control; porous tungsten vaporizer and housing; the lithium vaporizer experiment; a dry box for handling solid lithium; MPD thruster electrode modeling; engine lifetime definitions; cathode failure modeling; cathode erosion modeling; cathode thermal modeling; near cathode plasma model regions; cathode work function modeling; anode work function modeling; and radiation-cooled anodes.