Search NASA⌕ Search

SEARCH · Search NASA

Results for “Collisional Radiative Modeling”

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 91 records · Page 5

Collisional balance of the meteoritic complex

The present study has the objective to reevaluate the size distribution of interplanetary meteoroids on the basis of the most recent data, and to analyze the probable nature of the sinks and sources of meteoritic material. The flux of interplanetary meteorites at 1 AU is discussed, taking into account general characteristics, lunar crater distribution, flux curves, spatial densities, and cross-sectional distribution and light scattering. Collisional effects are examined, giving attention to catastrophic collisions, collision rate, and destroyed mass and generated fragments. The effect of radiation pressure on small particles is considered along with the difference between the lunar and interplanetary flux models, collisional evolution at 1 AU, potential sources for large meteoroids, and observational evidence of losses of small micrometeoroids.

Gruen, E.↗

A Collisional Algorithm for Modeling Circumstellar Debris Disks

Many planetary systems harbor circumstellar disks of dust and planetesimals thought to be debris left over from planet formation. These debris disks exhibit a range of morphological features which can arise from the gravitational perturbations of planets. Accurate models of these features, accounting for the interactions of the particles in a disk with each other and with whatever planets they contain, can act as signposts for planets in debris disks that otherwise could not be detected. Such models can also constrain the planet's mass and orbital parameters. Current models for many disks consider the gravitational and radiative effects of the star and planets on the disk, but neglect the morphological consequences of collisional interactions between the planetesimals. Many observed disk features are not satisfactorily explained by the current generation of models. I am developing a new kind of debris disk model that considers both the gravitational shaping of the disk by planets and the inelastic collisions between particles. I will use a hybrid N-body integrator to numerically solve the equations of motion for the particles and planets in the disk. To include the collisional effects, I begin with an algorithm that tests for collisions at each step of the orbit integration and readjusts the velocities of colliding particles. I am adapting this algorithm to the problem at hand by allowing each particle to represent a "swarm" of planetesimals with a range of masses. When the algorithm detects an encounter between swarms, two or three swarms are produced to approximate the range of possible trajectories of the daughter planetesimals. Here I present preliminary results from my collisional algorithm.

Nesvold, Erika↗

(abstract)Electron Impact Emission Cross Sections for Modeling UV Auroral and Dayglow Observations of the Upper Atmospheres of Planets

In the upper atmospheres of the Jovian and Terrestrial planets a dominant mechanism for energy transfer occurs through electron collisional processes with neutral species leading to UV radiation. In response to the need for accurate collision cross sections to model spectroscopic observations of planetary systems, JPL has measured in the laboratory emission cross sections and medium resolution spectra of H, H(sub 2), N(sub 2), SO(sub 2), and other important planetary gases.Voyager and International Ultraviolet Explorer (IUE) spacecraft have established that band systems of H(sub 2) and N(sub 2) are the dominant UV molecular emissions in the solar system produced by electron impact. Applications of our data to models of Voyager, IUE, Galileo, and Hubble Space Telescope observations of the planets will be described.

molecular emissions↗

Stellar Laboratories: 3. New Ba 5, Ba 6, and Ba 7 Oscillator Strengths and the Barium Abundance in the Hot White Dwarfs G191-B2B and RE 0503-289

Context. For the spectral analysis of high-resolution and high-signal-to-noise (S/N) spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. Aims. Reliable Ba 5-7 oscillator strengths are used to identify Ba lines in the spectra of the DA-type white dwarf G191−B2B and the DO-type white dwarf RE 0503−289 and to determine their photospheric Ba abundances. Methods. We newly calculated Ba v-vii oscillator strengths to consider their radiative and collisional bound-bound transitions in detail in our NLTE stellar-atmosphere models for the analysis of Ba lines exhibited in high-resolution and high-S/N UV observations of G191−B2B and RE 0503−289. Results. For the first time, we identified highly ionized Ba in the spectra of hot white dwarfs. We detected Ba vi and Ba vii lines in the Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum of RE 0503−289. The Ba vi/Ba vii ionization equilibrium is well reproduced with the previously determined effective temperature of 70 000 K and surface gravity of log g=7.5. The Ba abundance is 3.5 +/- 0.5 × 10(exp−4) (mass fraction, about 23 000 times the solar value). In the FUSE spectrum of G191−B2B, we identified the strongest Ba vii line (at 993.41 Å) only, and determined a Ba abundance of 4.0 +/- 0.5 × 10(exp−6) (about 265 times solar). Conclusions. Reliable measurements and calculations of atomic data are a pre-requisite for stellar-atmosphere modeling. Observed Ba vi-vii line profiles in two white dwarfs' (G191−B2B and RE 0503−289) far-ultraviolet spectra were well reproduced with our newly calculated oscillator strengths. This allowed to determine the photospheric Ba abundance of these two stars precisely.

