Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma global models”

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.

At least 145 records · Page 8

Quantification of the effect of uncertainty on impurity migration in PISCES-A simulated with GITR

A Bayesian inference strategy has been used to estimate uncertain inputs to global impurity transport code (GITR) modeling predictions of tungsten erosion and migration in the linear plasma device, PISCES-A. This allows quantification of GITR output uncertainty based on the uncertainties in measured PISCES-A plasma electron density and temperature profiles (n e , T e ) used as inputs to GITR. The technique has been applied for comparison to dedicated experiments performed for high (4 × 10 22 m –2 s –1 ) and low (5 × 10 21 m –2 s –1 ) flux 250 eV He–plasma exposed tungsten (W) targets designed to assess the net and gross erosion of tungsten, and corresponding W impurity transport. The W target design and orientation, impurity collector, and diagnostics, have been designed to eliminate complexities associated with tokamak divertor plasma exposures (inclined target, mixed plasma species, re-erosion, etc) to benchmark results against the trace impurity transport model simulated by GITR. The simulated results of the erosion, migration, and re-deposition of W during the experiment from the GITR code coupled to materials response models are presented. Specifically, the modeled and experimental W I emission spectroscopy data for a 429.4 nm line and net erosion through the target and collector mass difference measurements are compared. Furthermore, the methodology provides predictions of observable quantities of interest with quantified uncertainty, allowing estimation of moments, together with the sensitivities to plasma temperature and density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of the dayside boundary region of the magnetosphere

The region near the dayside boundary of the magnetosphere is of particular interest mainly because it is the most likely site for solar wind energy transfer. Exactly where and how this transfer occurs is one of the basic problems of magnetospheric physics. The global nature of the energy transfer is addressed in this research. In order to test whether energy transfer affects the large scale plasma flow in the magnetosheath, global properties of the magnetosheath were analyzed and compared to models with no energy transfer. Also, global models of merging sites on the magnetopause were modeled with the aid of the magnetosheath models. Since energy transfer and intrinsic magnetosheath properties depend upon the orientation of the interplanetary magnetic field (IMF), it was essential to have an IMF monitor for these studies. In the course of using ISEE 3 in this role, statistical studies were performed to test the accuracy of using data from an orbit so far upstream. Global properties of the magnetosheath energetic ion population, magnetic field strength and orientation, and flow pattern, determined with ISEE data are summarized.

Crooker, N.↗

A model of global convection in Jupiter's magnetosphere

Voyager observations of Jupiter's magnetosphere are compared with the planetary wind model in which corotation must break down outside some Alfven critical radius and a centrifugally driven wind outflow must develop. It is found that the model does not agree with the observations. A new global convection model for the Jovian magnetosphere is proposed, based on models of quasi-stationary plasma convection in the earth's magnetosphere. The model predicts a substantial dawn-dusk asymmetry in the structure, dynamics, and plasma composition of the magnetopause and magnetosheath. The model also predicts a region of cross-tail flow in the nightside plasma sheet containing a substantial admixture of solar wind plasma.

Cheng, A. F.↗

A numerical code for a three-dimensional magnetospheric MHD equilibrium model

Two dimensional and three dimensional MHD equilibrium models were begun for Earth's magnetosphere. The original proposal was motivated by realizing that global, purely data based models of Earth's magnetosphere are inadequate for studying the underlying plasma physical principles according to which the magnetosphere evolves on the quasi-static convection time scale. Complex numerical grid generation schemes were established for a 3-D Poisson solver, and a robust Grad-Shafranov solver was coded for high beta MHD equilibria. Thus, the effects were calculated of both the magnetopause geometry and boundary conditions on the magnetotail current distribution.

Voigt, G.-H.↗

A nonlinear dynamical analogue model of geomagnetic activity

Consideration is given to the solar wind-magnetosphere interaction within the framework of deterministic nonlinear dynamics. An earlier dripping faucet analog model of the low-dimensional solar wind-magnetosphere system is reviewed, and a plasma physical counterpart to that model is constructed. A Faraday loop in the magnetotail is considered, and the relationship of electric potentials on the loop to changes in the magnetic flux threading the loop is developed. This approach leads to a model of geomagnetic activity which is similar to the earlier mechanical model but described in terms of the geometry and plasma contents of the magnetotail. The model is characterized as an elementary time-dependent global convection model. The convection evolves within a magnetotail shape that varies in a prescribed manner in response to the dynamical evolution of the convection. The result is a nonlinear model capable of exhibiting a transition from regular to chaotic loading and unloading. The model's behavior under steady loading and also some elementary forms of time-dependent loading is discussed.

Klimas, A. J.↗

An MHD model of the earth's magnetosphere

