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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Specifications of a Plasmasphere Modeling Code for GGCM

The Dynamic Global Core Plasma Model (DGCPM) is a parameterized model for core or thermal plasma in the magnetosphere. The model accounts for dayside ionospheric outflow and nightside inflow. It accounts for the global pattern of convection and corotation. The model is capable of being coupled to ring current and superthermal electron models for the purpose of providing thermal plasma spatial distributions and for the purpose of accepting the dynamic influences of these plasma populations back upon the thermal plasma. The DGCPM is designed to operate alone or to operate as part of a larger integrated package. The convection electric field and magnetic field used within the DGCPM can be shared with models of other plasma populations, in addition to the exchange of parameters important to the collective modeling of whole plasma systems in the inner magnetosphere. This talk will present the features of the DGCPM model code and the various forms of information that can be exchanged with other cooperating codes.

Gallagher, D. L.↗
NASA NTRS2006

Comparison of Two IRI plasmasphere Extensions with GPS-TEC Observations

Comparisons of two model results with Global Positioning System GPS-TEC measurements have been carried out for different latitudinal, solar activity, magnetic activity, diurnal and seasonal conditions. The models evaluated are the Global Core Plasma Model (GCPM-2000) and the IRI extension with Russian plasmasphere model (IRI*).Data of 23 observatories providing GPS-TEC and ionosonde data have been used. It is shown that IRI* plasmasphere electron density is greater than GCPM results by an order of magnitude at 6370 km altitude (one Earth's radius) with this excess growing to 2-3 orders of magnitude towards the GPS satellite altitude of 20000 km. Another source of model and GPS-TEC differences is a way of selection of the F2 layer peak parameters driving the models either with ITU-R (former CCIR) maps or ionosonde observations. Plasmasphere amendment to IRI improves accuracy of TEC model predictions because the plasmasphere contribution to the total TEC varies from 10% by daytime under quiet magnetic conditions to more than 50% by night under stormy conditions.

Gulyaeva, T. L.↗
NASA NTRS2006

Analyzing Electric Field Morphology Through Data-Model Comparisons of the GEM IM/S Assessment Challenge Events

The storm-time inner magnetospheric electric field morphology and dynamics are assessed by comparing numerical modeling results of the plasmasphere and ring current with many in situ and remote sensing data sets. Two magnetic storms are analyzed, April 22,2001 and October 21-23,2001, which are the events selected for the Geospace Environment Modeling (GEM) Inner Magnetosphere/Storms (IM/S) Assessment Challenge (IMSAC). The IMSAC seeks to quantify the accuracy of inner magnetospheric models as well as synthesize our understanding of this region. For each storm, the ring current-atmosphere interaction model (RAM) and the dynamic global core plasma model (DGCPM) were run together with various settings for the large-scale convection electric field and the nightside ionospheric conductance. DGCPM plasmaspheric parameters were compared with IMAGE-EUV plasmapause extractions and LANL-MPA plume locations and velocities. RAM parameters were compared with Dst*, LANL-MPA fluxes and moments, IMAGE-MENA images, and IMAGE-HENA images. Both qualitative and quantitative comparisons were made to determine the electric field morphology that allows the model results to best fit the plasma data at various times during these events. The simulations with self-consistent electric fields were, in general, better than those with prescribed field choices. This indicates that the time-dependent modulation of the inner magnetospheric electric fields by the nightside ionosphere is quite significant for accurate determination of these fields (and their effects). It was determined that a shielded Volland-Stern field description driven by the 3-hour Kp index yields accurate results much of the time, but can be quite inconsistent. The modified Mcllwain field description clearly lagged in overall accuracy compared to the other fields, but matched some data sets (like Dst*) quite well. The rankings between the simulations varied depending on the storm and the individual data sets, indicating that each field description did well for some place, time, and energy range during the events, as well as doing less well in other places, times, and energies. Several unresolved issues regarding the storm-time inner magnetospheric electric field are discussed.

Liemohn, Michael W.↗
NASA NTRS1998

Convection of Plasmaspheric Plasma into the Outer Magnetosphere and Boundary Layer Region: Initial Results

We present initial results on the modeling of the circulation of plasmaspheric- origin plasma into the outer magnetosphere and low-latitude boundary layer (LLBL), using a dynamic global core plasma model (DGCPM). The DGCPM includes the influences of spatially and temporally varying convection and refilling processes to calculate the equatorial core plasma density distribution throughout the magnetosphere. We have developed an initial description of the electric and magnetic field structures in the outer magnetosphere region. The purpose of this paper is to examine both the losses of plasmaspheric-origin plasma into the magnetopause boundary layer and the convection of this plasma that remains trapped on closed magnetic field lines. For the LLBL electric and magnetic structures we have adopted here, the plasmaspheric plasma reaching the outer magnetosphere is diverted anti-sunward primarily along the dusk flank. These plasmas reach X = -15 R(sub E) in the LLBL approximately 3.2 hours after the initial enhancement of convection and continues to populate the LLBL for 12 hours as the convection electric field diminishes.

Ober, Daniel M.↗
NASA NTRS1998

Convection of Plasmaspheric Plasma into the Outer Magnetosphere and Boundary Layer Region: Initial Results

We present initial results on the modeling of the circulation of plasmaspheric-origin plasma into the outer magnetosphere and low-latitude boundary layer (LLBL), using a dynamic global core plasma model (DGCPM). The DGCPM includes the influences of spatially and temporally varying convection and refilling processes to calculate the equatorial core plasma density distribution throughout the magnetosphere. We have developed an initial description of the electric and magnetic field structures in the outer magnetosphere region. The purpose of this paper is to examine both the losses of plasmaspheric-origin plasma into the magnetopause boundary layer and the convection of this plasma that remains trapped on closed magnetic field lines. For the LLBL electric and magnetic structures we have adopted here, the plasmaspheric plasma reaching the outer magnetosphere is diverted anti-sunward primarily along the dusk flank. These plasmas reach X= -15 R(sub E) in the LLBL approximately 3.2 hours after the initial enhancement of convection and continues to populate the LLBL for 12 hours as the convection electric field diminishes.

Ober, Daniel M.↗
NASA NTRS2008

Multispectral Emission of the Sun during the First Whole Sun Month: Magnetohydrodynamic Simulations

We demonstrate that a three-dimensional magnetohydrodynamic (MHD) simulation of the corona can model its global plasma density and temperature structure with sufficient accuracy to reproduce many of the multispectral properties of the corona observed in extreme ultraviolet (EW) and X-ray emission. The key ingredient to this new type of global MHD model is the inclusion of energy transport processes (coronal heating, anisotropic thermal conduction, and radiative losses) in the energy equation. The calculation of these processes has previously been confined to one-dimensional loop models, idealized two-dimensional computations, and three-dimensional active region models. We refer to this as the thermodynamic MHD model, and we apply it to the time period of Carrington rotation 1913 (1996 August 22 to September 18). The form of the coronal heating term strongly affects the plasma density and temperature of the solutions. We perform our calculation for three different empirical heating models: (1) a heating function exponentially decreasing in radius; (2) the model of Schrijver et al.; and (3) a model reproducing the heating properties of the quiet Sun and active regions. We produce synthetic emission images from the density and temperature calculated with these three heating functions and quantitatively compare them with observations from E W Imaging Telescope on the Solar and Heliospheric Observatory and the soft X-ray telescope on Yohkoh. Although none of the heating models provide a perfect match, heating models 2 and 3 provide a reasonable match to the observations.

Lionello, Roberto↗
DOE OSTI2021

Turbulent field fluctuations in gyrokinetic and fluid plasmas

A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modeling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Turbulent field fluctuations in gyrokinetic and fluid plasmas

A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modelling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework is readily adaptable to experimental and astrophysical environments, and presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

AToM: Advanced Tokamak Modeling Environment

Stability constraints play an important role in integrated modeling. The global plasma stability and macro-instabilities such as internal kink modes, neoclassical tearing modes, and edge localized perturbations limit the plasma performance and can result in large-scale transient events and plasma disruptions. These macro-instabilities can be also beneficial to the plasma performance. For example, the peeling models often lead to the edge localized modes (ELMs). However, the formation of stationary edge harmonic oscillations (EHOs) due to nonlinear interaction of peeling modes can result in a transition to the Quiescent H-mode (QH-mode) without ELMs. The profiles predicted with transport models need to be refined using the MHD stability constraints. This is important especially for transient and nonlinear stages of discharges such as ramp-up access to hybrid and steady state operation, and L- to H-mode transition as well as transition to the QH-mode. Several stability and MHD codes are already included in the OMFIT framework. These codes include BALOO, ELITE, GATO, M3DC1, MARS, and NIMROD. However, the verifications of stability constraints are currently mostly excluded from transport modeling workflows. The only exception is the EPED module which includes the ELITE predictions to limit the pedestal height. Here, we improved the OMFIT workflow to include stability calculations. Including the stability conditions to the integrated modeling workflow improved the robustness of the predictive modeling discharges. Being implemented in the workflow, the stability conditions can be also used for the stability analysis of experimental data. This improved the physics understanding of various discharge scenarios and can be used in experiment planning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2005

Comparison of Two IRI Plasmasphere Extensions with GPS-TEC Observations

Two plasmasphere extensions of the International Reference Ionosphere are made available for the users. It is aimed to estimate the effect of charged particles on technical devices in the Earth's environment and to define the ionosphere-plasmasphere operational conditions compatible with existing and future systems of radio communication, radio navigation and other relevant radio technologies in the ranges of medium and higher frequencies. The Global Core Plasma Model (GCPM-2000) of Gallagher et al. (2000) is an empirical description of thermal plasma densities in the plasmasphere, plasmapause, magnetospheric trough, and polar cap. GCPM-2000 uses the Kp index and is coupled to IRI in the transition region 500-600 km. The IZMIRAN plasmasphere model (Chasovitin et al., 1998; Gulyaeva et al., 2002) is an empirical model based on whistler and satellite observations. It presents global vertical analytical profiles of electron density smoothly fitted to IRI electron density profile at 1000 km altitude and extended towards the plasmapause (up to 36,000 km). For the smooth fitting of the two models, the shape of the IRI topside electron density profile is improved using ISIS 1, ISIS 2, and IK19 satellite inputs (Gulyaeva, 2003). The plasmasphere model depends on solar activity and magnetic activity (kp-index). The two IRI plasmasphere extensions are compared in the present study with the total electron content derived from records of Global Positioning Satellites (GPS-TEC) observations for different latitudinal, solar activity, magnetic activity, diurnal and seasonal conditions. The differences of model TEC with observed TEC in the topside ionosphere and plasmasphere are discussed.

Gulyacva, Tamara↗
NASA NTRS1999

Empirical Modeling of the Plasmasphere

Over 40 years of ground and spacecraft plasmaspheric measurements have resulted in many statistical descriptions of plasmaspheric properties. In some cases, these properties have been represented as analytical descriptions that are valid for specific regions or conditions. For the most part, what has not been done is to extend regional empirical descriptions or models to the plasmasphere as a whole. In contrast, many related investigations depend on the use of representative plasmaspheric conditions throughout the inner magnetosphere. Wave propagation, involving the transport of energy through the magnetosphere, is strongly affected by thermal plasma density and its composition. Ring current collisional and wave particle losses also strongly depend on these quantities, Plasmaspheric also plays a secondary role in influencing radio signals from the Global Positioning System satellites. The Global Core Plasma Model (GCPM) is an attempt to assimilate previous empirical evidence and regional models for plasmaspheric density into a continuous, smooth model of then-nal plasma density in the inner magnetosphere. In that spirit, the International Reference Ionosphere is currently used to complete the low altitude description of density and composition in the model. The models and measurements on which the GCPM is currently based and its relationship to IRI will be discussed.

Gallagher, Dennis L.↗
NASA NTRS2001

Determining the Mass Density Along Magnetic Field Lines from Toroida Eigenfrequencies

Toroidal eigenfrequencies can be used to remotely sense the equatorial mass density rho(sub eq) and the density dependence along a magnetic field line. Here we present improvements to the method of Schulz [1996], which allows rho(sub eq) and the power law index alpha (for mass density along a field line proportional to R(sup -alpha), where R is the radial distance from the center of the Earth) to be determined from the y intercept and slope of a plot of toroidal frequency versus toroidal harmonic number n. Our modifications include a model form for eigenfrequencies with a fractional precision of 0.0005 for -6 less than or = alpha less than or = 6 and 2 less than or = L less than or = 8 (accuracy is doubtful beyond L = 5) and an iterative procedure for getting more accurate results than those found using Schulz's method. In addition, we do an analysis of the effect of random measurement errors. Observed frequencies need to be accurate to approx. 6% (3%) of the fundamental frequency in order to determine rho(sub eq) (alpha) to a precision of 30% (unity). We then apply our method to data generated using the Global Core Plasma Model for plasmaspheric mass density; our analysis demonstrates clearly bow the alpha index represents the mass density dependence on the outer part of the field line (R/(LR(sub E)) greater than or approx. 2/3).

Denton, R.↗
DOE OSTI2022

Design Considerations for Controlling Silicon Nanoparticle Nucleation and Growth in a Nonthermal Plasma

Controlling the nucleation and growth of nanoparticles in low temperature plasma systems is imperative for controlling nanoparticle size distributions; and for some applications such as deposition and etching in microelectronic processing, preventing particle contamination. In this work, silicon nanoparticle (NP) production from silane is used as a model system to investigate the nucleation process. Although the mechanisms responsible for silicon NP nucleation and growth have been studied, it is unclear how controllable system parameters (e.g., pressure, system geometry, and gas composition) can be used to inhibit or promote NP formation. For example, the transport of reactive silane species is expected to significantly affect the feed fraction of silane required to nucleate silicon NP (the nucleation onset fraction) due to losses at the reactor walls. In this work, NP mass density was determined as a function of system pressure, gas composition, and reactor diameter for a tubular flow-through radiofrequency (RF) plasma using Ar/H2/He/SiH4 gas mixtures. A quartz crystal microbalance (QCM) impactor was developed to measure the total aerosol mass density downstream of the plasma and thereby identify the nucleation onset and its dependence on process parameters. Finally, a reaction mechanism was developed and incorporated into a global plasma chemistry model to better understand the nucleation onset and NP growth.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A reaction mechanism for plasma electrolysis of AgNO 3 forming silver nanoclusters and nanoparticles

In plasma-driven solution electrolysis (PDSE), gas-phase plasma-produced species interact with an electrolytic solution to produce, for example, nanoparticles. An atmospheric pressure plasma jet (APPJ) directed onto a liquid solution containing a metallic salt will promote reduction of metallic ions in solution, generating metallic clusters that nucleate to form nanoparticles. In this article, results from a computational investigation are discussed of a PDSE process in which a radio-frequency APPJ sustained in helium impinges on a silver nitrate solution, resulting in growth of silver nanoparticles. A reaction mechanism was developed and implemented in a global plasma chemistry model to predict nanoparticle growth. To develop the reaction mechanism, density functional theory was used to generate probable silver growth pathways up to Ag 9 . Neutral clusters larger than Ag 9 were classified as nanoparticles. Kinetic reaction rate coefficients for thermodynamically favorable growth pathways were estimated based on an existing, empirically determined base reaction mechanism for smaller Ag particle interactions. These rates were used in conjunction with diffusion-controlled reaction rate coefficients that were calculated for other Ag species. The role of anions in reduction of Agn ions in forming nanoparticles is also discussed. Oxygen containing impurities or admixtures to the helium, air entrainment into the APPJ, and dissociation of saturated water vapor above the solution can produce additional reactive oxygen species in solution, resulting in the production of anions and O 2 – in particular. For a given molarity, delivering a sufficient fluence of reducing species will produce similar nanoparticle densities and sizes for all applied power levels. In conclusion, comparisons are made to alternate models for nanoparticle formation, including charged nanoparticles and use of direct current plasmas.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reaction mechanism for atmospheric pressure plasma treatment of cysteine in solution

Mechanisms for the cold atmospheric plasma (CAP) treatment of cells in solution are needed for more optimum design of plasma devices for wound healing, cancer treatment, and bacterial inactivation. However, the complexity of organic molecules on cell membranes makes understanding mechanisms that result in modifications (i.e. oxidation) of such compounds difficult. As a surrogate to these systems, a reaction mechanism for the oxidation of cysteine in CAP activated water was developed and implemented in a 0-dimensional (plug-flow) global plasma chemistry model with the capability of addressing plasma-liquid interactions. Reaction rate coefficients for organic reactions in water were estimated based on available data in the literature or by analogy to gas-phase reactions. The mechanism was validated by comparison to experimental mass-spectrometry data for COST-jets sustained in He/O 2 , He/H 2 O and He/N 2 /O 2 mixtures treating cysteine in water. Results from the model were used to determine the consequences of changing COST-jet operating parameters, such as distance from the substrate and inlet gas composition, on cysteine oxidation product formation. Results indicate that operating parameters can be adjusted to select for desired cysteine oxidation products, including nitrosylated products.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

A surface mechanism for O 3 production with N 2 addition in dielectric barrier discharges

Ozone, O 3 , is a strong oxidizing agent often used for water purification. O 3 is typically produced in dielectric barrier discharges (DBDs) by electron-impact dissociation of O 2 , followed by three-body association reactions between O and O 2 . Previous studies on O 3 formation in low-temperature plasma DBDs have shown that O 3 concentrations can drop to nearly zero after continued operation, termed the ozone-zero phenomenon (OZP). Including small (<4%) admixtures of N 2 can suppress this phenomenon and increase the O 3 production relative to using pure O 2 in spite of power deposition being diverted from O 2 to N 2 and the production of nitrogen oxides, N x O y . The OZP is hypothesized to occur because O 3 is destroyed on the surfaces in contact with the plasma. Including N 2 in the gas mixture enables N atoms to occupy surface sites that would otherwise participate in O 3 destruction. The effect of N 2 in ozone-producing DBDs was computationally investigated using a global plasma chemistry model. A general surface reaction mechanism is proposed to explain the increase in O 3 production with N 2 admixtures. The mechanism includes O 3 formation and destruction on the surfaces, adsorption and recombination of O and N, desorption of O 2 and N 2 , and NO x reactions. Without these reactions on the surface, the density of O 3 monotonically decreases with increasing N 2 admixture due to power absorption by N 2 leading to the formation of nitrogen oxides. With N-based surface chemistry, the concentrations of O 3 are maximum with a few tenths of percent of N 2 depending on the O 3 destruction probability on the surface. The consequences of the surface chemistry on ozone production are less than the effect of gas temperature without surface processes. An increase in the O 3 density with N-based surface chemistry occurs when the surface destruction probability of O 3 or the surface roughness was decreased.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

GloPW

This model predicts global plasma wave activity, with a focus on electromagnetic ion cyclotron (EMIC) waves in Earth’s magnetosphere.

Noh, Sung Jun↗
NASA NTRS1988

Modeling magnetospheric plasma; Proceedings of the First Huntsville Workshop on Magnetosphere/Ionosphere Plasma Models, Guntersville, AL, Oct. 14-16, 1987

The conference presents papers on the global modeling of magnetospheric plasma processes, the modeling of the midlatitude ionosphere and plasmasphere, the modeling of the auroral zone and boundary layer, the modeling of the polar magnetosphere and ionosphere, and the modeling of the plasma sheet and ring current. Particular attention is given to the kinetic approach in magnetospheric plasma transport modeling, self-consistent neutral point current and fields from single particle dynamics, preliminary statistical survey of plasmaspheric ion properties from observations by DE 1/RIMS, and a model of auroral potential structures based on dynamics explorer plasma data. Other topics include internal shear layers in auroral dynamics, quantitative parameterization of energetic ionospheric ion outflow, and open flux merging in an expanding polarcap model.

Moore, T. E.↗