Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “OPACITY”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Effective Tree Scattering at L-Band

For routine microwave Soil Moisture (SM) retrieval through vegetation, the tau-omega [1] model [zero-order Radiative Transfer (RT) solution] is attractive due to its simplicity and eases of inversion and implementation. It is the model used in baseline retrieval algorithms for several planned microwave space missions, such as ESA's Soil Moisture Ocean Salinity (SMOS) mission (launched November 2009) and NASA's Soil Moisture Active Passive (SMAP) mission (to be launched 2014/2015) [2 and 3]. These approaches are adapted for vegetated landscapes with effective vegetation parameters tau and omega by fitting experimental data or simulation outputs of a multiple scattering model [4-7]. The model has been validated over grasslands, agricultural crops, and generally light to moderate vegetation. As the density of vegetation increases, sensitivity to the underlying SM begins to degrade significantly and errors in the retrieved SM increase accordingly. The zero-order model also loses its validity when dense vegetation (i.e. forest, mature corn, etc.) includes scatterers, such as branches and trunks (or stalks in the case of corn), which are large with respect to the wavelength. The tau-omega model (when applied over moderately to densely vegetated landscapes) will need modification (in terms of form or effective parameterization) to enable accurate characterization of vegetation parameters with respect to specific tree types, anisotropic canopy structure, presence of leaves and/or understory. More scattering terms (at least up to first-order at L-band) should be included in the RT solutions for forest canopies [8]. Although not really suitable to forests, a zero-order tau-omega model might be applied to such vegetation canopies with large scatterers, but that equivalent or effective parameters would have to be used [4]. This requires that the effective values (vegetation opacity and single scattering albedo) need to be evaluated (compared) with theoretical definitions of these parameters. In a recent study [9], effective vegetation opacity of coniferous trees was compared with two independent estimates of the same parameter. First, a zero-order RT model was fitted to multiangular microwave emissivity data in a least-square sense to provide effective vegetation optical depth as done in spaceborne retrieval algorithms. Second, a ratio between radar backscatter measurements with a corner reflector under trees and in an open area was calculated to obtain measured tree propagation characteristics. Finally, the theoretical propagation constant was determined by forward scattering theorem using detailed measurements of size/angle distributions and dielectric constants of the tree constituents (trunk, branches, and needles). Results indicated that the effective attenuation values are smaller than but of similar magnitude to both the theoretical and measured values. This study will complement the previous work [9] and will focus on characterization of effective scattering albedo by assuming that effective vegetation opacity is same as theoretical opacity. The resultant effective albedo will not be the albedo of single forest canopy element anymore, but it becomes a global parameter, which depends on all the processes taking place within the canopy including multiple scattering as described.

Kurum, Mehmet↗
NASA NTRS2011

Interannual Variability of Dust and Ice in the Mars Atmosphere: Comparison of MRO Mars Climate Sounder Retrievals with MGS-TES Limb Sounding Retrievals

Dust and ice play important roles in Martian atmospheric dynamics on all time scales. Dust loading in particular exerts an important control on atmospheric temperatures and thereby on the strength of the atmospheric circulation in any given year. We present the first comparisons of MGS-TES aerosol opacity profiles with MRO-MCS aerosol opacity profiles. While the differences in vertical resolution are significant (a factor of 2), we find good agreement at particular seasons between nightside zonal average dust opacity profiles from the two instruments. Derived water ice opacities are likewise similar but show greater variability.

Mars Climate Sounder (MCS)↗
NASA NTRS2012

Titan's Tropopause Temperatures from CIRS: Implications for Stratospheric Methane Cloud Formation

Analysis of Cassini Composite Infrared Spectrometer (CIRS) far-IR spectra enable the construction of Titan's temperature profile in the altitude region containing the tropopause. Whereas the methane V4 band at 1306/cm (7.7 microns) is the primary opacity source for deducing thermal structure between 100 km and 500 km, N2-N2 collision-induced absorption between 70 and 140/cm (143 microns and 71 microns) is utilized to determine temperatures at Titan's tropopause. Additional opacity due to aerosol and nitrile ices must also be taken into account in this part of the far-IR spectral region. The spectral characteristics of these particulate opacities have been deduced from CIRS limb data at 58degS, 15degS, 15degN, and 85degN. Empirically, the spectral shapes of these opacities appear to be independent of both latitude and altitude below 300 km (Anderson and Samuelson, 2011, Icarus 212, 762-778), justifying the extension of these spectral properties to all latitudes. We find that Titan's tropopause temperature is cooler than the HAS! value of 70.5K by approx. 6K. This leads to the possibility that subsidence at high northern latitudes can cause methane condensation in the winter polar stratosphere. A search for methane clouds in this region is in progress.