It is pointed out that the earth's magnetosphere arises from the interaction of the solar wind with the earth's geomagnetic field. A global magnetohydrodynamics (MHD) model of the earth's magnetosphere has drawn much attention in recent years. In this model, MHD equations are used to describe the solar wind interaction with the magnetosphere. In the present paper, some numerical aspects of the model are considered. Attention is given to the ideal MHD equations, an equation of state for the plasma, the model as an initial- and boundary-value problem, the shock capturing technique, computational requirements and techniques for global MHD modeling, a three-dimensional mesh system employed in the global MHD model, and some computational results.

Wu, C. C.↗

The topside ionosphere - A region of dynamic transition

The review article concentrates on dynamic processes at work in the topside ionosphere (between the F2 peak and about 3000 km) where the H ion dominates and ionic reactions can be neglected. The history of ionosphere and plasmasphere research using radio waves is reviewed. Low-speed and high-speed multispecies plasma ion flow is studied with various models (13-moment approximation, 5-moment approximation, kinetic models of the polar wind). Experimental observations of the plasmapause, results of vertical soundings of the topside, and global pole-to-pole distributions of ion composition, plasma temperature, and electron density are reviewed.

Banks, P. M.↗

Modeling of Field-Aligned Guided Echoes in the Plasmasphere

The conditions under which high frequency (f>>f(sub uh)) long-range extraordinary-mode discrete field-aligned echoes observed by the Radio Plasma Imager (RPI) on board the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) satellite in the plasmasphere are investigated by ray tracing modeling. Field-aligned discrete echoes are most commonly observed by RPI in the plasmasphere although they are also observed over the polar cap region. The plasmasphere field-aligned echoes appearing as multiple echo traces at different virtual ranges are attributed to signals reflected successively between conjugate hemispheres that propagate along or nearly along closed geomagnetic field lines. The ray tracing simulations show that field-aligned ducts with as little as 1% density perturbations (depletions) and less than 10 wavelengths wide can guide nearly field-aligned propagating high frequency X mode waves. Effective guidance of wave at a given frequency and wave normal angle (Psi) depends on the cross-field density scale of the duct, such that ducts with stronger density depletions need to be wider in order to maintain the same gradient of refractive index across the magnetic field. While signal guidance by field aligned density gradient without ducting is possible only over the polar region, conjugate field-aligned echoes that have traversed through the equatorial region are most likely guided by ducting.

Fung, Shing F.↗

Dayside Proton Aurora: Comparisons between Global MHD Simulations and Image Observations

The IMAGE mission provides a unique opportunity to evaluate the accuracy of current global models of the solar wind interaction with the Earth's magnetosphere. In particular, images of proton auroras from the Far Ultraviolet Instrument (FUV) onboard the IMAGE spacecraft are well suited to support investigations of the response of the Earth's magnetosphere to interplanetary disturbances. Accordingly, we have modeled two events that occurred on June 8 and July 28, 2000, using plasma and magnetic field parameters measured upstream of the bow shock as input to three-dimensional magnetohydrodynamic (MHD) simulations. This paper begins with a discussion of images of proton auroras from the FUV SI-12 instrument in comparison with the simulation results. The comparison showed a very good agreement between intensifications in the auroral emissions measured by FUV SI-12 and the enhancement of plasma flows into the dayside ionosphere predicted by the global simulations. Subsequently, the IMAGE observations are analyzed in the context of the dayside magnetosphere's topological changes in magnetic field and plasma flows inferred from the simulation results. Finding include that the global dynamics of the auroral proton precipitation patterns observed by IMAGE are consistent with magnetic field reconnection occurring as a continuous process while the iMF changes in direction and the solar wind dynamic pressure varies. The global simulations also indicate that some of the transient patterns observed by IMAGE are consistent with sporadic reconnection processes. Global merging patterns found in the simulations agree with the antiparallel merging model. though locally component merging might broaden the merging region, especially in the region where shocked solar wind discontinuities first reach the magnetopause. Finally, the simulations predict the accretion of plasma near the bow shock in the regions threaded by newly open field lines on which plasma flows into the dayside ionosphere are enhanced. Overall the results of these initial comparisons between global MHD simulation results and IMAGE observations emphasize the interplay between reconnection and dynamic pressure processes at the dayside magnetopause. as well as the intricate connection between the bow shock and the auroral region.

Berchem, J.↗

SciDAC ISEP: Integrated Simulation of Energetic Particles in Burning Plasmas

