Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STAR MODELS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

New constraints on the cooling of the central compact object in CAS A

To examine the previously claimed fast cooling of the Central Compact Object (CCO) in the Cas A supernova remnant (SNR), we analyzed two Chandra observations of this CCO, taken in a setup minimizing instrumental spectral distortions. We fit the two CCO X-ray spectra from 2006 and 2012 with hydrogen and carbon neutron star atmosphere models. The temperature and flux changes in the 5.5 yr between the two epochs depend on the adopted constraints on the fitting parameters and the uncertainties of the effective area calibrations. If we allow a change of the equivalent emitting region size, R {sub Em}, the effective temperature remains essentially the same. If R {sub Em} is held constant, the best-fit temperature change is negative, but its statistical significance ranges from 0.8σ to 2.5σ, depending on the model. If we assume that the optical depth of the ACIS filter contaminant in 2012 was ±10% different from its default calibration value, the significance of the temperature drop becomes 0.8σ-3.1σ, for the carbon atmospheres with constant R {sub Em}. Thus, we do not see a statistically significant temperature drop in our data, but the involved uncertainties are too large to firmly exclude the previously reported fast cooling. Our analysis indicate a decrease of 4%-6% (1.9σ-2.9σ significance) for the absorbed flux in the energy range 0.6-6 keV between 2006 and 2012, most prominent in the ≈1.4-1.8 keV energy range. It could be caused by unaccounted changes of the detector response or contributions from unresolved SNR material along the line of sight to the CCO.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

RCSED—A Value-added Reference Catalog of Spectral Energy Distributions of 800,299 Galaxies in 11 Ultraviolet, Optical, and Near-infrared Bands: Morphologies, Colors, Ionized Gas, and Stellar Population Properties

We present RCSED, the value-added Reference Catalog of Spectral Energy Distributions of galaxies, which contains homogenized spectrophotometric data for 800,299 low- and intermediate-redshift galaxies (0.007<z<0.6) selected from the Sloan Digital Sky Survey spectroscopic sample. Accessible from the Virtual Observatory (VO) and complemented with detailed information on galaxy properties obtained with state-of-the-art data analysis, RCSED enables direct studies of galaxy formation and evolution over the last 5 Gyr. We provide tabulated color transformations for galaxies of different morphologies and luminosities, and analytic expressions for the red sequence shape in different colors. RCSED comprises integrated k-corrected photometry in up to 11 ultraviolet, optical, and near-infrared bands published by the GALEX, SDSS, and UKIDSS wide-field imaging surveys; results of the stellar population fitting of SDSS spectra including best-fitting templates, velocity dispersions, parameterized star formation histories, and stellar metallicities computed for instantaneous starburst and exponentially declining star formation models; parametric and non-parametric emission line fluxes and profiles; and gas phase metallicities. We link RCSED to the Galaxy Zoo morphological classification and galaxy bulge+disk decomposition results of Simard et al. We construct the color–magnitude, Faber–Jackson, and mass–metallicity relations; compare them with the literature; and discuss systematic errors of the galaxy properties presented in our catalog. RCSED is accessible from the project web site and via VO simple spectrum access and table access services using VO-compliant applications. We describe several examples of SQL queries to the database. Finally, we briefly discuss existing and future scientific applications of RCSED and prospective catalog extensions to higher redshifts and different wavelengths.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A COUPLED 2 × 2D BABCOCK–LEIGHTON SOLAR DYNAMO MODEL. I. SURFACE MAGNETIC FLUX EVOLUTION

