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 199 records · Page 11

Time-Dependent Photodissociation Regions

We present theoretical models of the time-dependent thermal and chemical structure of molecular gas suddenly exposed to far-ultraviolet (FUV) (6 eV less than hv less than 13.6 eV) radiation fields and the consequent time- dependent infrared emission of the gas. We focus on the response of molecular hydrogen for cloud densities ranging from n = 10(exp 3) to 10(exp 6)/cu cm and FUV fluxes G(sub 0) = 10(exp 3)-10(exp 6) times the local FUV interstellar flux. For G(sub 0)/n greater than 10(exp -2) cu cm, the emergent H(sub 2) vibrational line intensities are initially larger than the final equilibrium values. The H(sub 2) lines are excited by FUV fluorescence and by collisional excitation in warm gas. Most of the H(sub 2) intensity is generated at a characteristic hydrogen column density of N approximately 10(exp 21)/sq cm, which corresponds to an FUV optical depth of unity caused by dust opacity. The time dependence of the H(sub 2) intensities arises because the initial abundances of H(sub 2) at these depths is much higher than the equilibrium values, so that H(sub 2) initially competes more effectively with dust in absorbing FUV photons. Considerable column densities of warm (T approximately 1000) K H(sub 2) gas can be produced by the FUV pumping of H(sub 2) vibrational levels followed by collisional de-excitation, which transfers the energy to heat. In dense (n greater than or approximately 10(exp 5)/cu cm) gas exposed to high (G(sub 0) greater than or approximately 10(exp 4)) fluxes, this warm gas produces a 2-1 S(1)/1-0 S(l) H(sub 2) line ratio of approximately 0.1, which mimics the ratio found in shocked gas. In lower density regions, the FUV pumping produces a pure-fluorescent ratio of approximately 0.5. We also present calculations of the time dependence of the atomic hydrogen column densities and of the intensities of 0 I 6300 A, S II 6730 A, Fe II 1.64 microns, and rotational OH and H20 emission. Potential applications include star-forming regions, clouds near active galactic nuclei, and planetary nebulae. We apply our models to five planetary nebulae and conclude that only BD +30deg3639 shows evidence of enhanced H(sub 2) emission due to (high) nonequilibrium H(sub 2) abundances.

Hollenbach, David↗

A performant energy-conserving particle reweighting method for Particle-in-Cell simulations

A new particle-based reweighting method is developed and demonstrated in the Aleph Particle-in-Cell with Direct Simulation Monte Carlo (PIC-DSMC) program. Novel splitting and merging algorithms ensure that modified particles maintain physically consistent positions and velocities. This method allows a single reweighting simulation to efficiently model plasma evolution over orders of magnitude variation in density, while accurately preserving energy distribution functions (EDFs). Demonstrations on electrostatic sheath and collisional rate dynamics show that reweighting simulations achieve accuracy comparable to fixed weight simulations with substantial computational time savings. This highly performant reweighting method is recommended for modeling plasma applications that require accurate resolution of EDFs or exhibit significant density variations in time or space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simulation of plasma and neutral transport in PISCES-RF using SOLPS-ITER

In this research, the fluid plasma transport code SOLPS-ITER is applied and validated against experimental data from the plasma interaction surface component experimental station (PISCES)-RF linear plasma device to establish a physics basis for plasma and neutral transport in its two magnetic field (B-field) geometry setups-(1)the cusp and (2) non-cusp or linear B-field. The main focus of this study is to understand (1) radial plasma transport (2) heat and particle loads on the upstream dump and downstream target plate, and (3) the physics of plasma-neutral interactions in PISCES-RF. The simulation setup adheres to typical PISCES-RF experimental conditions, with a 2D helicon power deposition profile as an input heating source. SOLPS-ITER simulations reproduce experimental conditions with the Bohm diffusion model for both B-field configurations of the PISCES-RF experiment. Major energy loss channels include neutral radiation and power deposited on the wall and dump plate, with only 1% of the input power reaching the target. The ionization front is well confined near the dump plate due to the heating and puffing regions. Additionally, SOLPS-ITER simulation results are also found to be in very good agreement with the particle-in-cell calculations using the code-PICOS++ which supports the validity of SOLPS in low collisionality regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hot-gas cold-dust pumping for water masers associated with H II regions