Anderson, C. M.↗
DOE OSTI2018

Comparison of Shallow and Deep Learning Methods on Classifying the Regional Pattern of Diffuse Lung Disease

This study aimed to compare shallow and deep learning of classifying the patterns of interstitial lung diseases (ILDs). Using high-resolution computed tomography images, two experienced radiologists marked 1200 regions of interest (ROIs), in which 600 ROIs were each acquired using a GE or Siemens scanner and each group of 600 ROIs consisted of 100 ROIs for subregions that included normal and five regional pulmonary disease patterns (ground-glass opacity, consolidation, reticular opacity, emphysema, and honeycombing). We employed the convolution neural network (CNN) with six learnable layers that consisted of four convolution layers and two fully connected layers. The classification results were compared with the results classified by a shallow learning of a support vector machine (SVM). The CNN classifier showed significantly better performance for accuracy compared with that of the SVM classifier by 6–9%. As the convolution layer increases, the classification accuracy of the CNN showed better performance from 81.27 to 95.12%. Especially in the cases showing pathological ambiguity such as between normal and emphysema cases or between honeycombing and reticular opacity cases, the increment of the convolution layer greatly drops the misclassification rate between each case. Conclusively, the CNN classifier showed significantly greater accuracy than the SVM classifier, and the results implied structural characteristics that are inherent to the specific ILD patterns.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Mean-force scattering potential for calculating optical properties of dense plasmas

Here, we assess the relative importance of ionic structure on the opacity of dense plasmas by using the potential of mean force as a scattering potential within the Kubo-Greenwood formalism. We compare results from the potential of mean force, which includes realistic ionic structure, to results using an average atom potential, which includes only a crude treatment of ionic structure. Comparisons with less approximate but more expensive DFT-MD simulations for aluminum plasma reveal that the mean force generally improves agreement for DC conductivity. We also see improvement when applying the mean force to free-free transitions, whereas for bound-bound and bound-free transitions the mean force leads to poorer agreement on transition energies. Further, we assess the impact of accounting for correlations within the plasma at the temperature and density conditions relevant to iron opacity measurements at Sandia’s Z machine facility and find that these correlations do not account for the discrepancy between the measurements and leading opacity calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The STAG code: A fully relativistic super transition array calculation using Green’s functions

Calculating opacities for a wide range of plasma conditions (i.e. temperature, density, element) requires detailed knowledge of the plasma configuration space and electronic structure. For plasmas composed of heavier elements, relativistic effects are important in both the electronic structure and the details of opacity spectra. We extend our previously described superconfiguration and super transition array capabilities (Gill et al., 2023) to include a fully relativistic formalism. Further, the use of hybrid bound-continuum supershells in our superconfigurations demonstrates the importance of a consistent treatment of bound and continuum electrons in dense plasma opacities, and we expand the discussion of these consequences to include issues associated with equation of state and electron correlations between bound and continuum electrons.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Planck 2013 results. XI. All-sky model of thermal dust emission

This work introduces an all-sky model of dust emission from the Planck 353, 545, and 857 GHz, and IRAS 100 μm data. Using a modified blackbody fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a good representation of the IRAS and Planck data at 5' between 353 and 3000 GHz (850 and 100 μm). It shows variations of the order of 30% compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, down to an angular resolution of 5', providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions where the dust temperature varies strongly at small scales in response to dust evolution, extinction, and/or local production of heating photons. An increase of the dust opacity at 353 GHz, τ 353 / N H , from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T obs , that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at H I column densities lower than 10 20 cm -2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anticorrelation between τ 353 / N H and T obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant, confirming one of the Planck Early Results obtained on selected fields. This effect is compatible with the view that, in the diffuse ISM, T obs responds to spatial variations of the dust opacity, due to variations of dust properties, in addition to (small) variations of the radiation field strength. The indication is that in the diffuse high-latitude ISM τ 353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of τ 353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B -V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Investigation of mechanisms for He-I emission radial profile broadening in a weakly ionized cylindrical helium plasma with recombining edge