The need for reliable predictions of the solar activity cycle motivates the development of dynamo models incorporating a representation of surface processes sufficiently detailed to allow assimilation of magnetographic data. In this series of papers we present one such dynamo model, and document its behavior and properties. This first paper focuses on one of the model’s key components, namely surface magnetic flux evolution. Using a genetic algorithm, we obtain best-fit parameters of the transport model by least-squares minimization of the differences between the associated synthetic synoptic magnetogram and real magnetographic data for activity cycle 21. Our fitting procedure also returns Monte Carlo-like error estimates. We show that the range of acceptable surface meridional flow profiles is in good agreement with Doppler measurements, even though the latter are not used in the fitting process. Using a synthetic database of bipolar magnetic region (BMR) emergences reproducing the statistical properties of observed emergences, we also ascertain the sensitivity of global cycle properties, such as the strength of the dipole moment and timing of polarity reversal, to distinct realizations of BMR emergence, and on this basis argue that this stochasticity represents a primary source of uncertainty for predicting solar cycle characteristics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

PROBING THE ROTATION OF CORE-COLLAPSE SUPERNOVA WITH A CONCURRENT ANALYSIS OF GRAVITATIONAL WAVES AND NEUTRINOS

The next time a core-collapse supernova (SN) explodes in our galaxy, various detectors will be ready and waiting to detect its emissions of gravitational waves (GWs) and neutrinos. Current numerical simulations have successfully introduced multi-dimensional effects to produce exploding SN models, but thus far the explosion mechanism is not well understood. In this paper, we focus on an investigation of progenitor core rotation via comparison of the start time of GW emission and that of the neutronization burst. The GW and neutrino detectors are assumed to be, respectively, the KAGRA detector and a co-located gadolinium-loaded water Cherenkov detector, either EGADS or GADZOOKS!. Our detection simulation studies show that for a nearby SN (0.2 kpc) we can confirm the lack of core rotation close to 100% of the time, and the presence of core rotation about 90% of the time. Using this approach there is also the potential to confirm rotation for considerably more distant Milky Way SN explosions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

GAMMA-RAY BURSTS ARE OBSERVED OFF-AXIS

We constrain the jet opening angle and, for the first time, the off-axis observer angle for gamma-ray bursts in the Swift-XRT catalog by using the ScaleFit package to fit afterglow light curves directly to hydrodynamic simulations. The ScaleFit model uses scaling relations in the hydrodynamic and radiation equations to compute synthetic light curves directly from a set of high-resolution two-dimensional relativistic blast wave simulations. The data sample consists of all Swift-XRT afterglows from 2005 to 2012 with sufficient coverage and a known redshift, 226 bursts in total. We find that the jet half-opening angle varies widely but is commonly less than 0.1 rad. The distribution of the electron spectral index is also broad, with a median at 2.30. We find the observer angle to have a median value of 0.57 of the jet opening angle over our sample, which has profound consequences for the predicted rate of observed jet breaks and affects the beaming-corrected total energies of gamma-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Testing SgrA{sup *} with the spectrum of its accretion structure

SgrA{sup *} is the supermassive black hole candidate at the center of the Galaxy and an ideal laboratory to test general relativity. Following previous work by other authors, we use the Polish doughnut model to describe an optically thin and constant angular momentum ion torus in hydrodynamical equilibrium and model the accretion structure around SgrA{sup *}. The radiation mechanisms are bremsstrahlung, synchrotron emission, and inverse Compton scattering. We compute the spectrum as seen by a distant observer in Kerr and non-Kerr spacetimes and we study how an accurate measurement can constrain possible deviations form the Kerr solution. As in the case of emission from a thin accretion disk, we find a substantial degeneracy between the determination of the spin and of possible deviations from the Kerr geometry, even when the parameters of the ion torus are fixed. This means that this technique cannot independently test the nature of SgrA{sup *} even in the presence of good data and with the systematics under control. However, it might do it in combination with other measurements (black hole shadow, radio pulsar, etc.)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Hydrodynamic Photoevaporation of Protoplanetary Disks with Consistent Thermochemistry