A collisional pump with an internal sink is proposed for the water masers associated with H II regions, where the population inversion occurs due to the absorption by cold ice-mantle grains in a highly dusty cloud of the far-infrared line radiation of hot water vapor. A new escape probability method is developed to calculate the transfer of line radiation in dusty medium. The pump mechanism explains the power of usual maser sources associated with H II regions and the enormous power of the sources associated with W49 N and external galaxies. Models of maser clouds have a radius of 5 x 10 to the 15th-10 to the 16th cm, an H2 number density of 4 x 10 to the 9th/cu cm, an expansion velocity of 10-30 km/s, a kinetic temperature of 350 K, and a grain temperature of 100 K. Giant maser sources require grains of the size about 1 micron. The apparent size of the emission spots (approximately 10 to the 13th cm) observed by VLBI is interpreted as due to a fluctuation in the cloud, and the assembly of the spots is spread within a size of 10 to the 16th cm. The temperature difference between the dust and gas is due to a relaxation process after an infrared burst accompanying protostar formation.

Deguchi, S.↗

A hydrodynamic approach to cosmology - Texture-seeded cold dark matter and hot dark matter cosmogonies

Hydrodynamical simulations of galaxy formation in a texture-seeded cosmology are presented, with attention given to Omega = 1 galaxies dominated by both hot dark matter (HDM) and cold dark matter (CDM). The simulations include both gravitational and hydrodynamical physics with a detailed treatment of collisional and radiative thermal processes, and use a cooling criterion to estimate galaxy formation. Background radiation fields and Zel'dovich-Sunyaev fluctuations are explicitly computed. The derived galaxy mass function is well fitted by the observed Schechter luminosity function for a baryonic M/L of 3 and total M/L of 60 in galaxies. In both HDM and CDM texture scenarios, the 'galaxies' and 'clusters' are significantly more strongly correlated than the dark matter due to physical bias processes. The slope of the correlation function in both cases is consistent with observations. In contrast to Gaussian models, peaks in the dark matter density distributrion are less correlated than average.

Cen, R. Y.↗

Nuclear-pumped He/3/-Ar laser modeling

A first-order model of He(3)-Ar 1.79-micron laser is developed, compared to experimental results, and used to explain the qualitative features of this system. Results indicate that direct excitation of the argon upper level is at best very inefficient for population inversion. For argon concentrations which give the most efficient laser operation, the He-3(n,p)He-3(n,p)He-3 energy is used to produce atomic He ions that quickly convert into He molecular ions. These molecular ions subsequently form argon atomic ions through charge transfer. The dominant pumping mechanism is collisional-radiative recombination of the argon atomic ion and subsequent radiative cascading into the upper laser level.

Wilson, J. W.↗

Applications of a Venus thermospheric circulation model

A variety of Pioneer Venus observations suggest a global scale, day-to-night Venus thermospheric circulation. Model studies of the dynamics and energetics of the Venus thermosphere are presented in order to address new driving, mixing and cooling mechanisms for an improved model simulation. The adopted approach was to reexamine the circulation by first using a previous two dimensional code to quantify those physical processes which can be inferred from the Pioneer Venus observations. Specifically, the model was used to perform sensitivity studies to determine the degree to which eddy cooling, eddy or wave drag, eddy diffusion and 15 micrometer radiational cooling are necessary to bring the model temperature and composition fields into agreement with observations. Three EUV heating cases were isolated for study. Global temperature and composition fields in good agreement with Pioneer data were obtained. Large scale horizontal winds 220 m/s were found to be consistent with the observed cold nightside temperatures and dayside bulges of O, CO and CO2. Observed dayside temperatures were obtained by using a 7 to 19% EUV heating efficiency profile. The enhanced 15 micrometer cooling needed for thermal balance is obtained using the best rate coefficient available for atomic O collisional excitation of CO2(0,1,0). Eddy conduction was not found to be a viable cooling mechanism due to the weakened global circulation. The strong 15 micrometer damping and low EUV efficiency imply a very weak dependence of the general circulation to solar cycle variability. The NCAR terrestrial thermospheric general circulation model was adapted for Venus inputs using the above two dimensional model parameters, to give a three dimensional benchmark for future Venus modelling work.