New Ba V↗

The oxygen emission lines in Wolf-Rayet spectra

The statistical equilibrium of oxygen in the line-emitting regions of WR stars is studied using the one-representative point theory of Castor and van Blerkom (1970) for a wider range of the radiation temperature, electron temperature, and electron density. Radiative processes, collisional processes, and electronic recombination are studied for a 59-level model atom. The level populations depart by large factors from those expected for LTE at the electron temperature and electron density. When the parameters have the values found previously by Bhatia and Underhill to be representative for WR stars, the oxygen atoms are strongly concentrated in the lower levels of O(2+), O(3+), and O(4+). The predicted relative energies in typical lines of O III, O IV, and O V for a solar abundance ratio of O to H are compatible with what is observed in WC and WN spectra.

Bhatia, A. K.↗

Particle-in-cell simulations of expanding high energy density plasmas with laser ray tracing

The design and analysis of high energy density (HED) laser experiments typically rely on radiation hydrodynamics simulations. However, some laser–plasma interaction regimes are not collisional and cannot be adequately modeled with hydrodynamics. For example, strongly driven magnetic reconnection and magnetized collisionless shock experiments possess extended hydrodynamic or even kinetic properties, necessitating first-principles kinetic simulations. In this paper, we present the benchmarking and first results obtained with a laser-ray-tracing and inverse bremsstrahlung absorption module implemented in the particle-in-cell code PSC. The simulation results are compared to radiation hydrodynamic simulations using the FLASH code as well as analytical estimates. We successfully benchmark the energy deposition model and overall hydrodynamic evolution of the systems. We also consider possible kinetic effects that may be expected from laser-target ablation in the HED regime, including non-local transport and two-temperature effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stellar Laboratories II. New Zn Iv and Zn v Oscillator Strengths and Their Validation in the Hot White Dwarfs G191-B2B and RE0503-289

Context. For the spectral analysis of high-resolution and high-signal-to-noise (SN) spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. In a recent analysis of the ultraviolet (UV) spectrum of the DA-type white dwarf G191B2B,21 Zn iv lines were newly identified. Because of the lack of Zn iv data, transition probabilities of the isoelectronic Ge vi were adapted for a first, coarse determination of the photospheric Zn abundance.Aims. Reliable Zn iv and Zn v oscillator strengths are used to improve the Zn abundance determination and to identify more Zn lines in the spectra of G191B2B and the DO-type white dwarf RE 0503289. Methods. We performed new calculations of Zn iv and Zn v oscillator strengths to consider their radiative and collisional bound-bound transitions in detail in our NLTE stellar-atmosphere models for the analysis of the Zn iv v spectrum exhibited in high-resolution and high-SN UV observations of G191B2B and RE 0503289. Results. In the UV spectrum of G191B2B, we identify 31 Zn iv and 16 Zn v lines. Most of these are identified for the first time in any star. We can reproduce well almost all of them at log Zn 5.52 0.2 (mass fraction, about 1.7 times solar). In particular, the Zn iv Zn v ionization equilibrium, which is a very sensitive Teff indicator, is well reproduced with the previously determined Teff 60 000 2000 K and log g 7.60 0.05. In the spectrum of RE 0503289, we identified 128 Zn v lines for the first time and determined log Zn 3.57 0.2 (155 times solar). Conclusions. Reliable measurements and calculations of atomic data are a pre-requisite for stellar-atmosphere modeling. Observed Zn iv and Zn v line profiles in two white dwarf (G191B2B and RE 0503289) ultraviolet spectra were well reproduced with our newly calculated oscillator strengths. This allowed us to determine the photospheric Zn abundance of these two stars precisely.

strengths↗

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Atmospheric Response to High Nonthermal Electron Beam Fluxes in Solar Flares. I. Modeling the Brightest NUV Footpoints in the X1 Solar Flare of 2014 March 29