Photoevaporation is an important dispersal mechanism for protoplanetary disks. We conduct hydrodynamic simulations coupled with ray-tracing radiative transfer and consistent thermochemistry to study photoevaporative winds driven by ultraviolet and X-ray radiation from the host star. Most models have a three-layer structure: a cold midplane, warm intermediate layer, and hot wind, the last having typical speeds ∼40 km s{sup −1} and mass-loss rates ∼10{sup −9} M{sub ⊙} yr{sup −1} when driven primarily by ionizing UV radiation. Observable molecules, including CO, OH, and H{sub 2}O re-form in the intermediate layer and survive at relatively high wind temperatures due to reactions being out of equilibrium. Mass-loss rates are sensitive to the intensity of radiation in energy bands that interact directly with hydrogen. Comparison with previous works shows that mass-loss rates are also sensitive to the treatment of both the hydrodynamics and thermochemistry. Divergent results concerning the efficiency of X-ray photoevaporation are traced in part to differing assumptions about dust and other coolants.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A 5 × 10{sup 9}M{sub ⊙} BLACK HOLE IN NGC 1277 FROM ADAPTIVE OPTICS SPECTROSCOPY

The nearby lenticular galaxy NGC 1277 is thought to host one of the largest black holes known, however the black hole mass measurement is based on low spatial resolution spectroscopy. In this paper, we present Gemini Near-infrared Integral Field Spectrometer observations assisted by adaptive optics. We map out the galaxy's stellar kinematics within ∼440 pc of the nucleus with an angular resolution that allows us to probe well within the region where the potential from the black hole dominates. We find that the stellar velocity dispersion rises dramatically, reaching ∼550 km s{sup −1} at the center. Through orbit-based, stellar-dynamical models we obtain a black hole mass of (4.9 ± 1.6) × 10{sup 9} M{sub ⊙} (1σ uncertainties). Although the black hole mass measurement is smaller by a factor of ∼3 compared to previous claims based on large-scale kinematics, NGC 1277 does indeed contain one of the most massive black holes detected to date, and the black hole mass is an order of magnitude larger than expectations from the empirical relation between black hole mass and galaxy luminosity. Given the galaxy's similarities to the higher redshift (z ∼ 2) massive quiescent galaxies, NGC 1277 could be a relic, passively evolving since that period. A population of local analogs to the higher redshift quiescent galaxies that also contain over-massive black holes may suggest that black hole growth precedes that of the host galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

SOFIA/FORCAST IMAGING OF THE CIRCUMNUCLEAR RING AT THE GALACTIC CENTER

We present 19.7, 31.5, and 37.1 μm images of the inner 6 pc of the Galactic center of the Milky Way with a spatial resolution of 3.''2-4.''6 taken by the Faint Object Infrared Camera on the Stratospheric Observatory for Infrared Astronomy. The images reveal in detail the 'clumpy' structure of the circumnuclear ring (CNR)—the inner edge of the molecular torus orbiting the supermassive black hole at the Galactic center—and the prominent streamers of hot, ionized gas and dust within the CNR that compose the H II region Sgr A West. The CNR exhibits features of a classic H II region: the dust emission at 19.7 μm closely traces the ionized gas emission observed in the radio while the 31.5 and 37.1 μm emission traces the photo-dissociation region beyond the ionized gas. The 19.7/37.1 color temperature map reveals a radial temperature gradient across the CNR with temperatures ranging from 65 to 85 K, consistent with the prevailing paradigm in which the dust is centrally heated by the inner cluster of hot, young stars. We model the 37.1 μm intensity of the CNR as an inclined (θ{sub i} = 67°) ring with a thickness and radius of 0.34 pc and 1.4 pc, respectively, and find that it is consistent with the observed 37.1 μm map of the CNR. The 37.1 μm optical depth map also reveals the clumpy dust distribution of the CNR and implies a total gas mass of ∼610 M{sub ☉}. Dense (5-9 × 10{sup 4} cm{sup –3}) clumps with an FWHM of ∼0.15 pc exist along the inner edge of the CNR and shadow the material deeper into the ring. We find that the clumps are unlikely to be long-lived structures since they are not dense enough to be stable against tidal shear from the supermassive black hole and will be sheared out on a timescale of an orbital period (∼10{sup 5} yr)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE IMACS CLUSTER BUILDING SURVEY. IV. THE LOG-NORMAL STAR FORMATION HISTORY OF GALAXIES