Bougher, S. W.↗

Stellar Energetic Particle-driven Production of Biologically Relevant Molecules in Atmospheres of Young Earth-like Exoplanets Around Active G-K Stars

The chemistry of N2, CO2, CH4-rich atmospheres of terrestrial-type exoplanets around active G-K stars is a complex problem: the star's ionizing radiation in the form of X-ray and Extreme UV radiation and the precipitating energetic particles accelerated in coronal mass ejection driven shocks can drive complex chemistry in the exoplanetary atmosphere, and ignite the production of complex molecules (Airapetian et al., 2016). We developed a set of atmospheric models to simulate the formation of organic molecules that can be considered as atmospheric pre-biosignatures from exoplanets around active stars. The model has recently been extended to account for the production of ions, their role in the enhanced chemistry in the lower atmosphere, as well as the interaction of these ions with the aerosols. We applied the photochemical-collisional model to simulate the chemistry of a young terrestrial-type exoplanets resembling the Hadean Earth to highlight the creation of atmospheric pre-biosignatures, relevant greenhouse gases, and the presence of aerosols driven by stellar activity sources. We also discuss the impact of that modified atmospheric chemistry on the exoplanetary climate.

Guillaume Philippe Gronoff↗

New analysis of the Voyager UVS H Lyman-alpha emission of Saturn

The limb to limb Lyman-alpha reflectivities observed with the Voyager ultraviolet spectrometer (UVS) instruments during the fly-by of Saturn are reanalyzed using a revised H Lyman-alpha sensitivity for the Voyager 1 instrument. The new sensitivity reconciles the measured intensities to those of Voyager 2 and gives a coherent set of data. To fit the UV airglow observations, four sources are considered: (1) H resonance and H2 Rayleigh scattering of solar Lyman-alpha radiation, (2) the interplanetary Lyman-alpha radiation, (3) a possible internal source of unknown origin, (4) the possibility of atmospheric turbulence recently proposed to explain the Lyman-alpha bulge of Jupiter. The analysis supports neither a dominant collisional excitation source for the UV emissions nor the presence of strong atmospheric turbulence. The best fit, in terms of brightness but also in terms of shape of the limb to limb profile (that is to say independent on the absolute calibrations), is obtained for pure resonance and Rayleigh scattering of solar and interstellar wind line in an atmosphere enriched in atomic hydrogen up to three times the standard model. Influx of water from the rings of Saturn may provide a means for producing such enhanced H densities in the upper atmosphere.

Jaffel, L. Ben↗

Diagnostics: Chapter 8 of the special issue: on the path to tokamak burning plasma operation