The 2014 March 29 X1 solar flare (SOL20140329T17:48) produced bright continuum emission in the far- and near-ultraviolet (NUV) and highly asymmetric chromospheric emission lines, providing long-sought constraints on the heating mechanisms of the lower atmosphere in solar flares. We analyze the continuum and emission line data from the Interface Region Imaging Spectrograph (IRIS) of the brightest flaring magnetic footpoints in this flare. We compare the NUV spectra of the brightest pixels to new radiative-hydrodynamic predictions calculated with the RADYN code using constraints on a nonthermal electron beam inferred from the collisional thick-target modeling of hard X-ray data from Reuven Ramaty High Energy Solar Spectroscopic Imager. We show that the atmospheric response to a high beam flux density satisfactorily achieves the observed continuum brightness in the NUV. The NUV continuum emission in this flare is consistent with hydrogen (Balmer) recombination radiation that originates from low optical depth in a dense chromospheric condensation and from the stationary beam-heated layers just below the condensation. A model producing two flaring regions (a condensation and stationary layers) in the lower atmosphere is also consistent with the asymmetric Fe II chromospheric emission line profiles observed in the impulsive phase.

Kowalski, Adam F.↗

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λ c , that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λ c -based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan↗

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗

Dielectronic recombination rates, ionization equilibrium, and radiative emission rates for calcium and nickel ions in low-density high-temperature plasmas

The total dielectronic recombination rates for Ca and Ni ions were calculated taking into account autoionization to excited states of the recombining ion and stabilizing radiative transitions of the recombining electron. Radiative transitions of the recombining ion are found to be the dominant stabilizing processes, and the relative importance of Delta n = 0 and Delta n not equal to 0 transitions in the relevant temperature region is determined for each ion. The relative importance of the Delta n = 0 and Delta n not equal to 0 contributions is significantly altered for some ions by the inclusion of autoionization to excited levels. The relative abundance of the various ionization stages has been determined by using a corona equilibrium model in which collisional ionization and inner-shell excitation followed by autoionization are balanced by direct radiative and dielectronic recombination.

Jacobs, V. L.↗

VIII. New Zr IV-VII, Xe IV-V, and Xe VII Oscillator Strengths and the Al, Zr, and Xe Abundances in the Hot White Dwarfs G191-B2B and RE0503-289

Context: For the spectral analysis of high-resolution and high-signal-to-noise spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. Aims: To search for zirconium and xenon lines in the ultraviolet (UV) spectra of G191−B2B and RE 0503−289, new Zr iv-vii, Xe iv-v, and Xe vii oscillator strengths were calculated. This allows, for the first time, determination of the Zr abundance in white dwarf (WD) stars and improvement of the Xe abundance determinations.Methods: We calculated Zr iv-vii, Xe iv-v, and Xe vii oscillator strengths to consider radiative and collisional bound-bound transitions of Zr and Xe in our NLTE stellar-atmosphere models for the analysis of their lines exhibited in UV observations of the hot WDs G191−B2B and RE 0503−289. Results: We identified one new Zr iv, 14 new Zr v, and ten new Zr vi lines in the spectrum of RE 0503−289. Zr was detected for the first time in a WD. We measured a Zr abundance of −3.5 +/- 0.2 (logarithmic mass fraction, approx. 11 500 times solar). We identified five new Xe vi lines and determined a Xe abundance of −3.9 +/- 0.2 (approx. 7500 times solar). We determined a preliminary photospheric Al abundance of −4.3 +/- 0.2 (solar) in RE 0503−289. In the spectra of G191−B2B, no Zr line was identified. The strongest Zr iv line (1598.948 Å) in our model gave an upper limit of −5.6 +/- 0.3 (approx. 100 times solar). No Xe line was identified in the UV spectrum of G191−B2B and we confirmed the previously determined upper limit of −6.8 +/- 0.3 (ten times solar). Conclusions: Precise measurements and calculations of atomic data are a prerequisite for advanced NLTE stellar-atmosphere modeling. Observed Zr iv-vi and Xe vi-vii line profiles in the UV spectrum of RE 0503−289 were simultaneously well reproduced with our newly calculated oscillator strengths.

Rauch, T.↗

Collisional-radiative switching - A powerful technique for converging non-LTE calculations

A very simple technique has been developed to converge statistical equilibrium and model atmospheric calculations in extreme non-LTE conditions when the usual iterative methods fail to converge from an LTE starting model. The proposed technique is based on a smooth transition from a collision-dominated LTE situation to the desired non-LTE conditions in which radiation dominates, at least in the most important transitions. The proposed approach was used to successfully compute stellar models with He abundances of 0.20, 0.30, and 0.50; Teff = 30,000 K, and log g = 2.9.

Hummer, D. G.↗