The objective of the SciDAC Center for Integrated Simulation of Energetic Particles in Burning Plasmas (ISEP) is to improve physics understanding of energetic particle (EP) confinement and EP interactions with burning thermal plasmas through large-scale simulations. The ISEP center will develop a multiscale and multiphysics ISEP framework for a predictive capability of EP physics and deliver an EP module incorporating both first-principles simulations and high fidelity reduced transport models to the fusion whole device modeling (WDM) project. The ISEP framework will enable us to perform long time, global kinetic simulations of EP physics in burning plasmas, by utilizing the full power of the next generation supercomputers. Our research and development activities will build on fruitful collaborations with computer scientists and applied mathematicians to offer enabling technologies for performance scalability, portability, solvers, coupling for integration with the fusion WDM project, and long-term preservation of data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DE-1 observations of solar wind-magnetosphere coupling processes in the polar cusp

Measurements of hot plasma distribution functions in the midaltitude polar cusp by the Dynamics Explorer-1 spacecraft have provided new information on the processes of magnetosheath plasma injection and transport and the dependence of these phenomena on solar-wind parameters. Results on magnetosheath plasma injection have included: (1) localization of the region of plasma injection to the high-latitude cusp-magnetosheath interface, (2) movement of this region toward lower latitudes as B(z) becomes more southward, (3) maintenance of quasi-neutrality of suprathermal plasmas within the cusp, and (4) support for merging-type injection at the high-latitude cusp boundary for northward IMF. Plasma flow measurements from Dynamics Explorer-1 have led to a global IMF-dependent convection model which contains both merging and viscous driving forces.

Burch, J. L.↗

Global energy balance simulations of shattered pellet injection

The dynamics of shattered pellet injection (SPI) shutdowns are simulated using a time-dependent global energy balance model, based on a modification of the KPRAD framework. The new SPI particle source in the model calculates the ablation of individual pellet fragments that enter the plasma as a temporally resolved plume, thus capturing the effects of earlier fragments on the ablation of those that follow, which has a significant impact on the overall assimilation. Despite the reduced physics and the global averaging of all quantities, results from a large number of DIII-D, KSTAR, and JET experiments are well reproduced, including the plasma cooling timescales, particle assimilation, and current quench (CQ) rates. Cooling timescales and CQ rates are in good agreement for pellets containing as little as ∼1% neon, while particle assimilations are most accurate for neon fractions above ∼15% by number of atoms. Below this, the assimilation tends to be overestimated due to the lack of radial transport in the particle balance, which becomes important in the low-Z limit. Predictive simulations of mixed-composition dual-SPI shutdowns in ITER are compared against those with the 3D non-linear magnetohydrodynamic code JOREK and are found to reproduce overall trends observed in the higher-fidelity modeling across a range of injection scenarios. The general success of the model points to the critical role of energy balance in determining SPI particle assimilation and the subsequent disruption dynamics and highlights the value of these simulations for experimental interpretation and for optimizing the deployment of computationally expensive, higher-fidelity models.

energy balance model↗

High-Energy Pulsar Models: Developments and New Questions

The past few years have seen a major advance in observational knowledge of high-energy (HE) pulsars. The Fermi Large Area Telescope (LAT) and AGILE have increased the number of known gamma-ray pulsars by an order of magnitude, its members being divided roughly equally among millisecond pulsars (MSPs), young radio-loud pulsars, and young radio-quiet pulsars. Many new and diverse emission characteristics are being measured, while radio and X-ray follow-up observations increase the pulsar detection rate and enrich our multiwavelength picture of these extreme sources. The wealth of new data has provided impetus for further development and improvement of existing theoretical pulsar models. Geometric light curve (LC) modelling has uncovered three broad classes into which HE pulsars fall: those where the radio profile leads, is aligned with, or lags the gamma-ray profile. For example, the original MSP and original black widow system are members of the second class, requiring co-located emission regions and thereby breaking with traditional notions of radio emission origin. These models imply narrow accelerator gaps in the outer magnetosphere, indicating copious pair production even in MSP magnetospheres that were previously thought to be pair-starved. The increased quality and variety of the LCs necessitate construction of ever more sophisticated models. We will review progress in global magnetosphere solutions which specify a finite conductivity on field lines above the stellar surface, filling the gap between the standard vacuum and force-free (FF; plasma-filled) models. The possibility of deriving phase-resolved spectra for the brightest pulsars, coupled with the fact that the HE pulsar population is sizable enough to allow sampling of various pulsar geometries, will enable much more stringent testing of future radiation models. Reproduction of the observed phase-resolved behavior of this disparate group will be one of the next frontiers in pulsar science, impacting on our understanding of particle acceleration, emission, and magnetosphere geometry. One may now also study evolutionary trends of the measured or inferred quantities, and probe pulsar visibility and population properties such as radiation beam sizes of different pulsar classes, as well as the distribution of spin-down power, gamma-ray luminosity, conversion efficiency, spectral index, and cutoff energy across the population. Lastly, the recent detection of very-high-energy (VHE) pulsations from the Crab pulsar generated quite a few ideas to explain this emission, leading to an extension of standard models and possibly even a bridge between the physics of pulsars and pulsar wind nebulae (PWNe).