This chapter presents the activity conducted by the ITPA topical group (TG) on Diagnostics over about the last 15 years. Following a general introduction of the ITER Diagnostics led by their measurement roles, the document is organized in several subchapters detailing the design support, research and development activity conducted by each of the specialist working groups (WGs) of the TG. Please note that the magnetic diagnostics were supported at the TG without a specific WG. Their status is included in the general introduction. In the following some highlights of the subchapter’s contents are provided. Recent advances in ITER first wall (FW) diagnostics for the measurements of plasma-metallic wall interaction in support of the ITER research plan are reported. An InfraRed imaging Video Bolometer for ITER has been developed and tested on several tokamaks to measure the radiated power loss. A laser-induced breakdown spectroscopy (LIBS) technique which utilizes a pulsed laser beam to ablate locally by forming a crater, will measure local tritium inventory in the FW material. Real-time Residual Gas Analyzers will measure the neutral gas composition in a divertor port and an equatorial port during plasma operation. Due to the full metallic FW environment, the plasma-wall interaction in ITER will face several challenges such as the compromised radiated power and divertor heat flux measurements by reflection. Ray tracing and analysis codes have been developed to eliminate and correct the effects of reflection in the measurements. The characteristics of the reflecting surfaces depending on the roughness and angle of the incidence have been measured by dedicated experiments, and the results were applied to the reflection elimination. For the measurement of the metallic impurity radiation induced by eroded metallic atoms, a vacuum ultraviolet spectrometer has been developed and tested. An extensive thermonuclear diagnostic suite will be required to support the operation of ITER and the planned experimental program for future burning plasma experiments. Due to the harsh environmental conditions, the implementation of diagnostic systems in ITER is a major challenge. These conditions include high levels of neutron and gamma fluxes, neutron heating, particle bombardment. Therefore, the selection and design of diagnostic systems must take into account a number of phenomena previously unseen in diagnostic design. For this reason, the measurement of neutrons and confined or lost fast ions, with particular emphasis on alpha particles, is critical to ITER. The diagnostics associated with these measurements will be important for future plasma-burning experiments at ITER. The high neutron emission and very large plasma size in ITER make neutron diagnostics the main diagnostic method used to measure plasma parameters such as fusion power, fusion power density, ion temperature, energy of fast ions and their spatial distributions in the plasma core. Active spectroscopy techniques are methods where a neutral particle beam is injected into the plasma and information on plasma parameters is extracted from the measurement of line emission resulting from the beam-plasma interaction, either by plasma ions or by beam atoms. Spatial localization is achieved by crossing the beamline and multiple observation lines. The ITER plasma will be a high temperature, moderately dense, fully ionized collisional plasma. The plasma facing surfaces are principally metallic being fashioned from beryllium or tungsten but many other elements, arising from either structural or from operational needs, may enter this plasma. The energy range of the emitted photons range from meV (infra-red) to multi keV (x-rays) and originate from all areas of the plasma volume. The primary role of passive emission diagnostics is to identify what is in the plasma from spectral signatures. Extracting quantitative information from these measurements such as impurity content, ion temperature, rotation, degree of detachment and radiated power depends on calibrated instruments, a physics model of the atomic and molecular processes and plasma transport and an analysis workflow that takes into account environmental effects such as reflections. The particular needs for ITER have prompted a multi-machine, many-year effort to address all these aspects and this chapter reviews the work on diagnostic design, experiments and new analysis techniques. An overview of the laser diagnostics to be implemented on ITER is also provided in this paper. This includes descriptions of the Thomson scattering in the core, edge and divertor regions, polarimetry and interferometry diagnostics used for measuring plasma density and also measurements of helium density in the divertor using Laser Induced Flourescence. Techniques which can allow improvements on current measurements are also addressed in particular expanding poloidal polarimetry measurements to measure field fluctuations and proposed use of dispersion interferometery which has a number of advantages over existing methods. This paper identifies particular areas where further research and testing on existing tokamaks is useful even at this advanced stage to inform the design of diagnostics for ITER. Outstanding areas of concern for the implementation of laser diagnostics, in particular with a view to reliable operation are identified. An overview of the latest developments of microwave diagnostic systems and techniques is given. The primary focus is the contributions for ITER—the next step burning plasma experiment—which is supplemented by describing recent progress of techniques applicable for fusion experiments beyond ITER. The contributions are intentionally kept concise, and are being supplemented by a rich list of references for further studies. Radiation induced effects are receiving continuous and well-deserved attention of the ITER diagnostic community and they are in many cases one of the primary design drivers of the ITER diagnostic systems. The paper summarizes recent progress in this area focusing primarily on the ITER diagnostics but in some cases provides also outlook for the possible solutions for even more demanding radiation environment of fusion reactors beyond ITER. Despite advancements in the area of modeling and simulation of various radiation induced effects, experimental testing in a nuclear environment as close as possible to the target one is still seen as unavoidable for proper qualification of particular diagnostic functional elements. Recent advancement within three diagnostic areas: optical diagnostics, magnetics and bolometers is covered. Encouraging results on qualification of silica glass vacuum window assemblies are presented. In the area of magnetic sensors, progress of irradiation tests performed on ITER in-vessel LTCC inductive sensors is presented with outlook for novel technological approaches to inductive sensors utilizing thick printing and photolithography technologies being highlighted. Summary of advancements in the area of steady state magnetic field sensors based on Hall effect is given. New results of neutron irradiation test of the ITER borosilicate glass inserts for vacuum electrical feedthroughs are summarized finding negligible swelling at target level of neutron fluence. Off-line irradiation tests of fiber optic current sensors for plasma current measurement demonstrated that both for gamma doses up to 5 MGy and a total neutron fluence up to 10 15 cm −2 , radiation induced changes are still compatible with required measurement accuracy on ITER. The ITER bolometers are given as an example how considering radiation effects may influence the diagnostic design. Finally, outlook for future main R&D directions is outlined. All optical and laser-based diagnostics in ITER will be using mirrors to guide plasma radiation toward detectors, cameras and sensors. In the hostile plasma, radiation and particle environment the optical characteristics of diagnostic mirrors will degrade directly affecting the entire performance of involved diagnostic systems. An assessment of factors affecting mirror performance is provided. Among the prime adverse factors are deposition of plasma impurities, sputtering of mirror surface and steam ingress in the vicinity of mirrors. Within the International Tokamak Physics Activity with active support by ITER central team and domestic agencies, the structured research and development (R&D) program on mitigation of risks for diagnostic mirrors is underway. Within this program the mirror material development, the passive mitigation of mirror degradation by using diagnostic ducts and shutters along with an active mirror recovery program comprising the in-situ mirror cleaning and calibration is underway. Recent developments in diagnostic mirror R&D are described in this Chapter along with an example of their implementation of R&D solutions in ITER Infrared Thermography diagnostic. An assessment of still open engineering and physics questions, considerations on mirror risks during an early phase of ITER operation are given along with an overview of diagnostic mirror evolution in the late ITER operation stage toward the demonstration fusion power plant. Several crucial areas of diagnostic R&D outlined in ITER Research Plan are addressed. The basic control groups in a fusion reactor can be broken-down in five categories: (1) plasma position, magnetic configuration, and plasma current control, (2) profile control and confinement optimization, (3) MHD control and suppression, (4) edge dissipation control, radiation and plasma exhaust control and (5) break-down optimization. These categories are coupled via the physics (a control action in one domain will affect the other domains) and via shared actuators (e.g. ECRH for impurity accumulation avoidance, current density distribution control and MHD suppression). Consequently, a supervisory control system should determine the priority of the various control tasks, their couplings, and the interfaces with the safety and interlock system. For the systematic development of the various controllers taking the complexity of the plasma and the control system into account, a model-based approach is required. A short historical overview is given of the developments in systems and control theory and control engineering with special emphasis on those developments that are most relevant for Nuclear Fusion research and operation. An overview is given of the state of the field of fusion plasma control for the control categories. It will be shown how synthetic diagnostics are being developed in ITER and how they are used in diagnostic design and design validation and how they can be in model-based controller synthesis using relatively simple models. In modern control methods, multiple diagnostics are used to constrain relatively simple models. The constrained models provide an estimate for the state. This opens the route to state controllers, such as model predictive control. A major challenge in nuclear fusion research is the coherent combination of data from heterogeneous diagnostics and modeling codes for machine control and safety as well as physics studies. Measured data from different diagnostics often provide information about the same subset of physical parameters. Additionally, information provided by some diagnostics might be needed for the analysis of other diagnostics. A joint analysis of complementary and redundant data allows, e.g. to improve the reliability of parameter estimation, to increase the spatial and temporal resolution of profiles, to obtain synergistic effects, to consider diagnostics interdependencies and to find and resolve data inconsistencies. Physics-based modeling and parameter relationships provide additional information improving the treatment of ill-posed inversion problems. A coherent combination of all kind of available information within a probabilistic framework allows for improved data analysis results. The concept of integrated data analysis (IDA) in the framework of Bayesian probability theory is outlined and contrasted with conventional data analysis. Components of the probabilistic approach are summarized and specific ingredients beneficial for data analysis at fusion devices are discussed.