We present here a simple model for the star formation history (SFH) of galaxies that is successful in describing both the star formation rate density (SFRD) over cosmic time, as well as the distribution of specific star formation rates (sSFRs) of galaxies at the current epoch, and the evolution of this quantity in galaxy populations to a redshift of z = 1. We show first that the cosmic SFRD is remarkably well described by a simple log-normal in time. We next postulate that this functional form for the ensemble is also a reasonable description for the SFHs of individual galaxies. Using the measured sSFRs for galaxies at z {approx} 0 from Paper III in this series, we then construct a realization of a universe populated by such galaxies in which the parameters of the log-normal SFH of each galaxy are adjusted to match the sSFRs at z {approx} 0 as well as fitting, in ensemble, the cosmic SFRD from z = 0 to z = 8. This model predicts, with striking fidelity, the distribution of sSFRs in mass-limited galaxy samples to z = 1; this match is not achieved by other models with a different functional form for the SFHs of individual galaxies, but with the same number of degrees of freedom, suggesting that the log-normal form is well matched to the likely actual histories of individual galaxies. We also impose the sSFR versus mass distributions at higher redshifts from Paper III as constraints on the model, and show that, as previously suggested, some galaxies in the field, particularly low mass galaxies, are quite young at intermediate redshifts. As emphasized in Paper III, starbursts are insufficient to explain the enhanced sSFRs in intermediate redshift galaxies; we show here that a model using only smoothly varying log-normal SFHs for galaxies, which allows for some fraction of the population to have peak star formation at late times, does however fully explain the observations. Finally, we show that this model, constrained in detail only at redshifts z < 1, also produces the main sequence of star-formation observed at 1.5 < z < 2.5, again suggesting that the log-normal SFHs are a close approximation to the actual histories of typical galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

DUST AND GAS IN THE MAGELLANIC CLOUDS FROM THE HERITAGE HERSCHEL KEY PROJECT. I. DUST PROPERTIES AND INSIGHTS INTO THE ORIGIN OF THE SUBMILLIMETER EXCESS EMISSION

The dust properties in the Large and Small Magellanic clouds (LMC/SMC) are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 μm. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power-law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models, we investigate the origin of the submillimeter excess, defined as the submillimeter emission above that expected from SMBB models fit to observations <200 μm. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 μm submillimeter excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our fitting results show that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submillimeter excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 ± 1.7) × 10{sup 5} and (8.3 ± 2.1) × 10{sup 4} M {sub ☉} for the LMC and SMC, respectively. We find significant correlations between the submillimeter excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

MERIDIONAL CIRCULATION DYNAMICS FROM 3D MAGNETOHYDRODYNAMIC GLOBAL SIMULATIONS OF SOLAR CONVECTION

The form of solar meridional circulation is a very important ingredient for mean field flux transport dynamo models. However, a shroud of mystery still surrounds this large-scale flow, given that its measurement using current helioseismic techniques is challenging. In this work, we use results from three-dimensional global simulations of solar convection to infer the dynamical behavior of the established meridional circulation. We make a direct comparison between the meridional circulation that arises in these simulations and the latest observations. Based on our results, we argue that there should be an equatorward flow at the base of the convection zone at mid-latitudes, below the current maximum depth helioseismic measures can probe (0.75 R{sub ⊙}). We also provide physical arguments to justify this behavior. The simulations indicate that the meridional circulation undergoes substantial changes in morphology as the magnetic cycle unfolds. We close by discussing the importance of these dynamical changes for current methods of observation which involve long averaging periods of helioseismic data. Also noteworthy is the fact that these topological changes indicate a rich interaction between magnetic fields and plasma flows, which challenges the ubiquitous kinematic approach used in the vast majority of mean field dynamo simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Alfvén wave solar model (AWSoM): Coronal heating