gamma rays:↗

On the termination of the closed field line region of the magnetotail

Experimental models and global computer simulations are used to study the field structure and location of the site where the closed plasma sheet terminates. The site - also called the 'distant neutral line' - is described with references to several solar-wind parameters that can change its characteristic properties. A laboratory experiment is then conducted in which a flowing magnetized plasma beam interacts with a dipole field. High-latitude nightside reconnection is supported by the simulation results comparing the effects of northward and southward interplanetary fields. The results support the models by Dungey (1961, 1963) that predict high-latitude nightside reconnection and an open magnetosphere. A significant finding is that a magnetic neutral line at the closed-field line region is not a necessary component of the region.

Birn, J.↗

Helicity Injected Torus (HIT): Sustainment of Stable Equilibria Through Self-organization

This is the Final Progress Report of the Helicity Injected Torus (HIT) program covering May 2016 through December 2020. During this final period the HIT program continued to make advancements in computational modeling, theory, and experimental studies of global self-organization and inductively sustained spheromaks. This included many valuable contributions to the fusion, plasma and physics community, through peer-reviewed publications, PhD theses and conference presentations and posters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Potential Applications of Quantum Computing at Los Alamos National Laboratory, v0.3.0

Since the scientific revolution in the 16th and 17th centuries, the process of scientific discovery has followed an iterative feedback process of observation, hypothesis development and testing with physical experiments, which is widely referred to as the scientific method. This process remained largely unchanged until the middle of the 20th century, when the emergence of digital computers empowered scientist to build and inspect detailed simulations of physical phenomena. Over the last century, computational tools have transformed modern approaches to scientific discovery by enabling fast and affordable hypothesis testing before physical experiments are conducted, shown in Figure 1-1. Some notable examples include: global climate forecasts to understand how the environment may change over decades [130]; modeling the behavior of plasma to design fusion reactors [59]; and understanding the behavior of molecules in biological processes [161, 223].

36 MATERIALS SCIENCE↗

On the Importance of Using Event-Specific Wave Diffusion Rates in Modeling Diffuse Electron Precipitation

A few to tens of keV electron precipitation that carries substantial energy source down to the upper atmosphere to create aurora is manifested as an important magnetosphere-ionosphere coupling process. The precipitation is usually caused by scattering processes associated with plasma waves in the magnetosphere. The scattering process is often quantified by wave diffusion rates that indicate how fast an electron is scattered. Global models commonly use diffusion coefficients that are derived from statistical wave models. However, due to the statistical nature, many localized, transient features could be smeared out. In this study, we investigate electron precipitation using event-specific diffusion coefficients that are obtained based on simultaneous in-situ measured/inferred, rather than statistical, chorus wave dynamics. We find that the application of the event-specific diffusion coefficients associated with a more dynamic and intense chorus wave model leads more electrons, particularly at several to tens of keV in the dawn-to-noon sector at L > 3, to precipitate than using statistical coefficients. Here, the new simulation roughly captures both the intensity and variability of the precipitating flux as detected by the NOAA/POES satellites. Ionospheric electron density in the lower E region (100–120 km) observed by the mid-latitude Millstone Hill radar is also much better reproduced, while the case using statistical diffusion coefficients underestimates the ionization rate. This study implies the importance of using event-specific diffusion rates in simulating the diffuse electron precipitation and understanding the magnetosphere-ionosphere coupling.

79 ASTRONOMY AND ASTROPHYSICS↗

Reflection and transmission of GMIR shock at the heliopause and their relation to the 2- and 3-kHz radio emissions

We use Voyager 2 plasma and magnetic field data together with a one-fluid MHD model to study the interactions of the 1991 Global Merged Interaction Region (GMIR) shock with the heliopause. The 1991 GMIR is an extraordinarily large global solar wind structure in radial, longitudinal and latitudinal extents. It has a strong shock at the leading edge. After its penetration through the termination shock, the GMIR shock first propagates through the subsonic solar wind, then interacts with the heliopause. The interaction produces a transmitted shock propagating outward in the interstellar medium, and a reflected shock propagating backward in the subsonic solar wind. We identify the reflected shock and the transmitted shock as the possible source of the radio noise detected at Voyagers. The plasma frequency behind the reflected and the transmitted shock can be, respectively, responsible for the 2- and 3-kHz radio emissions. The two bands of radio noise are emitted from sources on both sides of the heliopause starting at about the same time. If the emission is generated by f(sub p)-radiation then the heliopause is located at R approximately 130 AU. If the emission is generated by 2f(sub p)-radiation the n R approximately 150 AU. Because the relative speed of the interstellar plasma with respect to the sun appears to be sub-Alfvenic, it is very unlikely there is a fast-bow shock of the heliosphere.

Whang, Y. C.↗