ITER↗

High-resolution Laboratory Measurements of K-shell X-ray Line Polarization and Excitation Cross Sections in Helium-like S XV Ions

We report measurements of electron-impact excitation cross sections for the strong K-shell n = 2 -> 1 transitions in S XV, using the LLNL EBIT-I electron beam ion trap, two crystal spectrometers, and the EBIT Calorimeter Spectrometer. The cross sections are determined by direct normalization to the well-known cross sections of radiative electron capture, measured simultaneously. Using contemporaneous polarization measurements with the two crystal spectrometers, whose dispersion planes are oriented parallel and perpendicular to the electron beam direction, the polarization of the direct excitation line emission is determined, and in turn the isotropic total cross sections are extracted. We further experimentally investigate various line-formation mechanisms, finding that radiative cascades and collisional inner-shell ionization dominate the degree of linear polarization and total line emission cross sections of the forbidden line, z.

atomic data↗

Detecting and Constraining N2 Abundances in Planetary Atmospheres Using Collisional Pairs

Characterizing the bulk atmosphere of a terrestrial planet is important for determining surface pressure and potential habitability. Molecular nitrogen (N2) constitutes the largest fraction of Earth's atmosphere and is likely to be a major constituent of many terrestrial exoplanet atmospheres. Due to its lack of significant absorption features, N2 is extremely difficult to remotely detect. However, N2 produces an N2-N2 collisional pair, (N2)2, which is spectrally active. Here we report the detection of (N2)2 in Earth's disk-integrated spectrum. By comparing spectra from NASA's EPOXI mission to synthetic spectra from the NASA Astrobiology Institute's Virtual Planetary Laboratory three-dimensional spectral Earth model, we find that (N2)2 absorption produces a ~35% decrease in flux at 4.15 μm. Quantifying N2 could provide a means of determining bulk atmospheric composition for terrestrial exoplanets and could rule out abiotic O2 generation, which is possible in rarefied atmospheres. To explore the potential effects of (N2)2 in exoplanet spectra, we used radiative transfer models to generate synthetic emission and transit transmission spectra of self-consistent N2-CO2-H2O atmospheres, and analytic N2-H2 and N2-H2-CO2 atmospheres. We show that (N2)2 absorption in the wings of the 4.3 μm CO2 band is strongly dependent on N2 partial pressures above 0.5 bar and can significantly widen this band in thick N2 atmospheres. The (N2)2 transit transmission signal is up to 10 ppm for an Earth-size planet with an N2-dominated atmosphere orbiting within the habitable zone of an M5V star and could be substantially larger for planets with significant H2 mixing ratios.