Stellar Laboratories: New GeV and Ge VI Oscillator Strengths and their Validation in the Hot White Dwarf RE0503-289

State-of-the-art spectral analysis of hot stars by means of non-LTE model-atmosphere techniques has arrived at a high level of sophistication. The analysis of high-resolution and high-S/N spectra, however, is strongly restricted by the lack of reliable atomic data for highly ionized species from intermediate-mass metals to trans-iron elements. Especially data for the latter has only been sparsely calculated. Many of their lines are identified in spectra of extremely hot, hydrogen-deficient post-AGB stars. A reliable determination of their abundances establishes crucial constraints for AGB nucleosynthesis simulations and, thus, for stellar evolutionary theory. Aims. In a previous analysis of the UV spectrum of RE 0503-289, spectral lines of highly ionized Ga, Ge, As, Se, Kr, Mo, Sn, Te, I, and Xe were identified. Individual abundance determinations are hampered by the lack of reliable oscillator strengths. Most of these identified lines stem from Ge V. In addition, we identified Ge VI lines for the first time. We calculated Ge V and Ge VI oscillator strengths in order to reproduce the observed spectrum. Methods. We newly calculated Ge V and Ge VI oscillator strengths to consider their radiative and collisional bound-bound transitions in detail in our non-LTE stellar-atmosphere models for the analysis of the Ge IV-VI spectrum exhibited in high-resolution and high-S/N FUV (FUSE) and UV (ORFEUS/BEFS, IUE) observations of RE 0503-289. Results. In the UV spectrum of RE 0503-289, we identify four Ge IV, 37 Ge V, and seven Ge VI lines. Most of these lines are identified for the first time in any star. We can reproduce almost all Ge IV, GeV, and Ge VI lines in the observed spectrum of RE 0503-289 (T(sub eff) = 70 kK, log g = 7.5) at log Ge = -3.8 +/- 0.3 (mass fraction, about 650 times solar). The Ge IV/V/VI ionization equilibrium, that is a very sensitive T(sub eff) indicator, is reproduced well. Conclusions. Reliable measurements and calculations of atomic data are a prerequisite for stellar-atmosphere modeling. Our oscillator-strength calculations have allowed, for the first time, Ge V and Ge VI lines to be successfully reproduced in a white dwarf s (RE 0503-289) spectrum and to determine its photospheric Ge abundance.

Rauch, T.↗

Vibrational and rotational excitation of CO in comets. Part 1: Non-equilibrium calculations. Part 2: Results of the calculation for standard bright comet, comet Iras-Araki-Alcock and comet Halley

The vibrational and rotational excitation of the CO molecule in cometary comae were investigated using a model which includes IR vibrational pumping by the solar flux, vibrational and rotational radiative decay, and collisional coupling among rotational states. Steady state was not assumed in solving the rate equations. The evolution of a shell of CO gas was monitored as it expanded from the nucleus into the outer coma. Collisional effects were treated using a kinetic temperature profile derived from theoretical work on the coma energy balance. The kinetic temperature was assumed to be extremely cold in the inner coma; this has significant consequences for the CO excitation. If optical depth effects are ignored, only low J transitions will be significantly excited in comets observed at high spatial resolution. Ground-based observations of CO co-vibrational and rotational transitions will be extremely difficult due to lack of sensitivity and/or terrestrial absorption. However, CO should be detectable from a large comet with favorable observing geometry if the CO is a parent molecule present at the 10% level (or greater) relative to H2O. Observations using cooled, spaceborne instruments should be capable of detecting CO emission from even moderately bright comets.

Chin, G.↗

A hydrodynamic treatment of the tilted cold dark matter cosmological scenario

A standard hydrodynamic code coupled with a particle-mesh code is used to compute the evolution of a tilted cold dark matter (TCDM) model containing both baryonic matter and dark matter. Six baryonic species are followed, with allowance for both collisional and radiative ionization in every cell. The mean final Zel'dovich-Sunyaev y parameter is estimated to be (5.4 +/- 2.7) x 10 exp -7, below currently attainable observations, with an rms fluctuation of about (6.0 +/- 3.0) x 10 exp -7 on arcmin scales. The rate of galaxy formation peaks at a relatively late epoch (z is about 0.5). In the case of mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum peaks in the vicinity of 10 exp 9.1 solar masses, with a reasonable fit to the Schechter luminosity function if the baryon mass to blue light ratio is about 4. It is shown that a bias factor of 2 required for the model to be consistent with COBE DMR signals is probably a natural outcome in the present multiple component simulations.

Cen, Renyue↗