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 accelerator compressible turbulent boundary layers numerical calculation taking into account MHD effects, electron thermal nonequilibrium and finite rate ionization
Coronal active stars of seven classes between spectral types F and M, single and double, are compared in their quiescent radio and X-ray luminosities L(R) and L(X). We find, largely independent of stellar class, log L(X) is less than about log L(R) + 15.5. This general relation points to an intimate connection between the nonthermal, energetic electrons causing the radio emission and the bulk plasma of the corona responsible for thermal X-rays. The relation, observed over six orders of magnitude, suggests that the heating mechanism necessarily involves particle acceleration. We derive requirements for simple models based on optically thin gyrosynchrotron emission of mildly relativistic electrons and thermal X-rays from the bulk plasma. We discuss the possibility that a portion of the accelerated particles heats the ambient plasma by collisions. More likely, plasma heating and particle acceleration may occur in parallel and in the same process, but at a fixed ratio.
Spectrally resolved passive Balmer- α (D- α, H- α) measurements from the DIII-D 16 channel edge main-ion charge exchange recombination system confirm the presence of higher energy neutrals (“thermal” neutrals) in addition to the cold neutrals that recycle off the walls in the edge region of DIII-D plasmas. Charge exchange between thermal ions and edge neutrals transfers energy and momentum between the populations giving rise to thermal neutrals with energies approximating the ions in the pedestal region. Multiple charge exchange events in succession allow an electron to effectively take a random walk, transferring from ion to ion, providing a pathway of increasing energy and velocity, permitting a neutral to get deeper into the plasma before a final ionization event that contributes to the ion and electron particle fueling. Spectrally resolved measurements provide information about the density and velocity distribution of these neutrals, which has been historically valuable for validating Monte Carlo neutral models, which include the multi stage charge exchange dynamics. Here, in this study, a multi-channel set of such measurements is used to specifically isolate the details of the thermal neutrals that are responsible for fueling inside the pedestal top. Being able to separate the thermal from the cold emission overcomes several challenges associated with optical filter-based neutral density measurements. The neutral dynamics, deeper fueling by the thermal neutrals, and spectral measurement are modeled with the FIDASIM Monte Carlo collisional radiative code, which also produces synthetic spectra with a shape that is in close agreement with the measurements. By scaling the number of neutrals in the simulation to match the intensity of the thermal emission, we show it is possible to obtain local neutral densities and ionization source rates.
Abstract Nonlinear 3D MHD simulations of shattered-pellet injection (SPI) in JET show prototypical SPI-driven disruptions using the M3D-C1 and NIMROD extended-MHD codes. Initially, radiation-driven thermal quenches are accelerated by MHD activity as the pellet crosses rational surfaces, leading to a radiation spike, global stochasticization of the magnetic field, and a complete thermal quench. Eventually, current quenches, preceded by a current spike are seen as the Ohmic heating becomes equal to the radiative cooling. The results are qualitatively similar for both a single monolithic pellet, pencil-beam model, and a realistic shatter to represent the SPI plume. A scan in viscosity from 500 to 2000 m 2 s −1 for MHD simulations finds that reducing viscosity increases MHD activity and decreases thermal quench time slightly. A realistic cloud of fragments modeling shows that mixed-D–Ne pellet travels deeper into the plasma core before the thermal quench. At the slow pellet speeds, the pellet is found to be moving slowly enough inward that even the 5% neon in the mixed pellet is enough to effectively radiate the thermal energy available. Radiation toroidal peaking is predicted to be at levels consistent with experimental observations and reduced as the pellet travels deeper into the plasma. These simulations lay the ground work for more-sophisticated validative and predictive modeling of SPI in JET using both M3D-C1 and NIMROD.
The influences of ion trapping and fluctuations of electron temperature and plasma flow on cross-beam energy transfer (CBET) are examined using two- and three-dimensional particle-in-cell simulations in parameter regimes relevant to recent CBET experiments at the OMEGA laser facility. In mid-Z plasma irradiated by an intense pump beam and weaker probe beam, ion trapping, collisional de-trapping, and plasma flow induced by thermal effects are shown to affect the CBET gain. Ion trapping can enhance or detune the CBET resonance [Nguyen et al., Phys. Plasmas 28, 082705 (2021)]. Collisional de-trapping can affect the CBET gain at low seed beam intensity near the onset threshold for ion trapping. Thermal-effects-induced flow can also detune the CBET resonance at a level comparable to that from trapping at low seed beam intensity. As a consequence, the CBET gain is sensitive to collisions and dimensionality at low seed beam intensity where ion trapping is weak but is insensitive to collisions and dimensionality at high seed beam intensity where ion trapping is strong.
Plasma-filled loop structures are common in the solar corona. Because detailed modeling of the dynamical evolution of these structures is computationally costly, an efficient method for computing approximate but quick physics-based solutions is to rely on space-integrated 0D simulations. The enthalpy-based thermal evolution of loops (EBTEL) framework is a commonly used method to study the exchange of mass and energy between the corona and transition region. EBTEL solves for density, temperature, and pressure, averaged over the coronal part of the loop, velocity at coronal base, and the instantaneous differential emission measure distribution in the transition region. The current single-fluid version of the code, EBTEL2, assumes that at all stages the flows are subsonic. However, sometimes the solutions show the presence of supersonic flows during the impulsive phase of heat input. It is thus necessary to account for this effect. Here, we upgrade EBTEL2 to EBTEL3 by including the kinetic energy term in the Navier–Stokes equation. We compare the solutions from EBTEL3 with those obtained using EBTEL2, as well as the state-of-the-art field-aligned hydrodynamics code HYDRAD. We find that the match in pressure between EBTEL3 and HYDRAD is better than that between EBTEL2 and HYDRAD. Additionally, the velocities predicted by EBTEL3 are in close agreement with those obtained with HYDRAD when the flows are subsonic. However, EBTEL3 solutions deviate substantially from HYDRAD’s when the latter predicts supersonic flows. Using the mismatches in the solution, we propose a criterion to determine the conditions under which EBTEL can be used to study flows in the system.
The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.
The directions and energy distributions of neutral atoms and molecules ejected from Io are modeled analytically. The particles are emitted from the Io exosphere by collisions of incident plasma ions. Some of the Io atmosphere is attributed to particles sputtered off the surface by incoming plasma ions, while other departing particles originate at the top of volcanic plumes. Account is taken of the thermal motion of incident plasma and Io's gravitational field when calculating the distribution of the ejected particles. The distributions associated with each type of ejection force are projected to be different and highly dependent on the inclusion of realistic cross-sections. Attention is focused on Na-rich species and defining ejection processes which can be verified with ground-based observations. It is, however, concluded that no particles are ejected directly from the surface by Jean's escape.
A thermomechanical study of the intermittent rubbing of a rotor blade tip and a casing seal is carried out taking into account the existence of thermal contact resistance. The effect of the thermal properties of a plasma sprayed coating on the blade tip is studied. The influence of a variable heat generation and variable thermal contact resistance at the blade tip as it passes along the path of rubbing is also discussed.
Plasma loops in solar active regions have been observed in EUV and soft X-rays for decades. Their formation mechanism and properties, however, are still not fully understood. Predictions by early models, based on 1D hydrostatic equilibria with uniform plasma heating, are not consistent with high-resolution measurements. In this Letter, we demonstrate, via 3D simulations, that a class of heating models can lead to the dynamic formation of plasma loops provided the plasma is heated sufficiently to match SXT soft X-ray measurements. We show that individual flux tubes in a 3D magnetic structure tend to stand out against their neighbors. The loops have large aspect ratios and nearly uniform cross sections in the corona, similar to those observed by EIT and TRACE. The coronal EUV emission from these thermally unstable solutions is roughly consistent with EIT measurements. The solution oscillates in time through a large-amplitude, nonlinear cycle, leading to repeated brightening and fading of the loops.
This two-part paper reports the thermophysical properties of several cold and vacuum plasma sprayed monolithic Cu and Ni-based alloy coatings. Part I presents the electrical and thermal conductivity, thermal diffusivity, and total hemispherical emissivity data while Part II reports the specific heat capacity data for these coatings. Metallic copper alloys, stoichiometric NiAl and NiCrAlY coatings were fabricated by either the cold sprayed or the vacuum plasma spray deposition processes for thermal property measurements between 77 and 1223 K. The temperature dependencies of the thermal conductivities, thermal diffusivities, electrical conductivities and total hemispherical emissivities of these cold and vacuum sprayed monolithic coatings are reported in this paper. The electrical and thermal conductivity data correlate reasonably well for Cu-8%Cr-1%Al, Cu-23%Cr-5%Al and NiAl in accordance with the Wiedemann-Franz (WF) law although a better fit is obtained using the Smith-Palmer relationship. The Lorentz numbers determined from the WF law are close to the theoretical value.
The heating of He(+) ions observed by the GEOS satellite at L of about seven starting from an initial anisotropic proton distribution is studied. A simulation of electromagnetic ion cyclotron waves in a plasma consisting of hot anisotropic protons, a dominant thermal hydrogen plasma, and a minority species of helium ions is utilized. It is shown that the thermal ions are heated to 100 times their initial temperature due to the growth of large-amplitude ion cyclotron waves. He(+) ions are heated mainly in the perpendicular direction, and they are heated more than cold H(+) ions. The He(+) ions are heated in a two-step process. First, they are set into oscillations by the growing wave until some of the ions reach a parallel velocity on the order of the resonant velocity; then, strong heating occurs.
We study the structure of reverse shocks formed by the collision of supersonic, magnetized plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminum wire arrays, the shock is sharp and well-defined, consistent with magneto-hydrodynamic theory. In contrast, we do not observe a well-defined shock using tungsten wires, and instead we see a broad region dominated by density fluctuations on a wide range of spatial scales. We diagnose these two very different interactions using interferometry, Thomson scattering, shadowgraphy, and a newly developed imaging refractometer that is sensitive to small deflections of the probing laser corresponding to small-scale density perturbations. We conclude that the differences in shock structure are most likely due to radiative cooling instabilities, which create small-scale density perturbations elongated along magnetic field lines in the tungsten plasma. These instabilities grow more slowly and are smoothed by thermal conduction in the aluminum plasma.
Triboelectrically-charged objects may create threshold sparks, electrostatic discharge (ESD) events, to equilibrate charge between themselves and other relatively charged objects. ESD events exhibit many complex physical phenomena. They are a nexus of several fields of physics with disparate characteristic scales: plasma physics, chemical kinetics, hydrodynamics, circuit models, etc. These scales can span many orders of magnitude from the varied collisions thermalizing information within a plasma on the $\mathcal{O}(fs/ps)$ to the physical size of the plasma channel on the $\mathcal{O}(100µm)$, to the speed of a nonlinear hydrodynamic wave propagating at $\mathcal{O}(µm, ns)$. These threshold ESD events may occur in situations of programmatic importance, delivering energy and power profiles to a “victim load” generating deleterious consequences. To predict and mitigate these consequences we must answer questions about the spark’s energy budget: how much energy goes into producing the spark channel; how much gets radiated away; how much energy is advected away into the hydrodynamics; and how much energy is delivered to a victim load. An ESD simulation toolset has been created and evolved in order to answer these questions. An appropriate, physics guided implementation for simulation can be done by gaining insight into its constituent physics and leveraging that intuition to choose a suitable numerical operator. We examine in detail the chemical kinetics, circuit discharge, and hydrodynamics to deter mine dominant regimes, values, timescales, and interactions to uncover the underlying physical dynamics. We also examine and propose model reduction schemes for high-dimensional chemical kinetics. We use past and current work with experimentally validated and theoretically-verified hydrodynamics to calculate applicability limits of the non-ionizing strong shock limit. We quantify the energy budget from a hydrodynamic perspective and demonstrate that a significant fraction of the stored energy is “earmarked” for hydrodynamic advection as an energy terminus. Lastly, we combine the constituent physics of an ESD event (chemical kinetics, circuit model, and hydrodynamics) into a cohesive, actionable toolset and obtain promising results from an isothermal test case. We then propose a viable, modular evolution of the ESD toolset based upon the performed examination of the physics uncovering dominant physical scales and the stiffness of the compositional differential system.
Measurements of ions from three different instruments on the IMP-7 and 8 spacecraft are combined to yield the differential energy spectra of ions over the entire energy range of 100 eV to 4 MeV in the earth's distant (30 to 40 earth radii) plasma sheet. These spectra, obtained during times of relatively small bulk flow velocities, span the intensity range from 10 to the -5th to 10 to the 5th per sq cm-sec-sr-keV, varying smoothly over the entire energy range both when the plasma is cold (about 1 keV) and hot (about 9 keV). Overall, the shape of the spectrum resembles a Maxwellian but with a high energy (not less than 50 keV) tail described well by a power law. The high energy tail is displaced in a parallel fashion to higher or lower intensities when the plasma is hot or cold, respectively. It is found that the energetic particle populations in the plasma sheet appear to be directly related to the mean thermal energies of the corresponding plasmas.
Thermal barrier and environmental barrier coatings (TBCs and EBCs) will play a crucial role in future advanced gas turbine engines because of their ability to significantly extend the temperature capability of the ceramic matrix composite (CMC) engine components in harsh combustion environments. In order to develop high performance, robust coating systems for effective thermal and environmental protection of the engine components, appropriate test approaches for evaluating the critical coating properties must be established. In this paper, a laser high-heat-flux, thermal gradient approach for testing the coatings will be described. Thermal cyclic behavior of plasma-sprayed coating systems, consisting of ZrO2-8wt%Y2O3 thermal barrier and NASA Enabling Propulsion Materials (EPM) Program developed mullite+BSAS/Si type environmental barrier coatings on SiC/SiC ceramic matrix composites, was investigated under thermal gradients using the laser heat-flux rig in conjunction with the furnace thermal cyclic tests in water-vapor environments. The coating sintering and interface damage were assessed by monitoring the real-time thermal conductivity changes during the laser heat-flux tests and by examining the microstructural changes after the tests. The coating failure mechanisms are discussed based on the cyclic test results and are correlated to the sintering, creep, and thermal stress behavior under simulated engine temperature and heat flux conditions.
An important question in solar physics is whether solar microflares, the smallest currently observable flare events in X-rays, possess the same energetic properties as large flares. Recent surveys have suggested that microflares may be less efficient particle accelerators than large flares, and hence contribute less non-thermal energy, which may have implications for coronal heating mechanisms. We therefore explore the energetic properties of microflares by combining EUV and X-ray measurements. We present forward-fitting differential emission measure (DEM) analysis of 10 microflares. The fitting is constrained by combining, for the first time, high-temperature Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) observations and flux data from the Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA). Two fitting models are tested for the DEM; a Gaussian distribution and a uniform DEM profile. A Gaussian fit proved unable to explain the observations for any of the studied microflares. However, 8 of 10 events studied were reasonably fit by a uniform DEM profile. Hence microflare plasma can be considered to be significantly multi-thermal, and may not be significantly peaked or contain resolvable fine structure, within the uncertainties of the observational instruments. The thermal and non-thermal energy is estimated for each microflare, comparing the energy budget with an isothermal plasma assumption. From the multi-thermal fits the minimum non-thermal energy content was found to average approximately 30% of the estimated thermal energy. By comparison, under an isothermal model the non-thermal and thermal energy estimates were generally comparable. Hence, multi-thermal plasma is an important consideration for solar microflares that substantially alters their thermal and non-thermal energy content.
We developed a technique to measure the thermal release profiles of a suite of labile elements (Zn, As, Se, Cd, In, Sn, Sb, Te, Pt, Hg, Au, Tl, Pb, Bi). Conclusions are reached about the behavior of each element during parent-body alteration. Additional information is contained in the original extended abstract.