Schwieterman, Edward W.↗

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane-parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of nonthermal electrons and X-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick H I, Ca II and Mg II transitions using the non-local thermodynamic equilibrium (non-LTE) prescription in the program MULTI. Hydrogen ionization rates from X-ray photoionization and collisional ionization by nonthermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in the pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring. We have computed the chromospheric radiation that results from a range of coronal heating rates, with particular emphasis on the widely observed diagnostic H(alpha). Our conclusion is that the H(alpha) fluxes and profiles actually observed in flares can only be produced under conditions of a low-pressure corona with strong beam heating. Therefore we suggest that H(alpha) in flares is produced primarily at the footprints of newly heated loops where significant evaporation has not yet occurred. As a single loop evolves in time, no matter how strong the heating rate may become, the H(alpha) flux will diminish as the corona becomes denser and hence more effective at stopping the beam. This prediction leads to several observable consequences regarding the spatial and temporal signatures of the X-ray and H(alpha) radiation during flares.

Hawley, Suzanne L.↗

Multiple outflows in the bipolar planetary nebula M1-16: A molecular line study

Extensive observations of the molecular gas in the young, compact planetary nebula M1-16 have been made, using the Swedish-ESO-Submillimeter Telescope. A map of the CO J = 2-1 emission shows that the molecular envelope contains both a slow and a fast outflow with expansion velocities of 19 km/s and greater than 34 km/s, respectively. The slow outflow is mildly elliptical, while the fast molecular outflow is bipolar. This fast outflow is roughly aligned with the very fast outflows recently found in the optical, while the long axis of the slow elliptical outflow is roughly orthogonal to the optical outflow axis. The kinematic timescales for the CO fast outflow and the optical very fast outflow agree closely, supporting the view that the former represents material in the slow outflow accelerated by the very fast outflow. The kinematic signature of a disk expanding with about 15.5 km/s can also be seen in the CO J = 2-1 data. The mass-loss rate (a) for the slow outflow is greater than or equal to 2.8 x 10(exp -5) solar mass/yr and possibly as large as 9 x 10(exp -5) solar mass/yr, (b) for the fast outflow is greater than or equal to 5 x 10(exp -6) solar mass/yr, and (c) for the very fast optically visible outflow is approximately equal 5 x 10(exp -7) solar mass/yr. The disk mass is approximately equal 6 x 10(exp -3) solar mass. Grain photoelectric heating results in temperatures of 20-70 K in molecular gas of the slow outflow. The (13)C/(12)C abundance ratio in M1-16 is found to be 0.33, quite possibly the highest found for any evolved object. Upper limits for the (18)O/(16)O and (17)O/(16)O ratios were found to be consistent with the values found in AGB stars. A search for other molecular species in M1-16 resulted in the detection of the high-excitation species HCN, CN, (13)CN, HCO(+), and H(13)CO(+) and possibly N2H(+). Both the HCO(+)/HCN and CN/HCN line-intensity ratios are enhanced, the former by a very large factor, over the values found in the envelopes of AGB stars, probably as a result of enhancement of the CN and HCO(+) abundances due to photochemistry induced by the stellar UV. The CS J = 2-1, SiO J = 2-1 (v = 0), and SiS J = 6-5 lines were not detected to low levels. For the high-excitation molecules, adequate collisional excitation of rotational levels and survival against photodissociation by the UV radiation requires significant clumping of the molecular gas into clumps with H2 densities approximately 10(exp 5)/cu cm. The IRAS fluxes of M1-16, assuming negligible contribution from line emission, imply the presence of about (1.7-0.4) x 10(exp -3) solar mass of cool dust (temperature around 50 K) and a smaller quantity, (2.7-3.1) x 10(exp -6) solar mass, of warmer dust (temperature around 125 K) for a power-law emissivity index p = 1-2. The evolutionary nature of M1-16 cannot be explained by existing single-star models of post-AGB evolution. The very high (13)C/(12)C abundance ratio in M1-16 suggests a possible evolutionary connection between M1-16 and the rare class of J-type silicate-carbon stars which also have high (13)C/(12)C ratios and are thought to be binary systems with accretion disks.