Spatially resolved spectroscopic measurements of He-I line emission are used to study the causes of emission profile broadening radially across the cylinder of a weakly ionized helium plasma. The plasma consists of an ionizing core (r < 2 cm) surrounded by a recombining edge (r > 2 cm) plasma. The brightness profiles of low-n EUV He-I resonance lines are shown to be strongly radially broadened due to opacity. The brightness profiles of high-n visible lines are also found to be strongly radially broadened, but dominantly due to edge recombination. Visible low-n lines are less strongly radially broadened, apparently by a combination of both recombination and EUV opacity. The low-n visible He-I line ratio method with central opacity correction is found to calculate central electron density and temperature well, with poor agreement at the edge, as expected for recombining plasma. In the recombining edge, high-n Boltzmann analysis is found to accurately measure the cold (T{sub e} < 0.2 eV) edge temperature. Near the core, however, high-n Boltzmann analysis can be complicated by electron-impact excitation, giving incorrect (T{sub e}≈ 0.1 eV) apparent temperatures. Probe measurements were not able to capture the cold edge temperature accurately, probably due to large potential fluctuations, even when using fast triple probe measurements. Fast spectroscopic measurements show that this discrepancy is not explained by recombining plasma alternating with ionizing plasma in the edge region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Vector correlation between the alignment of reactant N{sub 2} (A {sup 3}Σ{sub u}{sup +}) and the alignment of product NO (A {sup 2}Σ{sup +}) rotation in the energy transfer reaction of aligned N{sub 2} (A {sup 3}Σ{sub u}{sup +}) + NO (X {sup 2}Π) → NO (A {sup 2}Σ{sup +}) + N{sub 2} (X {sup 1}Σ{sub g}{sup +})

