skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A Reduced-order NLTE Kinetic Model for Radiating Plasmas of Outer Envelopes of Stellar Atmospheres

Abstract

The present work proposes a self-consistent reduced-order NLTE kinetic model for radiating plasmas found in the outer layers of stellar atmospheres. A detailed collisional-radiative kinetic mechanism is constructed by leveraging the most up-to-date set of ab initio and experimental data available in the literature. This constitutes the starting point for the derivation of a reduced-order model, obtained by lumping the bound energy states into groups. In order to determine the needed thermo-physical group properties, uniform and Maxwell–Boltzmann energy distributions are used to reconstruct the energy population of each group. Finally, the reduced set of governing equations for the material gas and the radiation field is obtained based on the moment method. Applications consider the steady flow across a shock wave in partially ionized hydrogen. The results clearly demonstrate that adopting a Maxwell–Boltzmann grouping allows, on the one hand, for a substantial reduction of the number of unknowns and, on the other, to maintain accuracy for both gas and radiation quantities. Also, it is observed that, when neglecting line radiation, the use of two groups already leads to a very accurate resolution of the photo-ionization precursor, internal relaxation, and radiative cooling regions. The inclusion of line radiation requires adopting just onemore » additional group to account for optically thin losses in the α , β , and γ lines of the Balmer and Paschen series. This trend has been observed for a wide range of shock wave velocities.« less

Authors:
 [1];  [2];  [3]
  1. Aerospace Engineering Department, University of Illinois at Urbana-Champaign, 206A Talbot Lab., 104 S. Wright Street, Urbana, IL 61801 (United States)
  2. NASA Ames Research Center, Moffett Field, 94035 CA (United States)
  3. Aerospace Engineering Department, University of Illinois at Urbana-Champaign, 306 Talbot Lab., 104 S. Wright Street, Urbana, IL 61801 (United States)
Publication Date:
OSTI Identifier:
22661195
Resource Type:
Journal Article
Resource Relation:
Journal Name: Astrophysical Journal; Journal Volume: 838; Journal Issue: 2; Other Information: Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; ACCURACY; ENERGY SPECTRA; EQUATIONS; HYDROGEN; IONIZATION; LAYERS; LOSSES; PLASMA; RADIANT HEAT TRANSFER; RADIATIVE COOLING; REDUCTION; RELAXATION; RESOLUTION; SHOCK WAVES; STARS; STEADY FLOW; STELLAR ATMOSPHERES

Citation Formats