Sahai, Raghvendra↗

The heating of interstellar clouds by vibrationally excited molecular hydrogen

The possibility that vibrationally excited H2 may be collisionally de-excited, so providing a heating mechanism for interstellar clouds which operates by coupling the stellar radiation to the gas, is discussed. The majority of excitations in the Lyman and Werner bands of H2 return the molecules to the ground electronic state in a vibrationally excited level, the most favored level being 7. The heating rate obtained in this way is compared with other mechanisms which have been postulated, and the results of calculations of temperature as a function of depth into clouds of different densities are presented. It appears that this mechanism is a significant one, which should be taken into account in detailed models of dense clouds.

Stecher, T. P.↗

The High-Excitation Planetary Nebula NCG 246: Optical and Near-Ultraviolet Observations and Two-dimensional Numerical Models

We have imaged the planetary nebula (PN) NGC 246 in the near-ultraviolet wavelengths [Ne v] 342.6 nm, the Bowen fluorescence line of 0 111 at 344.4 nm, and a nearby line-free region centered on 338.6 nm, as well as H(alpha), [O III] 500.7 nm, and [S II] 673.0 and 671.5 nm. Imaging in the 344.4 nm line is necessary to deconvolve contamination of the [Ne v] images by O III 342.9 nm. The emission from the shell and inner parts of the nebula is detected in [Ne v]. The radial profiles of the [Ne v] brightness decrease with radius from the exciting star, indicating that the bulk of the emission from this ion is due to the hard UV stellar radiation field, with a (probably) marginal contribution from collisional ionization in a shock between the PN shell and the interstellar medium (ISM). In contrast, the radial profiles of the emission in H(alpha), [0 III] 500.7 nm, and [S II] are flatter and peak at the location of the shell. The emission of [S II] probably traces the interaction of the PN with the ambient ISM. We also present two-dimensional numerical simulations for this PN-ISM interaction. The simulations consider the stellar motion with respect to the ambient ISM, with a velocity of 85 km/s , and include the time evolution of the wind parameters and UV radiation field from the progenitor star.