The vector correlation between the alignment of reactant N{sub 2} (A {sup 3}Σ{sub u}{sup +}) and the alignment of product NO (A {sup 2}Σ{sup +}) rotation has been studied in the energy transfer reaction of aligned N{sub 2} (A {sup 3}Σ{sub u}{sup +}) + NO (X {sup 2}Π) → NO (A {sup 2}Σ{sup +}) + N{sub 2} (X {sup 1}Σ{sub g}{sup +}) under the crossed beam condition at a collision energy of ∼0.07 eV. NO (A {sup 2}Σ{sup +}) emission in the two linear polarization directions (i.e., parallel and perpendicular with respect to the relative velocity vector v{sub R}) has been measured as a function of the alignment of N{sub 2} (A {sup 3}Σ{sub u}{sup +}) along its molecular axis in the collision frame. The degree of polarization of NO (A {sup 2}Σ{sup +}) emission is found to depend on the alignment angle (θ{sub v{sub R}}) of N{sub 2} (A {sup 3}Σ{sub u}{sup +}) in the collision frame. The shape of the steric opacity function at the two polarization conditions turns out to be extremely different from each other: The steric opacity function at the parallel polarization condition is more favorable for the oblique configuration of N{sub 2} (A {sup 3}Σ{sub u}{sup +}) at an alignment angle of θ{sub v{sub R}} ∼ 45° as compared with that at the perpendicular polarization condition. The alignment of N{sub 2} (A {sup 3}Σ{sub u}{sup +}) is found to give a significant effect on the alignment of NO (A {sup 2}Σ{sup +}) rotation in the collision frame: The N{sub 2} (A {sup 3}Σ{sub u}{sup +}) configuration at an oblique alignment angle θ{sub v{sub R}} ∼ 45° leads to a parallel alignment of NO (A {sup 2}Σ{sup +}) rotation (J-vector) with respect to v{sub R}, while the axial and sideways configurations of N{sub 2} (A {sup 3}Σ{sub u}{sup +}) lead to a perpendicular alignment of NO (A {sup 2}Σ{sup +}) rotation with respect to v{sub R}. These stereocorrelated alignments of the product rotation have a good correlation with the stereocorrelated reactivity observed in the multi-dimensional steric opacity function [H. Ohoyama and S. Maruyama, J. Chem. Phys. 137, 064311 (2012)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

INFRARED TRANSMISSION SPECTROSCOPY OF THE EXOPLANETS HD 209458b AND XO-1b USING THE WIDE FIELD CAMERA-3 ON THE HUBBLE SPACE TELESCOPE

Exoplanetary transmission spectroscopy in the near-infrared using the Hubble Space Telescope (HST) NICMOS is currently ambiguous because different observational groups claim different results from the same data, depending on their analysis methodologies. Spatial scanning with HST/WFC3 provides an opportunity to resolve this ambiguity. We here report WFC3 spectroscopy of the giant planets HD 209458b and XO-1b in transit, using spatial scanning mode for maximum photon-collecting efficiency. We introduce an analysis technique that derives the exoplanetary transmission spectrum without the necessity of explicitly decorrelating instrumental effects, and achieves nearly photon-limited precision even at the high flux levels collected in spatial scan mode. Our errors are within 6% (XO-1) and 26% (HD 209458b) of the photon-limit at a resolving power of {lambda}/{delta}{lambda} {approx} 70, and are better than 0.01% per spectral channel. Both planets exhibit water absorption of approximately 200 ppm at the water peak near 1.38 {mu}m. Our result for XO-1b contradicts the much larger absorption derived from NICMOS spectroscopy. The weak water absorption we measure for HD 209458b is reminiscent of the weakness of sodium absorption in the first transmission spectroscopy of an exoplanet atmosphere by Charbonneau et al. Model atmospheres having uniformly distributed extra opacity of 0.012 cm{sup 2} g{sup -1} account approximately for both our water measurement and the sodium absorption. Our results for HD 209458b support the picture advocated by Pont et al. in which weak molecular absorptions are superposed on a transmission spectrum that is dominated by continuous opacity due to haze and/or dust. However, the extra opacity needed for HD 209458b is grayer than for HD 189733b, with a weaker Rayleigh component.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The influence of differential irradiation and circulation on the thermal evolution of gas giant planets. I. Upper limits from radiative equilibrium

As a planet ages, it cools and its radius shrinks at a rate set by the efficiency with which heat is transported from the interior out to space. The bottleneck for this transport is at the boundary between the convective interior and the radiative atmosphere; the opacity there sets the global cooling rate. Models of planetary evolution are often one dimensional (1D), such that the radiative-convective boundary (RCB) is defined by a single temperature, pressure, and opacity. In reality the spatially inhomogeneous stellar heating pattern and circulation in the atmosphere could deform the RCB, allowing heat from the interior to escape more efficiently through regions with lower opacity. We present an analysis of the degree to which the RCB could be deformed and the resultant change in the evolutionary cooling rate. In this initial work we calculate the upper limit for this effect by comparing an atmospheric structure in local radiative equilibrium to its 1D equivalent. We find that the cooling through an uneven RCB could be enhanced over cooling through a uniform RCB by as much as 10%-50%. We also show that the deformation of the RCB (and the enhancement of the cooling rate) increases with a greater incident stellar flux or a lower inner entropy. Our results indicate that this mechanism could significantly change a planet's thermal evolution, causing it to cool and shrink more quickly than would otherwise be expected. This may exacerbate the well-known difficulty in explaining the very large radii observed for some hot Jupiters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The atmospheric circulation of the super Earth GJ 1214b: Dependence on composition and metallicity

We present three-dimensional atmospheric circulation models of GJ 1214b, a 2.7 Earth-radius, 6.5 Earth-mass super Earth detected by the MEarth survey. Here we explore the planet's circulation as a function of atmospheric metallicity and atmospheric composition, modeling atmospheres with a low mean molecular weight (MMW; i.e., H{sub 2}-dominated) and a high MMW (i.e., water- and CO{sub 2}-dominated). We find that atmospheres with a low MMW have strong day-night temperature variations at pressures above the infrared photosphere that lead to equatorial superrotation. For these atmospheres, the enhancement of atmospheric opacities with increasing metallicity lead to shallower atmospheric heating, larger day-night temperature variations, and hence stronger superrotation. In comparison, atmospheres with a high MMW have larger day-night and equator-to-pole temperature variations than low MMW atmospheres, but differences in opacity structure and energy budget lead to differences in jet structure. The circulation of a water-dominated atmosphere is dominated by equatorial superrotation, while the circulation of a CO{sub 2}-dominated atmosphere is instead dominated by high-latitude jets. By comparing emergent flux spectra and light curves for 50× solar and water-dominated compositions, we show that observations in emission can break the degeneracy in determining the atmospheric composition of GJ 1214b. The variation in opacity with wavelength for the water-dominated atmosphere leads to large phase variations within water bands and small phase variations outside of water bands. The 50× solar atmosphere, however, yields small variations within water bands and large phase variations at other characteristic wavelengths. These observations would be much less sensitive to clouds, condensates, and hazes than transit observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

On the minimum core mass for giant planet formation at wide separations

In the core accretion hypothesis, giant planets form by gas accretion onto solid protoplanetary cores. The minimum (or critical) core mass to form a gas giant is typically quoted as 10 M {sub ⊕}. The actual value depends on several factors: the location in the protoplanetary disk, atmospheric opacity, and the accretion rate of solids. Motivated by ongoing direct imaging searches for giant planets, this study investigates core mass requirements in the outer disk. To determine the fastest allowed rates of gas accretion, we consider solid cores that no longer accrete planetesimals, as this would heat the gaseous envelope. Our spherical, two-layer atmospheric cooling model includes an inner convective region and an outer radiative zone that matches onto the disk. We determine the minimum core mass for a giant planet to form within a typical disk lifetime of 3 Myr. The minimum core mass declines with disk radius, from ∼8.5 M {sub ⊕} at 5 AU to ∼3.5 M {sub ⊕} at 100 AU, with standard interstellar grain opacities. Lower temperatures in the outer disk explain this trend, while variations in disk density are less influential. At all distances, a lower dust opacity or higher mean molecular weight reduces the critical core mass. Our non-self-gravitating, analytic cooling model reveals that self-gravity significantly affects early atmospheric evolution, starting when the atmosphere is only ∼10% as massive as the core.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

How gravitational lensing helps γ-ray photons avoid γ – γ absorption

We investigate potential γ – γ absorption of γ-ray emission from blazars arising from inhomogeneities along the line of sight, beyond the diffuse Extragalactic Background Light (EBL). As plausible sources of excess γ – γ opacity, we consider (1) foreground galaxies, including cases in which this configuration leads to strong gravitational lensing, (2) individual stars within these foreground galaxies, and (3) individual stars within our own galaxy, which may act as lenses for microlensing events. We found that intervening galaxies close to the line of sight are unlikely to lead to significant excess γ – γ absorption. This opens up the prospect of detecting lensed gamma-ray blazars at energies above 10 GeV with their gamma-ray spectra effectively only affected by the EBL. The most luminous stars located either in intervening galaxies or in our galaxy provide an environment in which these gamma-rays could, in principle, be significantly absorbed. However, despite a large microlensing probability due to stars located in intervening galaxies, γ-rays avoid absorption by being deflected by the gravitational potentials of such intervening stars to projected distances ({sup i}mpact parameters{sup )} where the resulting γ – γ opacities are negligible. Thus, neither of the intervening excess photon fields considered here, provide a substantial source of excess γ – γ opacity beyond the EBL, even in the case of very close alignments between the background blazar and a foreground star or galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

REVERSE SHOCK EMISSION AND IONIZATION BREAKOUT POWERED BY POST-MERGER MILLISECOND MAGNETARS

There is accumulating evidence that at least a fraction of binary neutron star (BNS) mergers result in rapidly spinning magnetars, with subrelativistic neutron-rich ejecta as massive as a small fraction of a solar mass. The ejecta could be heated continuously by the Poynting flux emanated from the central magnetars. Such a Poynting flux could become lepton dominated so that a reverse shock develops. It was demonstrated that such a picture is capable of accounting for the optical transient PTF11agg. In this paper we investigate the X-ray and ultraviolet (UV) radiation, as well as the optical and radio radiation studied by Wang and Dai. UV emission is particularly important because it has the right energy to ionize the hot ejecta at times t ≲ 600 s. It is thought that the ejecta of BNS mergers are a remarkably pure sample of r-process material, about which our understanding is still incomplete. In this paper we evaluate the possibility of observationally determining the bound-bound and bound-free opacities of the r-process material by timing the X-ray, UV, and optical radiation. It is found that these timings depend on the opacities weakly, and therefore only loose constraints on the opacities can be obtained.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SPECTROSCOPIC INFRARED EXTINCTION MAPPING AS A PROBE OF GRAIN GROWTH IN IRDCs

We present spectroscopic tests of MIR to FIR extinction laws in IRDC G028.36+00.07, a potential site of massive star and star cluster formation. Lim and Tan developed methods of FIR extinction mapping of this source using Spitzer-MIPS 24 μm and Herschel-PACS 70 μm images, and by comparing to MIR Spitzer-IRAC 3–8 μm extinction maps, found tentative evidence for grain growth in the highest mass surface density regions. Here we present results of spectroscopic infrared extinction mapping using Spitzer-IRS (14–38 μm) data of the same Infrared dark cloud (IRDC). These methods allow us to first measure the SED of the diffuse Galactic interstellar medium that is in the foreground of the IRDC. We then carry out our primary investigation of measuring the MIR to FIR opacity law and searching for potential variations as a function of mass surface density within the IRDC. We find relatively flat, featureless MIR–FIR opacity laws that lack the ∼12 and ∼35 μm features associated with the thick water ice mantle models of Ossenkopf and Henning. Their thin ice mantle models and the coagulating aggregate dust models of Ormel et al. are a generally better match to the observed opacity laws. We also find evidence for generally flatter MIR to FIR extinction laws as mass surface density increases, strengthening the evidence for grain and ice mantle growth in higher density regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

AN ATMOSPHERIC STRUCTURE EQUATION FOR GRAIN GROWTH

We present a method to include the evolution of the grain size and grain opacity κ{sub gr} in the equations describing the structure of protoplanetary atmospheres. The key assumption of this method is that a single grain size dominates the grain size distribution at any height r. In addition to following grain growth, the method accounts for mass deposition by planetesimals and grain porosity. We illustrate this method by computation of a simplified atmosphere structure model. In agreement with previous works, grain coagulation is seen to be very efficient. The opacity drops to values much below the often-used ''interstellar medium opacities'' (∼1 cm{sup 2} g{sup –1}) and the atmosphere structure profiles for temperature and density resemble that of the grain-free case. Deposition of planetesimals in the radiative part of the atmosphere hardly influences this outcome as the added surface is quickly coagulated away. We observe a modest dependence on the internal structure (porosity), but show that filling factors cannot become too large because of compression by gas drag.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

PRINCIPAL COMPONENT ANALYSIS STUDIES OF TURBULENCE IN OPTICALLY THICK GAS

In this work we investigate the sensitivity of principal component analysis (PCA) to the velocity power spectrum in high-opacity regimes of the interstellar medium (ISM). For our analysis we use synthetic position–position–velocity (PPV) cubes of fractional Brownian motion and magnetohydrodynamics (MHD) simulations, post-processed to include radiative transfer effects from CO. We find that PCA analysis is very different from the tools based on the traditional power spectrum of PPV data cubes. Our major finding is that PCA is also sensitive to the phase information of PPV cubes and this allows PCA to detect the changes of the underlying velocity and density spectra at high opacities, where the spectral analysis of the maps provides the universal −3 spectrum in accordance with the predictions of the Lazarian and Pogosyan theory. This makes PCA a potentially valuable tool for studies of turbulence at high opacities, provided that proper gauging of the PCA index is made. However, we found the latter to not be easy, as the PCA results change in an irregular way for data with high sonic Mach numbers. This is in contrast to synthetic Brownian noise data used for velocity and density fields that show monotonic PCA behavior. We attribute this difference to the PCA's sensitivity to Fourier phase information.

79 ASTRONOMY AND ASTROPHYSICS↗