Munafò, Alessandro, Mansour, Nagi N., and Panesi, Marco, E-mail: munafo@illinois.edu, E-mail: nagi.n.mansour@nasa.gov, E-mail: m.panesi@illinois.edu. A Reduced-order NLTE Kinetic Model for Radiating Plasmas of Outer Envelopes of Stellar Atmospheres. United States: N. p., 2017. Web. doi:10.3847/1538-4357/AA602E.
Munafò, Alessandro, Mansour, Nagi N., & Panesi, Marco, E-mail: munafo@illinois.edu, E-mail: nagi.n.mansour@nasa.gov, E-mail: m.panesi@illinois.edu. A Reduced-order NLTE Kinetic Model for Radiating Plasmas of Outer Envelopes of Stellar Atmospheres. United States. doi:10.3847/1538-4357/AA602E.
Munafò, Alessandro, Mansour, Nagi N., and Panesi, Marco, E-mail: munafo@illinois.edu, E-mail: nagi.n.mansour@nasa.gov, E-mail: m.panesi@illinois.edu. Sat . "A Reduced-order NLTE Kinetic Model for Radiating Plasmas of Outer Envelopes of Stellar Atmospheres". United States. doi:10.3847/1538-4357/AA602E.
@article{osti_22661195,
title = {A Reduced-order NLTE Kinetic Model for Radiating Plasmas of Outer Envelopes of Stellar Atmospheres},
author = {Munafò, Alessandro and Mansour, Nagi N. and Panesi, Marco, E-mail: munafo@illinois.edu, E-mail: nagi.n.mansour@nasa.gov, E-mail: m.panesi@illinois.edu},
abstractNote = {The present work proposes a self-consistent reduced-order NLTE kinetic model for radiating plasmas found in the outer layers of stellar atmospheres. A detailed collisional-radiative kinetic mechanism is constructed by leveraging the most up-to-date set of ab initio and experimental data available in the literature. This constitutes the starting point for the derivation of a reduced-order model, obtained by lumping the bound energy states into groups. In order to determine the needed thermo-physical group properties, uniform and Maxwell–Boltzmann energy distributions are used to reconstruct the energy population of each group. Finally, the reduced set of governing equations for the material gas and the radiation field is obtained based on the moment method. Applications consider the steady flow across a shock wave in partially ionized hydrogen. The results clearly demonstrate that adopting a Maxwell–Boltzmann grouping allows, on the one hand, for a substantial reduction of the number of unknowns and, on the other, to maintain accuracy for both gas and radiation quantities. Also, it is observed that, when neglecting line radiation, the use of two groups already leads to a very accurate resolution of the photo-ionization precursor, internal relaxation, and radiative cooling regions. The inclusion of line radiation requires adopting just one additional group to account for optically thin losses in the α , β , and γ lines of the Balmer and Paschen series. This trend has been observed for a wide range of shock wave velocities.},
doi = {10.3847/1538-4357/AA602E},
journal = {Astrophysical Journal},
number = 2,
volume = 838,
place = {United States},
year = {Sat Apr 01 00:00:00 EDT 2017},
month = {Sat Apr 01 00:00:00 EDT 2017}
}
  • We present results of our calculations of NLTE model stellar atmospheres for hot Population III stars composed of hydrogen and helium. We use our own computer code for the calculation of spherically symmetric NLTE model atmospheres in hydrostatic and radiative equilibrium. The model atmospheres are then used for the calculation of emergent fluxes. These fluxes serve to evaluate the flow of high-energy photons for energies higher than ionization energies of hydrogen and helium, the so-called ionizing photon fluxes. We also present the time evolution of the ionizing photon fluxes.
  • The purpose of this work is to continue development of a model to provide a fast and accurate in-line NLTE capability for calculating plasma spectral properties in large-scale radiation-transport hydrodynamic simulations. A method has recently been developed to transform the large detailed atomic models into very small models that can be used for fast in-line calculations. The reduced model is more accurate than the average-atom models conventionally used in such simulations. In the present work, spectra calculated with the reduced model are compared to the original detailed model and the average-atom model. The spectra of iron and gold plasmas undermore » various plasma conditions are compared.« less
  • We present a high order kinetic flux-vector splitting (KFVS) scheme for the numerical solution of a conservative interface-capturing five-equation model of compressible two-fluid flows. This model was initially introduced by Wackers and Koren (2004) . The flow equations are the bulk equations, combined with mass and energy equations for one of the two fluids. The latter equation contains a source term in order to account for the energy exchange. We numerically investigate both one- and two-dimensional flow models. The proposed numerical scheme is based on the direct splitting of macroscopic flux functions of the system of equations. In two spacemore » dimensions the scheme is derived in a usual dimensionally split manner. The second order accuracy of the scheme is achieved by using MUSCL-type initial reconstruction and Runge-Kutta time stepping method. For validation, the results of our scheme are compared with those from the high resolution central scheme of Nessyahu and Tadmor . The accuracy, efficiency and simplicity of the KFVS scheme demonstrate its potential for modeling two-phase flows.« less
  • We investigate possible wind models for late-type stars which appear not to have hot coronae and transition regions. Taking Arcturus as our prototypical star, we consider wind models with T> or approx. =20,000 K and search for solutions with mass loss rates approx.10/sup -9/ M/sub sun/ yr/sup -1/. Thermally driven models which are spherically symmetric or include widely diverging geometries predict mass loss rates orders of magnitude less than 10/sup -9/ M/sub sun/ yr/sup -1/. We find that the radiation pressure of L..cap alpha.. resonance scattering can exceed the force of gravity in the chromosphere and initiate a net outflow,more » but it is insufficient to sustain a wind. If an additional momentum input term, such as Alfven wave pressure, is also present, then L..cap alpha.. radiation pressure may play a crucial role in turning on strong winds in stars like Arcturus by producing a critical point in the chromosphere and a locally supersonic flow which the additional mechanism can further accelerate. We conclude that L..cap alpha.. radiation-pressure-initiated winds can occur in stars to the right of the Linsky-Haisch dividing line in the H-R diagram between stars with and without transition regions and presumably hot coronae, and that the existence of these winds may explain energetically the absence of hot coronae in these stars.« less
  • The heavy element compositions of the solar corona, solar wind (SW), solar energetic particles (SEP, ''mass-unbiasedd baseline''; Paper I) and galactic cosmic-ray sources (GCRS) are remarkably similar. They all show the same pattern as compared to standard local galactic (or solar photospheric) composition: an underabundance of heavy elements with first ionization potential (FIP) > or approx. =9 eV relative to elements with lower FIP, by factors of approx.4-6. Only C is clearly more abundant in GCR sources than in SEP (by a factor of approx.2-3), as well as possibly O.