Szentgyorgyi, Andrew↗

Sensitivity of Arctic Clouds to Ice Microphysical Processes in the NorESM2 Climate Model

Abstract Ice formation remains one of the most poorly represented microphysical processes in climate models. While primary ice production (PIP) parameterizations are known to have a large influence on the modeled cloud properties, the representation of secondary ice production (SIP) is incomplete and its corresponding impact is therefore largely unquantified. Furthermore, ice aggregation is another important process for the total cloud ice budget, which also remains largely unconstrained. In this study, we examine the impact of PIP, SIP, and ice aggregation on Arctic clouds, using the Norwegian Earth System Model, version 2 (NorESM2). Simulations with both prognostic and diagnostic PIP show that heterogeneous freezing alone cannot reproduce the observed cloud ice content. The implementation of missing SIP mechanisms (collisional breakup, drop shattering, and sublimation breakup) in NorESM2 improves the modeled ice properties, while improvements in liquid content occur only in simulations with prognostic PIP. However, results are sensitive to the description of collisional breakup. This mechanism, which dominates SIP in the examined conditions, is very sensitive to the treatment of the sublimation correction factor, a poorly constrained parameter that is included in the utilized parameterization. Finally, variations in ice aggregation treatment can also significantly impact cloud properties, mainly through their impact on collisional breakup efficiency. Overall, enhancement in ice production through the addition of SIP mechanisms and the reduction in ice aggregation (in line with radar observations of shallow Arctic clouds) result in enhanced cloud cover and decreased TOA radiation biases, compared to satellite measurements, especially during the cold months. Significance Statement Arctic clouds remain a large source of uncertainty in projections of the future climate due to the poor representation of the microphysical processes that govern their life cycle. Ice formation is among the least understood processes. While it is widely recognized that better constraints on primary ice production (PIP) are needed to improve existing parameterizations, we show that secondary ice production (SIP) and ice aggregation can have also a significant impact on ice number concentrations. Constraining ice formation through the addition of missing SIP mechanisms and reducing ice aggregation can improve the representation of the cloud macrophysical properties and enhance total cloud cover in the Arctic region, which in turn contributes to decreased TOA radiation biases in the cold months.

Meteorology & Atmospheric Sciences↗

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