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.
Anomalous transport of multi-species plasma is considered with the generalized Hasegawa–Wakatani model. It is shown that the transport of all plasma species is described by fractional diffusion equations with the same effective diffusion coefficient. Strongly enhanced perturbations of heavy impurity density are found in long-living plasma flow vortices.
A theory of the three-wave coupling of triplets of tearing modes in toroidal pinches [i.e., either reversed field pinches (RFPs) or tokamaks] was proposed by R. Fitzpatrick [Phys. Plasmas 6, 1168 (1999)]. However, this theory only applies to toroidal pinches with negligible equilibrium plasma pressure gradients. Such a limitation is particularly inappropriate to RFPs. This paper generalizes the analysis of R. Fitzpatrick [Phys. Plasmas 6, 1168 (1999)] in order to take the equilibrium pressure gradient into account. However, for the sake of simplicity, a stepped pressure profile, rather than a continuous profile, is employed. In the limit in which the number of steps becomes very large, the results obtained from the generalized theory are presumably equivalent to those that would have been achieved using a continuous pressure profile. The generalized theory is used to investigate the formation of the characteristic toroidally localized pattern of phase-locked m = 1 and m = 0 tearing modes in RFP plasmas that is known as the “slinky” pattern. The incorporation of the equilibrium plasma pressure into the analysis is found to be of crucial importance when determining the properties of the pattern. This is the case because the plasma pressure controls the number of unstable m = 1 and m = 0 tearing modes, and also significantly affects the strength of three-wave coupling, as well as the phase relation between the phase-locked m = 1 and m = 0 modes.
The next-generation magnetic recoil spectrometer (MRSnext) is being designed to replace the current MRS at the National Ignition Facility and OMEGA for measurements of the neutron spectrum from an inertial confinement fusion implosion. The MRSnext will provide a far-superior performance and faster data turnaround than the current MRS systems, i.e., a 2× and 6× improvement in energy resolution at the NIF and OMEGA, respectively, and 20× improvement in data turnaround time. The substantially improved performance of the MRSnext is enabled by using electromagnets that provide a short focal plane (12–16 cm) and unprecedented flexibility for a wide range of applications. In addition to being able to measure neutron yield, apparent ion temperature, areal density, and plasma-flow velocity over a wide range of yields, the NIF MRSnext will be able to directly, uniquely assess the alpha heating of the fuel ions through measurements of the alpha knock-on tail in the neutron spectrum. The goal is to implement a radiation-hard electronic detection system capable of providing rapid data acquisition and analysis. The development of the MRSnext will also set the foundation for the more advanced, time-resolving MRSt and serve as a testbed for its implementation on the NIF.
The tearing mode stability of an inverse aspect-ratio expanded tokamak plasma equilibrium of general shape is investigated using asymptotic matching techniques. Particular emphasis is placed on the conservation of toroidal electromagnetic angular momentum. The TJ code, which is a specific implementation of the results of the investigation, is described.
Two-dimensional particle-in-cell simulations in the (radial-azimuthal) plane perpendicular to the axial direction of a cylindrical $E$ x $B$ Penning discharge are presented. The low-pressure discharge is self-consistently supported by plasma ionization from the electron beam injected axially, along the direction of the external magnetic field. It is shown that with the increasing strength of the external magnetic field, the discharge undergoes a sequence of transitions between several azimuthal modes. Azimuthal m > 1 spiral arm structures are excited at low magnetic field values as plasma confinement improves and the radial density profile becomes peaked. With a larger field, spiral arms with m > 1 are replaced by the m = 1 spoke mode, most clearly seen in plasma density. A transition from spiral arms to the spoke regime occurs when the plasma potential in the center changes from weakly positive (or zero) to negative. Further increase in the magnetic field results in a well-developed m = 1 spoke mode with additional small-scale higher-frequency m > 1 structures inside and around the spoke. It is shown that while ionization and collisions affect some characteristics of the observed fluctuations, the basic features of the spoke and m > 1 spiral structure remained similar without ionization. The role of energy conservation in small-scale high-frequency modes and spoke dynamics is discussed. It is demonstrated that in regimes with the m = 1 spoke mode, additional m = 4 harmonics of the ion and electron fluxes to the wall appear due to the square boundary. The frequency of the m = 1 mode is weakly affected by the geometry of the boundary.
Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.
Here a technique to identify intrinsic error fields (EFs) in tokamaks with minimized risk of disruption is demonstrated on the DIII-D tokamak. The method extends the conventional driven magnetic island ‘compass scan’ approach by modifying asynchronous control waveforms to enable prompt healing of the island instability. Healing of the island is achieved by reducing the imposed non-axisymmetric coil current and raising the density (here via gas fueling). The method is also shown to support multiple island threshold measurements per pulse, thus reducing the number of dedicated pulses necessary to conduct an EF identification. Non-linear modeling with the TM1 code reproduces the experimental results and approximately recovers the critical density required for island healing. Island healing is explained in the non-linear modeling by an increase in the viscous coupling between the static island and the nearby flowing plasma, thus healing the island as it accelerates into the plasma frame. Due to both simplicity and risk minimization, this technique is suitable for plasma-based EF identification in the early commissioning stages of future disruption-averse tokamaks such as ITER and SPARC.
Abstract In this work, the evolution of a highly unstable m = 1 resistive tearing mode, leading to plasmoid formation in a Harris sheet, is studied in the framework of full MHD model using the Non-Ideal Magnetohydrodynamics with Rotation, Open Discussion simulation. Following the initial nonlinear growth of the primary m = 1 island, the X-point develops into a secondary elongated current sheet that eventually breaks into plasmoids. Two distinctive viscous regimes are found for the plasmoid formation and saturation. In the low viscosity regime (i.e. P r ≲ 1), the plasmoid width increases sharply with viscosity, whereas in the viscosity dominant regime (i.e. P r ≳ 1), the plasmoid size gradually decreases with viscosity. Such a finding quantifies the role of viscosity in modulating the plasmoid formation process through its effects on the plasma flow and the reconnection itself.
In this work, we perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (~1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the HandOff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to harm3d, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the HandOff approach.
We present the first statistical study of magnetic structures and associated energy dissipation observed during a single period of turbulent magnetic reconnection, by using the in situ measurements of the Magnetospheric Multiscale mission in the Earth's magnetotail on 26 July 2017. The structures are selected by identifying a bipolar signature in the magnetic field and categorized as plasmoids or current sheets via an automated algorithm which examines current density and plasma flow. The size of the plasmoids forms a decaying exponential distribution ranging from subelectron up to ion scales. The presence of substantial number of current sheets is consistent with a physical picture of dynamic production and merging of plasmoids during turbulent reconnection. The magnetic structures are locations of significant energy dissipation via electric field parallel to the local magnetic field, while dissipation via perpendicular electric field dominates outside of the structures. Significant energy also returns from particles to fields.
Predicting, responding to, or interpreting the chemical record preserved in debris derived from nuclear events can be challenging due to chemical fractionation. Chemical fractionation is where different species of the evolving radionuclide inventory segregate and/or are lost from the system over the timescales of debris formation. Both historic data and recent research suggest that the interaction and character of the local environment may exert controls on chemical fractionation by influencing the cooling and evolution of the associated fireball as well as the composition of the vapor term and resultant speciation. Prior to this work, an integrated platform permitting dynamic and concurrent consideration of physical and chemical evolution of early time post-detonation event environments did not exist. Our work merged historic data and experimental approaches to support development of a computational framework able to simulate fundamental processes (e.g., entrainment of local environment, oxidation chemistry, and cooling time scales) that may perturb the radionuclide inventory captured in post-detonation debris. Work with historic debris confirmed that entrained environmental material affect debris composition, structure, and radionuclide incorporation. Complementary work utilizing a readily controllable and tunable benchtop setup (a plasma flow reactor) simulated the late cooling of a nuclear fireball (e.g., T < 6000 K) and bounded the sensitivity of actinide speciation and particle size distribution to variations in oxygen concentration and cooling rates. Concurrent laser ablation and laser heating experiments were used to investigate the chemistry and physics of processes occurring in vaporized and/or rapidly heated actinides and other elements in the presence of oxygen. A more computationally efficient microphysical model was developed for predicting and evolving size distributions of particles forming from mixed vapor terms and simulating particle formation processes under a variety of extreme conditions. Continued study of historic nuclear event film confirmed that shockwave data and physics codes agree to within the uncertainty of the data. Good agreement was achieved for thermal emission from an airburst, however the paucity of low-temperature molecular opacity data for mixtures of air, bomb debris, entrained dirt, and water vapor complicate agreement for more elaborate scenarios. A multiphysics code (ALE3D) was modified to bring the necessary physics and chemistry, including these new data and insights, onto a single platform. Code development included improved initialization of large physical systems, modernization of chemistry capabilities, and modifications to enable inclusion of particle transport.
Fallout analysis is a crucial field of study as it provides insight into the structural interactions that occur between radiological components and their surroundings following a nuclear event. Due to the complexity of actual fallout samples, setups like the plasma flow reactor (PFR) at LLNL have been developed to synthesize fallout in a controlled environment. While the PFR allows for control of parameters such as initial temperature, cooling rate, and concentration of radiological species, controlling the concentration of water has remained a challenge as samples are introduced into the setup via aqueous solution. To that end, our goal with this project was to develop methodologies to (1) monitor and (2) control the presence of water in the PFR. We found that a commercially available FTIR and MCT detector was sufficient to monitor the presence of both condensed and vapor-phase water, while a desolvating nebulizer could be used to control the amount of water injected into the system. In total, the results of this feasibility study have expanded the ability of the flow reactor to investigate key sensitivities that influence fallout formation.
Neutral hydrogen atoms in solar plasma flow as energy source of ring current belt and auroral substorms
Effects of interplanetary magnetic fields on flow of plasma around moon
Imp-I observations of solar wind properties, including interaction with geomagnetic field and flow characteristics of positive ion component
Plasma stream conductivity and velocity measurement, using systems of primary coil and several secondary coils
Perturbations of interplanetary magnetic field by lunar wake compared with disturbed solar plasma flow
Variations in monopole and dipole antenna reactance related to ionospheric sounding rocket spin period and plasma flow distribution