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.
A solar oven heated by concentrated solar radiation as the heat source of a liquid metal magnetohydrodynamic (LMMHD) power generation system is proposed. The design allows the production of electric power in space, as well as on Earth, at high rates of efficiency. Two types of the solar oven suitable for the system are discussed.
Laser and optical emission-based measurements have been developed and implemented for use on low-power hydrogen arcjet thrusters and xenon-propelled electric thrusters. In the case of low power hydrogen arcjets, these laser induce fluorescence measurements constitute the first complete set of data that characterize the velocity and temperature field of such a device. The research performed under the auspices of this NASA grant includes laser-based measurements of atomic hydrogen velocity and translational temperature, ultraviolet absorption measurements of ground state atomic hydrogen, Raman scattering measurements of the electronic ground state of molecular hydrogen, and optical emission based measurements of electronically excited atomic hydrogen, electron number density, and electron temperature. In addition, we have developed a collisional-radiative model of atomic hydrogen for use in conjunction with magnetohydrodynamic models to predict the plasma radiative spectrum, and near-electrode plasma models to better understand current transfer from the electrodes to the plasma. In the final year of the grant, a new program aimed at developing diagnostics for xenon plasma thrusters was initiated, and results on the use of diode lasers for interrogating Hall accelerator plasmas has been presented at recent conferences.
Forward modeling is used to interpret inversion patterns of the pedestal–Scrape of Layer (SOL) Electron Cyclotron Emission (ECE) in DIII-D H-mode experiments. The modeling not only significantly improves the ECE data interpretation quality but also leads to the potential measurements of (1) the magnetic field strength |B| at the separatrix, (2) the pedestal |B| evolution during an inter-Edge Localized Mode (ELM) period, and (3) the pedestal Magnetohydrodynamics (MHD) radial structure. The ECE shine-through effect leads to three types of pedestal–SOL radiation inversions that are discussed in this paper. The first type of inversion is the non-monotonic T e,rad profile with respect to the major radius. Using the ECE frequency at the minimum T e,rad , the inversion can be applied to measure the magnetic field |B| at the separatrix and calibrate the mapping of the ECE channels with respect to the separatrix. The second type of inversion refers to the opposite phase between the radiation fluctuations δT e,rad at the pedestal and SOL. This δT e,rad phase inversion is sensitive to density and temperature fluctuations at the pedestal foot and, thus, can be used to qualitatively measure the MHD radial structure. The third type of inversion appears when the pedestal and SOL T e,rad evolve in an opposite trend, which can be used to infer the pedestal |B| field change during an inter-ELM period. Finally, the bandwidth effect on measuring δT e,rad due to pedestal MHD is also investigated in the radiation modeling.
Here we use a 3D radiative non-ideal magnetohydrodynamic simulation to investigate the formation and evolution of a young protostellar disc from a magnetized pre-stellar core. The simulation covers the first ${\sim }10\, {\rm kyr}$ after protostar formation and shows a massive, weakly magnetized disc with radius that initially grows and then saturates at ${\sim }30\, {\rm au}$. The disc is gravitationally unstable with prominent large-amplitude spiral arms. We use our simulation results and a series of physical arguments to construct a predictive and quantitative physical picture of Class 0/I protostellar disc evolution from several aspects, including (i) the angular-momentum redistribution in the disc, self-regulated by gravitational instability to make most of the disc marginally unstable; (ii) the thermal profile of the disc, well-approximated by a balance between radiative cooling and accretion heating; and (iii) the magnetic-field strength and magnetic-braking rate inside the disc, regulated by non-ideal magnetic diffusion. Using these physical insights, we build a simple 1D semi-analytic model of disc evolution. We show that this 1D model, when coupled to a computationally inexpensive simulation for the evolution of the surrounding pseudo-disc, can be used reliably to predict disc evolution in the Class 0/I phase. The predicted long-term evolution of disc size, which saturates at ${\sim }30\, {\rm au}$ and eventually shrinks, is consistent with a recent observational survey of Class 0/I discs. Such hierarchical modelling of disc evolution circumvents the computational difficulty of tracing disc evolution through Class 0/I phase with direct, numerically converged simulations.
Both near-field (<1 m away toroidally from the pellet) and far-field (>1 m away toroidally from the pellet) poloidal and toroidal impurity (carbon ion) flows are measured using visible imaging and fast bolometry during a polypropylene pellet-induced rapid plasma shutdown in the DIII-D tokamak. In the near field, the pellet appears to increase poloidal flow in the ion diamagnetic direction, possibly due to the strong radial temperature gradient caused by the pellet ablation. In the far field, the poloidal impurity flow typically appears slower and in the opposite direction. Toroidal impurity flow appears to be strongly influenced by the plasma initial toroidal rotation, especially in the far field. These results demonstrate that rapid shutdown impurity flows are not necessarily global in structure but can be quite different close to and far from the injected pellet.
We present the first end-to-end calculation connecting the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf to observable kilonova signatures, combining two-dimensional (2D) general-relativistic neutrino-magnetohydrodynamic simulations, followed by radiation hydrodynamics with in-situ nuclear network and 2D Monte Carlo radiative transfer with spatially resolved heating rates. Unlike all previous unmagnetized AIC models – which predicted proton-rich, $^{56}$Ni-dominated ejecta – strong magnetic fields eject ${\approx }\, 0.2\, \mathrm{ M}_\odot$ of neutron-rich material ($\langle Y_e \rangle \sim 0.24$) on dynamical time-scales, before neutrino irradiation can raise the electron fraction, enabling strong r-process nucleosynthesis up to and beyond the third peak. The resulting kilonova is lanthanide-rich ($X_{\rm lan} \approx 8~{{\ \rm per\ cent}}$) and dominated by near-infrared emission. We compute synthetic light curves in the Large Synoptic Survey Telescope and J ames Webb Space Telescope bands and find striking agreement, without parameter tuning, between the observations of AT 2023vfi/GRB 230307A and our broadband light curves for polar viewing angles. These results establish magnetized AIC as a viable channel for heavy r-process element production and a compelling progenitor candidate for long-duration gamma-ray bursts with kilonova signatures.
Periodic modulations of energetic electron fluxes observed throughout distant radiation zone by Geiger counters aboard satellites suggest MHD WAVE activity
Solar flare phenomena are examined in an introduction for advanced undergraduate and graduate physics students. Chapters are devoted to the history of observations, flare spectroscopy, flare magnetohydrodynamics, flare plasma physics, radiative processes in the solar plasma, preflare conditions, the impulsive phase, the gradual phase, and coronal mass ejections. Diagrams, graphs, and photographs are provided.
Theoretical models, observations and measurements have preoccupied astrophysicists for many centuries. Only in recent years, has the theory of relativity as applied to astrophysical flows met the challenges of how the governing equations can be solved numerically with accuracy and efficiency. Even without the effects of relativity, the physics of magnetohydrodynamic flow instability, turbulence, radiation, and enhanced transport in accretion disks has not been completely resolved. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks and also in the study of Gamma-Ray bursts (GRB). Thus, our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flowfield-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for Computational Relativistic Astrophysics (CRA) are demonstrated.
A study is made of the effect of a temperature-dependent Coulomb-like resistivity on the planar tearing mode. The local evolution of the temperature is described by an energy equation which includes Joule heating and optically thin radiation. The resulting system of coupled linear magnetohydrodynamic equations is solved numerically, and eigenfunctions and growth rates are obtained. In the absence of radiation, there are two distinct solutions above a critical value of the magnetic Reynolds number S, a tearing-like mode and a Joule-heating mode. Below this point, the growth rates coalesce into a conjugate-complex pair. When stable radiation (dR/dT greater than 0) is added, the heating mode disappears and a modified tearing excitation exists to much lower values of S before its growth is cut off by Ohmic heating. Examples are given for solar coronal parameters, and for those characteristic of fusion-research devices. The introduction of an effective value for the resistivity, in the presence of energy transport, allows a simple qualitative discussion of the different modes.
Discharge- and laser-produced plasma (DLPP) devices are being used as light sources for extreme ultraviolet (EUV) generation. A key challenge for both, DPP and LPP, is achieving sufficient brightness to support the throughput requirements of nanometrology tools. To simulate the environment of a hybrid DLPP device and optimize the EUV output, we have developed an integrated HEIGHTS-DLPP computer simulation package. The package integrates simulation of two evolving plasmas (DPP and LPP) and includes modeling of a set of integrated self-consistent processes: external power source and plasma energy balance, plasma resistive magnetohydrodynamics, plasma heat conduction, detailed radiation transport (RT), and laser absorption and refraction. We simulated and optimized DLPP devices using Xe gas as a target material. We synchronized the external circuit parameters, chamber gas parameters, and laser beam temporal and spatial profiles to achieve maximum EUV output. The full 3D Monte Carlo scheme was integrated for detailed RT and EUV output calculations in Xe using more than 3600 spectral groups. The modeling results are in good agreement with Julich Forschungszentrum experimental data. Theoretical models, developed and integrated into the HEIGHTS package, showed wide capabilities and flexibility. In conclusion, the models and package can be used for optimization of the experimental parameters and settings, investigation of DLPP devices with complex design, analyzing the impact of integrated spatial effects and working timeline arrangement on the final EUV output, and EUV source size, shape, and angular distribution.
A challenging aspect of the description of a tokamak disruption is evaluating the hot tail runaway electron seed that emerges during the thermal quench. This problem is made challenging due to the requirement of describing a strongly non-thermal electron distribution, together with the need to incorporate a diverse range of multiphysics processes, including magnetohydrodynamic instabilities, impurity transport, and radiative losses. Here this work develops a physics-informed neural network (PINN) tailored to the solution of the hot tail seed during an axisymmetric thermal quench. Here, a PINN is developed to identify solutions to the adjoint relativistic Fokker–Planck equation in the presence of a rapid quench of the plasma's thermal energy. It is shown that the PINN is able to accurately predict the hot tail seed across a range of parameters, including the thermal quench timescale, initial plasma temperature, and local current density, in the absence of experimental or simulation data. The hot tail PINN is verified by comparison with direct Monte Carlo simulations, with excellent agreement found across a broad range of thermal quench conditions.
We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic, and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers, we describe (1) a new solver for the Einstein equations based on the Z4c formalism, and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging Advanced RISC Machine processors. AthenaK shows excellent performance and weak scaling, achieving over 1 billion cell updates per second for hydrodynamics in three dimensions on a single NVIDIA Grace Hopper processor. It does this with a typical parallel efficiency of 80% on 65,536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics.
Shock wave structure in ionized hydrogen, using Eddington approximation to solve Rankine- Hugoniot jump equation
Gaseous breakdown of simulated planetary atmospheres and illustration of Greenhouse effect, thermal radiation from ionized argon, liquid magnetohydrodynamic power conversion
High temporal resolution x-ray streak camera studies of micropinch formation in Cu hybrid x pinches reveal key plasma conditions. Analysis of Ne-like Cu lines indicate an average electron temperature of about 200 eV and 4.5×10 28 m -3 electron density. Here, the spectra suggest that the electron temperature jumps to about 1 keV, inferred from the continuum and the postcontinuum line emission that includes Li-like Cu lines. There is no sign of a rapid temperature change or a substantial surge in radiation emission during the 200 ps precontinuum x-ray burst, suggesting that the radiative collapse process does not play a major role in micropinch formation. Two-dimensional extended Magnetohydrodynamic (MHD) simulations, coupled to a collisional-radiative spectral analysis code, suggest the significance of the rapid radial implosion of high-temperature, low-density plasma, the axial outflow, and the dynamic plasma pressure in micropinch formation.
In many black hole (BH) systems, the accretion disk is expected to be misaligned with respect to the BH spin axis. If the scale height of the disk is much smaller than the misalignment angle, the spin of the BH can tear the disk into multiple, independently precessing "sub-disks." This is most likely to happen during outbursts in black hole X-Ray binaries (BHXRBs) and in active galactic nuclei (AGNs) accreting above a few percent of the Eddington limit, because the disk becomes razor-thin. Disk tearing has the potential to explain variability phenomena including quasi-periodic oscillations in BHXRBs and changing-look phenomena in AGNs. Here, we present the first radiative two-temperature general relativistic magnetohydrodynamic (GRMHD) simulation of a strongly tilted (65°) accretion disk around an M BH = 10 M ⊙ BH, which tears and precesses. This leads to luminosity swings between a few percent and 50% of the Eddington limit on sub-viscous timescales. Surprisingly, even where the disk is radiation-pressure-dominated, the accretion disk is thermally stable over t ≳ 14,000 r g /c. This suggests warps play an important role in stabilizing the disk against thermal collapse. The disk forms two nozzle shocks perpendicular to the line of nodes where the scale height of the disk decreases tenfold and the electron temperature reaches T e ~ 10 8 –10 9 K. In addition, optically thin gas crossing the tear between the inner and outer disk gets heated to T e ~ 108 K. This suggests that warped disks may emit a Comptonized spectrum that deviates substantially from idealized models.
The composition and magnetic field morphology of relativistic jets can be studied using circular polarization (CP) measurement. Recent three-dimensional relativistic magnetohydrodynamic (3D RMHD) simulations coupled with radiative transfer (RT) calculations make strong predictions about the level (and morphology) of the jet’s CP emission. These simulations show that the sign of CP and the electric vector position angle (EVPA) are both sensitive to the jet’s magnetic field morphology within the radio core. We probe this theory by exploring whether the jet’s radio core EVPA orientation is consistent with the observed sign of the core CP in deep full-track polarimetric observations. Based on a selection of sources from earlier MOJAVE observations, we aim to probe the nature of linear polarization (LP) and CP in the innermost regions of jets from a small sample of nine blazars. This sample includes sources that have exhibited: (i) positive CP; (ii) negative CP; or (iii) positive and negative CP simultaneously in the radio core region. By coupling deep polarimetric observations of a carefully selected sample of blazars with state-of-the-art RMHD and RT calculations, we hope to gain a deeper understanding of the physics of blazar jets.