We present a new version of the Alfvén wave solar model, a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfvén wave turbulence. The injection of Alfvén wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimensional magnetic field topology is simulated using data from photospheric magnetic field measurements. This model does not impose open-closed magnetic field boundaries; those develop self-consistently. The physics include the following. (1) The model employs three different temperatures, namely the isotropic electron temperature and the parallel and perpendicular ion temperatures. The firehose, mirror, and ion-cyclotron instabilities due to the developing ion temperature anisotropy are accounted for. (2) The Alfvén waves are partially reflected by the Alfvén speed gradient and the vorticity along the field lines. The resulting counter-propagating waves are responsible for the nonlinear turbulent cascade. The balanced turbulence due to uncorrelated waves near the apex of the closed field lines and the resulting elevated temperatures are addressed. (3) To apportion the wave dissipation to the three temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating. (4) We have incorporated the collisional and collisionless electron heat conduction. We compare the simulated multi-wavelength extreme ultraviolet images of CR2107 with the observations from STEREO/EUVI and the Solar Dynamics Observatory/AIA instruments. We demonstrate that the reflection due to strong magnetic fields in the proximity of active regions sufficiently intensifies the dissipation and observable emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Observation of N 15 + α resonant structures in F 19 using the thick target in inverse kinematics scattering method

The spectroscopy of F 19 is of interest for nuclear astrophysics and nuclear structure. In astrophysics, fluorine and the reactions producing and destroying it play a key role in constraining models of stars in different evolutionary stages, such as the asymptotic giant branch (AGB) stars, responsible of the production of about half of the elements heavier than Fe. In nuclear structure, F 19 has been subject to many investigations aiming at the identification of α and more exotic cluster structures. Also, its spectroscopy is very useful to constrain the nuclear properties of the 19 Ne mirror nucleus. Here in this work, we report on the measurement of the 15 N α elastic scattering using the thick target inverse kinematics approach, allowing us to span a very large fluorine excitation energy range (~6-10 MeV). The use of 15 N α scattering proves very useful to study α clustering in F 19 thanks to the likelihood for populating states with such a structure. Indeed, the R-matrix analysis of the measured differential cross sections shows the occurrence of many candidate α-cluster states of F 19 . It also calls for the redefinition of the spin-parity and widths of a number of F 19 states with respect to what reported in the literature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

A GRAVITATIONAL DOUBLE-SCATTERING MECHANISM FOR GENERATING HIGH-VELOCITY OBJECTS DURING HALO MERGERS

We present a dynamical model that describes how halo particles can receive a significant energy kick from the merger between their own host halo and a target halo. This could provide a possible explanation for some high-velocity objects, including extended systems like globular clusters (GCs). In the model we especially introduce a double-scattering mechanism, where a halo particle receives a significant part of its total energy kick by first undergoing a gravitational deflection by the target halo and subsequently by its original host halo. This generates an energy kick that is due to the relative velocity between the halos during the deflections. We derive analytically the total kick energy of the particle, which is composed of energy from the double-scattering mechanism and tidal fields, as a function of its position in its original host halo just before merger. In the case of a 1:10 merger, we find that the presented mechanisms can easily generate particles with a velocity approximately two times the virial velocity of the target halo. This motivates us to suggest that the high velocity of the recently discovered GC HVGC-1 can be explained by a head-on halo merger. Finally, we illustrate the orbital evolution of high-velocity particles outside the virial sphere of the target halo by solving the equation of motion in an expanding universe. We find a sweet spot around a scale factor of 0.3-0.5 for ejecting particles into large orbits, which can easily reach beyond approximately five virial radii.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

Radiography in High Mass X-Ray Binaries: Micro-Structure of the Stellar Wind Through Variability of the Column Density

