SEARCH · Search NASA
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.
In recent years, the Global Positioning System (GPS) has been exploited via radio occultation techniques to obtain profiles of refractivity, temperature, pressure and water vapor in the neutral atmosphere and electron density in the ionosphere.
In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.
In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.
Direct sublimation of a comet nucleus surface is usually considered to be the main source of gas in the coma of a comet. However, evidence from a number of comets including the recent spectacular images of Comet 103P/Hartley 2 by the EPOXI mission indicates that the nucleus alone may not be responsible for all, or possibly at times even most, of the total amount of gas seen in the coma. Indeed, the sublimation of icy grains, which have been injected into the coma, appears to constitute an important source. We use the fully-kinetic Direct Simulation Monte Carlo model of Tenishev et al. to reproduce the measurements of column density and rotational temperature of water in Comet 73P-B/Schwassmann-Wachmann 3 obtained with a very high spatial resolution of approx. 30 km using IRCS/Subaru in May 2006. For gas released solely from the cometary nucleus at a heliocentric distance of 1 AU, modeled rotational temperatures start at 110 K close to the surface and decrease to only several tens of degrees by 10-20 nucleus radii. However, the measured decay of both rotational temperature and column density with distance from the nucleus is much slower than predicted by this simple model. The addition of a substantial (distributed) source of gas from icy grains in the model slows the decay in rotational temperature and provides a more gradual drop in column density profiles. Together with a contribution of rotational heating of water molecules by electrons, the combined effects allow a much better match to the IRCS/Subaru observations. From the spatial distributions of water abundance and temperature measured in 73P/SW3-B, we have identified and quantified multiple mechanisms of release. The application of this tool to other comets may permit such studies over a range of heliocentric and geocentric distances.
New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.
Vibrational and magnetic properties of single-crystal uranium-thorium dioxide (U 𝑥 Th 1−𝑥 O 2 ) with a full range of 0 < 𝑥 < 1 are investigated. Thorium dioxide is a diamagnet whose thermal properties are governed by lattice vibration. The addition of paramagnetic uranium ions leads to the emergence of magnetic effects that alter the thermophysical properties noticeably even at room temperature. The interaction of phonons with magnetic moments of uranium 5𝑓 electrons mediated by magnetoelastic coupling results in an anomalous low-temperature thermal conductivity profile. Analysis of the magnetic susceptibility measurements indicates a departure from the Curie-Weiss relationship characteristic of noninteracting paramagnetic ions, previously associated with the dynamic Jahn-Teller (DJT) effect characterized by coupling between spin and the oxygen sublattice. The T 2g Raman peak position follows a nonlinear trend as a function of uranium concentration and hints that these Raman active optical modes play a role in either DJT or in mediating quadrupole-quadrupole interactions. A first-principle-based thermal transport model is implemented to explain the low-temperature transport measurements, where the anomalous reduction is attributed to phonon-spin resonant scattering. The interplay between spins and phonons is also captured using high-resolution inelastic x-ray scattering (IXS) measurements of phonon linewidths. Our results provide insights into the phonon interactions with the magnetic excitations governing DJT effect and impacting the low-temperature thermal transport processes in this material system. Furthermore, these findings have implications for understanding low-temperature thermal transport and magnetic properties in advanced materials for information processing and energy applications.
Using a recently installed impurity powder dropper (IPD), boron powder (< 150 μm) was injected into lower single null (LSN) L-mode discharges in WEST. IPDs possibly enable real-time wall conditioning of the plasma-facing components and may help to facilitate H-mode access in the full-tungsten environment of WEST. The discharges in this experiment featured Ip = 0.5 MA, BT = 3.7 T, q95 = 4.3, tpulse = 12–30 s, ne,0 ~ 4×1019 m-2, and PLHCD ~ 4.5 MW. Estimates of the deuterium and impurity particle fluxes, derived from a combination of visible spectroscopy measurements and their corresponding S/XB coefficients, showed decreases of ~ 50% in O+, N+, and C+ populations during powder injection and a moderate reduction of these low-Z impurities (~ 50%) and W (~ 10%) in the discharges that followed powder injection. Along with the improved wall conditions, WEST discharges with B powder injection observed improved confinement, as the stored energy WMHD, neutron rate, and electron temperature Te increased significantly (10–25% for WMHD and 60–200% for the neutron rate) at constant input power. These increases in confinement scale up with the powder drop rate and are likely due to the suppression of ion temperature gradient (ITG) turbulence from changes in Zeff and/or modifications to the electron density profile.
This paper describes new modeling capabilities for predicting H-mode pedestal profiles in spherical tokamaks. Temperature profiles for NSTX discharges 132543 and 132588 are modeled by coupling the \textsc{astra} transport solver with neoclassical transport and gyrokinetic-based reduced models for electron temperature gradient (ETG) and kinetic ballooning mode (KBM) instabilities. A quasi-linear surrogate model for ion-scale transport is developed using linear \textsc{gene} simulations, requiring only a single free parameter calibrated to one discharge. Time-evolving the temperatures with fixed density yields good agreement with experiments for both discharges. Systematic analysis of the transport mechanisms reveals that neoclassical transport is huge across the entire pedestal region for the ion channel. ETG turbulence is large in the plasma edge and low density gradient region, contributing substantially to the electron channel. However, KBM/MHD-like modes also drive significant transport in both the ion and electron thermal channels, making them essential for accurate pedestal modeling. Further refinements, including explicit E×B shear suppression and scaled ETG transport, produce quantitative but not qualitative improvements. This work lays the foundation for predictive modeling of future devices. This paper is on arxiv and has been submitted to Nuclear Fusion.
Measured and theoretical stark-broadened line profiles and asymptotic wing approximation of hydrogen Balmer lines at specific electron density and temperature
The importance of the spectral range from about 80 to 800 A for determining physical conditions in different regions of the solar atmosphere is discussed. Examples are given of line ratios that may be used to determine electron densities in quiet-sun regions, active regions, and flares. The possibility of determining electron temperatures from line ratios in the EUV is considered. It is shown that profiles, as well as intensities, of spectral lines must be obtained for a proper interpretation of the spectra. Approximate parameters are provided for a solar grazing-incidence spectrograph suitable for the study of the 80-800-A wavelength region.
Self-consistent transonic multifluid solutions for the dynamical and thermal structure of the two-phase dusty gas atmosphere of a comet containing a dirty H2O-clathrate ice nucleus have been obtained. This was done by solving a simultaneous set of differential equations representing conservation of number density, momentum, and energy, as well as the transfer of solar radiation responsible for the major photolytic processes and nuclear heating. The results presented include velocity and Mach number profiles of heavy species and the velocity of dust within the subsonic and supersonic regions, temperature profiles of the dust and gas in the supersonic region, and radial profiles of electron and ion densities in the collision-dominated region.
One of the most convenient and effective methods of stressing a part or structure for shearographic evaluation is thermal excitation. This technique involves heating the part, often convectively with a heat gun, and then monitoring with a shearography device the deformation during cooling. For a composite specimen, unbonds, delaminations, inclusions, or matrix cracking will deform during cooling differently than other more structurally sound regions and thus will appear as anomalies in the deformation field. However, one of the difficulties that cause this inspection to be dependent on the operator experience is the conventional heating process. Fanning the part with a heat gun by hand introduces a wide range of variability from person to person and from one inspection to the next. The goal of this research effort was to conduct research in the methods of thermal excitation for shearography inspection. A computerized heating system was developed for inspection of 0.61 m (24 in.) square panels. The Thermal Excitation System for Shearography (TESS) provides radiant heating with continuous digital measurement of the surface temperature profile to ensure repeatability. The TESS device functions as an accessory to any electronic shearography device.
Increasingly long-duration space missions require a better understanding of the failure mechanisms of hollow cathodes to enable prediction of lifetime. Energetic ions have been measured in the plume of hollow cathodes but the energies are higher than can be explained classically. One of the mechanisms that appears to be contributing to the ion heating is ion acoustic turbulence (IAT), which has been investigated with plasma probes recently at JPL. In this work the ion velocity distribution has been measured using laser-induced fluorescence (LIF) at several high-current operating conditions on the centerline of a 100 A class hollow cathode at JPL. A retarding potential analyzer was (RPA) was used to measure the ion energy profile in the far plume and a Langmuir probe was used to measure electron temperature and plasma potential fluctuations. RPA energies are shown to correlate with the values measured by LIF. The ion temperature is shown to increase with distance downstream of the cathode. A parametric study of the ion heating as a function of discharge current and cathode flow rate was performed. It was found that for increasing current the ion temperature evolution did not vary even though the wave IAT wave energy increased. The effect of flow rate was more pronounced, with higher values reducing both ion temperature and wave energy linearly. These results indicate that IAT is at least partially responsible for energetic ion production in the hollow cathode plume.
To validate nonlinear gyrokinetic simulations of electron temperature turbulence, the experimental correlation electron cyclotron emission (CECE) measurements are to be compared using a synthetic CECE diagnostic, which generates modeled CECE measurement quantities by implementing realistic measurement parameters (e.g., spatial and wavenumber resolutions, radial location, etc.) to nonlinear gyrokinetic simulations. In this work, we calculate the radial and vertical spatial and wavenumber transfer functions, which are defined by the electron cyclotron emission emissivity radial profile and vertical probing antenna pattern, respectively. These transfer functions are applied to nonlinear gyrokinetic simulations of electron temperature turbulence using the continuum gyrokinetic code. A simultaneous comparison of the experimental electron temperature turbulence power spectrum and root-mean-square (RMS) level, as well as the radial correlation length with the new synthetic CECE diagnostic at a core location ρ ∼ 0.75 in an L-mode DIII-D tokamak plasma, is presented. Here, the preliminary result shows that the synthetic CECE output underestimates the RMS level by ∼42% and overestimates the radial correlation length by ∼40%.
The goal of the proposed research was to increase the understanding of coronal plasma phenomena by making use of different observational approaches and combine the observations with the necessary theoretical considerations. During the funding period we studied the formation of spectral lines in the corona and transition region. We compared the spectral line ratios that would arise from the type of temperature profile commonly used to explain in situ ion fractions with the actual observed line ratios. We also carried out a theoretical study to investigate how large the electron temperatures can be in the near sun region. We carried out more detailed studies to show that differential flow speeds between ions of the same element can not bridge the above gap between low coronal electron temperatures and high in situ ion fractions. To investigate the drift between core and halo in the electron distribution function, which is observed in situ in the solar wind, we developed the first solar wind model with two electron populations. It was found that the anomalous frictional forces acting on the halo electrons are the dominant factor inhibiting the core-halo drift in the East solar wind. We used kinetic modeling to investigate the Alfvenic turbulence in the extended corona. In these studies we also included the effects of proton heating.
The flow properties of the solar wind from 1 R(sub s) to 1 AU were obtained using a two fluid model constrained by density and scale height temperatures derived from white light observations, as well as knowledge of the electron temperature in coronal holes. The observations were obtained with the white light coronographs on SPARTAN 201-1 and at Mauna Loa (Hawaii), in a north polar coronal hole from 1.16 to 5.5 R(sub s) on 11 Apr. 1993. By specifying the density, temperature, Alfven wave velocity amplitude and heating function at the coronal base, it was found that the model parameters fit well the constraints of the empirical density profiles and temperatures. The optimal range of the input parameters was found to yield a higher proton temperature than electron temperature in the inner corona. The results indicate that no preferential heating of the protons at larger distances is needed to produce higher proton than electron temperatures at 1 AU, as observed in the high speed solar wind.
We present an application of the Portable Diagnostic Package (PDP) on the Wisconsin High-temperature-superconducting Axisymmetric Mirror (WHAM), which integrates an optical emission spectroscopy (OES) system and an active Thomson scattering (TS) system. The OES system facilitates a comprehensive impurity line survey and enables flow measurements through the Doppler effect observed on impurity lines. Plasma rotation profiles were successfully derived from doubly charged carbon lines. The TS system enabled the first measurements of the electron temperature in commissioning plasmas on WHAM. Notably, the PDP was installed, commissioned, and used to obtain OES and TS data within ∼6 months from project start, enabled by its designed portability and standardized interfaces. These results demonstrate the PDP’s potential to accelerate diagnostic readiness and advance experimental plasma studies.
Photoionization rates have been calculated from recent solar EUV flux measurements and are compared with the electron-ion recombination rates estimated from electron density profiles and current reaction rates. These photoionization rates also provided the first step in calculating the photoelectron heating rates of the ambient electrons, which are compared with the cooling rates deduced from ion and electron temperature measurements at Arecibo. The effects of changes in the assumed neutral models, the rate coefficients, and the temperature measurements are discussed. Comparisons of the profiles of energy input and loss with similar profiles of ion production and recombination suggest that the solar EUV fluxes have been underestimated. These comparisons also imply lower molecular neutral densities in winter than in summer.