Context. In high mass X-ray binaries, an accreting compact object orbits a high mass star, which loses mass through a dense and in homogeneous wind. Aims. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass-loss rates for high mass stars. Methods. We developed a simplified representation of the stellar wind where all the matter is gathered in spherical “clumps” that are radially advected away from the star. This model enables us to explore the connections between the stochastic properties of the wind and the variability of the column density for a comprehensive set of parameters related to the orbit and to the wind micro-structure, such as the size of the clumps and their individual mass. In particular, we focus on the evolution with the orbital phase of the standard deviation of the column density and of the characteristic duration of enhanced absorption episodes. Using the porosity length, we derive analytical predictions and compare them to the standard deviations and coherence time scales that were obtained. Results. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a single clump in front of the compact object. We thus provide a procedure to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. Conclusions. The coherence time scale grants direct access to the size of the clumps, while their mass can be deduced separately from the amplitude of the variability. We further show how monitoring the variability at superior conjunctions can probe the onset of the clump-forming region above the stellar photosphere. If the high column density variations in some high mass X-ray binaries are due to unaccreted clumps which are passing by the line-of-sight, this would require high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are marginally compatible with the ones derived from radiative-hydrodynamics simulations. Alternatively, the following components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism that has yet to be identified or a dense environment in the immediate vicinity of the accretor, such as an accretion disk, an outflow, or a spherical shell surrounding the magnetosphere of the accreting neutron star.

stars: mass loss↗
DOE OSTI2017

SYNTHETIC OBSERVATIONS OF MAGNETIC FIELDS IN PROTOSTELLAR CORES

The role of magnetic fields in the early stages of star formation is not well constrained. In order to discriminate between different star formation models, we analyze 3D magnetohydrodynamic simulations of low-mass cores and explore the correlation between magnetic field orientation and outflow orientation over time. We produce synthetic observations of dust polarization at resolutions comparable to millimeter-wave dust polarization maps observed by the Combined Array for Research in Millimeter-wave Astronomy and compare these with 2D visualizations of projected magnetic field and column density. Cumulative distribution functions of the projected angle between the magnetic field and outflow show different degrees of alignment in simulations with differing mass-to-flux ratios. The distribution function for the less magnetized core agrees with observations finding random alignment between outflow and field orientations, while the more magnetized core exhibits stronger alignment. We find that fractional polarization increases when the system is viewed such that the magnetic field is close to the plane of the sky, and the values of fractional polarization are consistent with observational measurements. The simulation outflow, which reflects the underlying angular momentum of the accreted gas, changes direction significantly over over the first ∼0.1 Myr of evolution. This movement could lead to the observed random alignment between outflows and the magnetic fields in protostellar cores.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SUPERNOVA DRIVING. II. COMPRESSIVE RATIO IN MOLECULAR-CLOUD TURBULENCE

The compressibility of molecular cloud (MC) turbulence plays a crucial role in star formation models, because it controls the amplitude and distribution of density fluctuations. The relation between the compressive ratio (the ratio of powers in compressive and solenoidal motions) and the statistics of turbulence has been previously studied systematically only in idealized simulations with random external forces. In this work, we analyze a simulation of large-scale turbulence (250 pc) driven by supernova (SN) explosions that has been shown to yield realistic MC properties. We demonstrate that SN driving results in MC turbulence with a broad lognormal distribution of the compressive ratio, with a mean value ≈0.3, lower than the equilibrium value of ≈0.5 found in the inertial range of isothermal simulations with random solenoidal driving. We also find that the compressibility of the turbulence is not noticeably affected by gravity, nor are the mean cloud radial (expansion or contraction) and solid-body rotation velocities. Furthermore, the clouds follow a general relation between the rms density and the rms Mach number similar to that of supersonic isothermal turbulence, though with a large scatter, and their average gas density probability density function is described well by a lognormal distribution, with the addition of a high-density power-law tail when self-gravity is included.

79 ASTRONOMY AND ASTROPHYSICS↗