Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stochastic processes”

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 1,405 records · Page 78

NEAR-INFRARED PERIODIC AND OTHER VARIABLE FIELD STARS IN THE FIELD OF THE CYGNUS OB7 STAR-FORMING REGION

We present a subset of the results of a three-season, 124 night, near-infrared monitoring campaign of the dark clouds Lynds 1003 and Lynds 1004 in the Cygnus OB7 star-forming region. In this paper, we focus on the field star population. Using three seasons of UKIRT J, H, and K-band observations spanning 1.5 years, we obtained high-quality photometry on 9200 stars down to J = 17 mag, with photometric uncertainty better than 0.04 mag. After excluding known disk-bearing stars we identify 149 variables-1.6% of the sample. Of these, about 60 are strictly periodic, with periods predominantly <2 days. We conclude this group is dominated by eclipsing binaries. A few stars have long period signals of between 20 and 60 days. About 25 stars have weak modulated signals, but it was not clear if these were periodic. Some of the stars in this group may be diskless young stellar objects with relatively large variability due to cool starspots. The remaining {approx}60 stars showed variations which appear to be purely stochastic.

79 ASTRONOMY AND ASTROPHYSICS↗

Andromeda (M31) optical and infrared disk survey. I. Insights in wide-field near-IR surface photometry

We present wide-field near-infrared J and K{sub s} images of the Andromeda Galaxy (M31) taken with WIRCam at the Canada-France-Hawaii Telescope as part of the Andromeda Optical and Infrared Disk Survey. This data set allows simultaneous observations of resolved stars and near-infrared (NIR) surface brightness across M31's entire bulge and disk (within R = 22 kpc), permitting a direct test of the stellar composition of near-infrared light in a nearby galaxy. Here we develop NIR observation and reduction methods to recover a uniform surface brightness map across the 3° × 1° disk of M31 with 27 WIRCam fields. Two sky-target nodding strategies are tested, and we find that strictly minimizing sky sampling latency cannot improve background subtraction accuracy to better than 2% of the background level due to spatio-temporal variations in the NIR skyglow. We fully describe our WIRCam reduction pipeline and advocate using flats built from night-sky images over a single night, rather than dome flats that do not capture the WIRCam illumination field. Contamination from scattered light and thermal background in sky flats has a negligible effect on the surface brightness shape compared to the stochastic differences in background shape between sky and galaxy disk fields, which are ∼0.3% of the background level. The most dramatic calibration step is the introduction of scalar sky offsets to each image that optimizes surface brightness continuity. Sky offsets reduce the mean surface brightness difference between observation blocks from 1% to <0.1% of the background level, though the absolute background level remains statistically uncertain to 0.15% of the background level. We present our WIRCam reduction pipeline and performance analysis to give specific recommendations for the improvement of NIR wide-field imaging methods.

79 ASTRONOMY AND ASTROPHYSICS↗

REGULARITY OF THE NORTH–SOUTH ASYMMETRY OF SOLAR ACTIVITY: REVISITED

Extended time series of Solar Activity Indices (ESAI) extended the Greenwich series of sunspot area from the year 1874 back to 1821. The ESAI's yearly sunspot area in the northern and southern hemispheres from 1821 to 2013 is utilized to investigate characteristics of the north–south hemispherical asymmetry of sunspot activity. Periodical behavior of about 12 solar cycles is also confirmed from the ESAI data set to exist in dominant hemispheres, linear regression lines of yearly asymmetry values, and cumulative counts of yearly sunspot areas in the hemispheres for solar cycles. The period is also inferred to appear in both the cumulative difference in the yearly sunspot areas in the hemispheres over the entire time interval and in its statistical Student's t-test. The hemispherical bias of sunspot activity should be regarded as an impossible stochastic phenomenon over a long time period.

79 ASTRONOMY AND ASTROPHYSICS↗

EFFECTS OF ROTATION ON STOCHASTICITY OF GRAVITATIONAL WAVES IN THE NONLINEAR PHASE OF CORE-COLLAPSE SUPERNOVAE

By performing three-dimensional (3D) simulations that demonstrate the neutrino-driven core-collapse supernovae aided by the standing accretion shock instability (SASI), we study how the spiral modes of the SASI can impact the properties of the gravitational-wave (GW) emission. To see the effects of rotation in the nonlinear postbounce phase, we give a uniform rotation on the flow advecting from the outer boundary of the iron core, the specific angular momentum of which is assumed to agree with recent stellar evolution models. We compute fifteen 3D models in which the initial angular momentum and the input neutrino luminosities from the protoneutron star are changed in a systematic manner. By performing a ray-tracing analysis, we accurately estimate the GW amplitudes generated by anisotropic neutrino emission. Our results show that the gravitational waveforms from neutrinos in models that include rotation exhibit a common feature; otherwise, they vary much more stochastically in the absence of rotation. The breaking of the stochasticity stems from the excess of the neutrino emission parallel to the spin axis. This is because the compression of matter is more enhanced in the vicinity of the equatorial plane due to the growth of the spiral SASI modes, leading to the formation of the spiral flows circulating around the spin axis with higher temperatures. We point out that recently proposed future space interferometers like Fabry-Perot-type DECIGO would permit the detection of these signals for a Galactic supernova.

79 ASTRONOMY AND ASTROPHYSICS↗

MODELING OF THE INITIAL MASS FUNCTION USING THE METROPOLIS-HASTINGS ALGORITHM

A stochastic model has been developed for the hierarchical fragmentation of a molecular cloud. Here, the number of fragments, time between successive fragmentation steps, and mass of a fragment are considered as random variables, and fragment masses are generated using the Metropolis-Hastings algorithm. The resulting mass spectra, computed at different projected distances and taking opacity into consideration, show a signature of mass segregation. The critical mass, mass spectrum, and mass segregation are consistent with the observations of young massive clusters in our Galaxy as well as in external galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

STOCHASTIC GRAVITATIONAL WAVE BACKGROUND FROM COALESCING BINARY BLACK HOLES

We estimate the stochastic gravitational wave (GW) background signal from the field population of coalescing binary stellar mass black holes (BHs) throughout the universe. This study is motivated by recent observations of BH-Wolf-Rayet (WR) star systems and by new estimates in the metallicity abundances of star-forming galaxies that imply BH-BH systems are more common than previously assumed. Using recent analytical results of the inspiral-merger-ringdown waveforms for coalescing binary BH systems, we estimate the resulting stochastic GW background signal. Assuming average quantities for the single source energy emissions, we explore the parameter space of chirp mass and local rate density required for detection by advanced and third-generation interferometric GW detectors. For an average chirp mass of 8.7 M{sub sun}, we find that detection through 3 years of cross-correlation by two advanced detectors will require a rate density, r{sub 0} {>=} 0.5 Mpc{sup -3} Myr{sup -1}. Combining data from multiple pairs of detectors can reduce this limit by up to 40%. Investigating the full parameter space we find that detection could be achieved at rates r{sub 0} {approx} 0.1 Mpc{sup -3} Myr{sup -1} for populations of coalescing binary BH systems with average chirp masses of {approx}15 M{sub sun} which are predicted by recent studies of BH-WR star systems. While this scenario is at the high end of theoretical estimates, cross-correlation of data by two Einstein Telescopes could detect this signal under the condition r{sub 0} {>=} 10{sup -3}Mpc{sup -3} Myr{sup -1}. Such a signal could potentially mask a primordial GW background signal of dimensionless energy density, {Omega}{sub GW} {approx} 10{sup -10}, around the (1-500) Hz frequency range.

79 ASTRONOMY AND ASTROPHYSICS↗

INCORPORATING KINETIC PHYSICS INTO A TWO-FLUID SOLAR-WIND MODEL WITH TEMPERATURE ANISOTROPY AND LOW-FREQUENCY ALFVEN-WAVE TURBULENCE

We develop a one-dimensional solar-wind model that includes separate energy equations for the electrons and protons, proton temperature anisotropy, collisional and collisionless heat flux, and an analytical treatment of low-frequency, reflection-driven, Alfven-wave (AW) turbulence. To partition the turbulent heating between electron heating, parallel proton heating, and perpendicular proton heating, we employ results from the theories of linear wave damping and nonlinear stochastic heating. We account for mirror and oblique firehose instabilities by increasing the proton pitch-angle scattering rate when the proton temperature anisotropy exceeds the threshold for either instability. We numerically integrate the equations of the model forward in time until a steady state is reached, focusing on two fast-solar-wind-like solutions. These solutions are consistent with a number of observations, supporting the idea that AW turbulence plays an important role in the origin of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

MODELING THE EFFECTS OF STAR FORMATION HISTORIES ON H{alpha} AND ULTRAVIOLET FLUXES IN NEARBY DWARF GALAXIES

We consider the effects of non-constant star formation histories (SFHs) on H{alpha} and GALEX far-ultraviolet (FUV) star formation rate (SFR) indicators. Under the assumption of a fully populated Chabrier initial mass function (IMF), we compare the distribution of H{alpha}-to-FUV flux ratios from {approx}1500 simple, periodic model SFHs with observations of 185 galaxies from the Spitzer Local Volume Legacy survey. We find a set of SFH models that are well matched to the data, such that more massive galaxies are best characterized by nearly constant SFHs, while low-mass systems experience burst amplitudes of {approx}30 (i.e., an increase in the SFR by a factor of 30 over the SFR during the inter-burst period), burst durations of tens of Myr, and periods of {approx}250 Myr; these SFHs are broadly consistent with the increased stochastic star formation expected in systems with lower SFRs. We analyze the predicted temporal evolution of galaxy stellar mass, R-band surface brightness, H{alpha}-derived SFR, and blue luminosity, and find that they provide a reasonable match to observed flux distributions. We find that our model SFHs are generally able to reproduce both the observed systematic decline and increased scatter in H{alpha}-to-FUV ratios toward low-mass systems, without invoking other physical mechanisms. We also compare our predictions with those from the Integrated Galactic IMF theory with a constant SFR. We find that while both predict a systematic decline in the observed ratios, only the time variable SFH models are capable of producing the observed population of low-mass galaxies (M{sub *} {approx}< 10{sup 7} M{sub Sun }) with normal H{alpha}-to-FUV ratios. These results demonstrate that a variable IMF alone has difficulty explaining the observed scatter in the H{alpha}-to-FUV ratios. We conclude by considering the limitations of the model SFHs and discuss the use of additional empirical constraints to improve future SFH modeling efforts.

79 ASTRONOMY AND ASTROPHYSICS↗

MAGNETIC FIELD STRENGTH FLUCTUATIONS IN THE HELIOSHEATH: VOYAGER 1 OBSERVATIONS DURING 2009

We analyze the ''microscale fluctuations'' of the magnetic field strength B on a scale of several hours observed by Voyager1 (V1) in the heliosheath during 2009. The microscale fluctuations of B range from coherent to stochastic structures. The amplitude of microscale fluctuations of B during 1 day is measured by the standard deviation (SD) of 48 s averages of B. The distribution of the daily values of SD is lognormal. SD(t) from day of year (DOY) 1 to 331, 2009, is very intermittent. SD(t) has a 1/f or 'pink noise' spectrum on scales from 1 to 100 days, and it has a broad multifractal spectrum f({alpha}) with 0.57 {<=} {alpha} {<=} 1.39. The time series of increments SD(t + {tau}) - SD(t) has a pink noise spectrum with {alpha}' = 0.88 {+-} 0.14 on scales from 1 to 100 days. The increments have a Tsallis (q-Gaussian) distribution on scales from 1 to 165 days, with an average q = 1.75 {+-} 0.12. The skewness S and kurtosis K have Gaussian and lognormal distributions, respectively. The largest spikes in K(t) and S(t) are often associated with a change in B across a data gap and with identifiable physical structures. The 'turbulence' observed by V1 during 2009 was weakly compressible on average but still very intermittent, highly variable, and highly compressible at times. The turbulence observed just behind the termination shock by Voyager 2 was twice as strong. These observations place strong constraints on any model of 'turbulence' in the heliosheath.

79 ASTRONOMY AND ASTROPHYSICS↗

THE PARALLEL DIFFUSION COEFFICIENT AND ADIABATIC FOCUSING OF COSMIC-RAY PARTICLES

In this paper, the problem of focused particle transport is revisited. A description in terms of a system of stochastic differential equations, completely equivalent to the Fokker-Planck equation, is suggested. The coefficient for spatial diffusion parallel to the mean magnetic field is calculated. The case of isotropic pitch angle scattering and weak focusing is analyzed in detail. The disagreement between a recent analysis by Shalchi and other treatments of the same problem is discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

SLUG-STOCHASTICALLY LIGHTING UP GALAXIES. I. METHODS AND VALIDATING TESTS

The effects of stochasticity on the luminosities of stellar populations are an often neglected but crucial element for understanding populations in the low-mass or the low star formation rate regime. To address this issue, we present SLUG, a new code to 'Stochastically Light Up Galaxies'. SLUG synthesizes stellar populations using a Monte Carlo technique that properly treats stochastic sampling including the effects of clustering, the stellar initial mass function, star formation history, stellar evolution, and cluster disruption. This code produces many useful outputs, including (1) catalogs of star clusters and their properties such as their stellar initial mass distributions and their photometric properties in a variety of filters, (2) two dimensional histograms of color-magnitude diagrams of every star in the simulation, and (3) the photometric properties of field stars and the integrated photometry of the entire simulated galaxy. After presenting the SLUG algorithm in detail, we validate the code through comparisons with STARBURST99 in the well-sampled regime, and with observed photometry of Milky Way clusters. Finally, we demonstrate SLUG's capabilities by presenting outputs in the stochastic regime. SLUG is publicly distributed through the Web site http://sites.google.com/site/runslug/.

79 ASTRONOMY AND ASTROPHYSICS↗

THE METALLICITY EVOLUTION OF INTERACTING GALAXIES

Nuclear inflows of metal-poor interstellar gas triggered by galaxy interactions can account for the systematically lower central oxygen abundances observed in local interacting galaxies. Here, we investigate the metallicity evolution of a large set of simulations of colliding galaxies. Our models include cooling, star formation, feedback, and a new stochastic method for tracking the mass recycled back to the interstellar medium from stellar winds and supernovae. We study the influence of merger-induced inflows, enrichment, gas consumption, and galactic winds in determining the nuclear metallicity. The central metallicity is primarily a competition between the inflow of low-metallicity gas and enrichment from star formation. An average depression in the nuclear metallicity of {approx}0.07 is found for gas-poor disk-disk interactions. Gas-rich disk-disk interactions, on the other hand, typically have an enhancement in the central metallicity that is positively correlated with the gas content. The simulations fare reasonably well when compared to the observed mass-metallicity and separation-metallicity relationships, but further study is warranted.

79 ASTRONOMY AND ASTROPHYSICS↗

CHALLENGES IN FORMING PLANETS BY GRAVITATIONAL INSTABILITY: DISK IRRADIATION AND CLUMP MIGRATION, ACCRETION, AND TIDAL DESTRUCTION

We present two-dimensional hydrodynamic simulations of self-gravitating protostellar disks subject to axisymmetric, continuing mass loading from an infalling envelope and irradiation from the central star to explore the growth of gravitational instability (GI) and disk fragmentation. We assume that the disk is built gradually and smoothly by the infall, resulting in good numerical convergence. We confirm that for disks around solar-mass stars, infall at high rates at radii beyond {approx}50 AU leads to disk fragmentation. At lower infall rates, however, irradiation suppresses fragmentation. We find that, once formed, the fragments or clumps migrate inward on typical type I timescales of {approx}2 Multiplication-Sign 10{sup 3} yr initially, but with a stochastic component superimposed due to their interaction with the GI-induced spiral arms. Migration begins to deviate from the type I timescale when the clump becomes more massive than the local disk mass, and/or when the clump begins to form a gap in the disk. As they migrate, clumps accrete from the disk at a rate between 10{sup -3} and 10{sup -1} M{sub J} yr{sup -1}, consistent with analytic estimates that assume a 1-2 Hill radii cross section. The eventual fates of these clumps, however, diverge depending on the migration speed: 3 out of 13 clumps in our simulations become massive enough (brown dwarf mass range) to open gaps in the disk and essentially stop migrating; 4 out of 13 are tidally destroyed during inward migration; and 6 out of 13 migrate across the inner boundary of the simulated disks. A simple analytic model for clump evolution explains the different fates of the clumps and reveals some limitations of two-dimensional simulations. Overall, our results indicate that fast migration, accretion, and tidal destruction of the clumps pose challenges to the scenario of giant planet formation by GI in situ, although we cannot address whether or not remnant solid cores can survive after tidal stripping. The models where the massive clumps are not disrupted and open gaps may be relevant to the formation of close binary systems similar to Kepler 16. A clump formed by GI-induced fragmentation can be as large as 10 AU and as luminous as 2 Multiplication-Sign 10{sup -3} L{sub Sun }, which will be easily detectable with ALMA, but its lifetime before dynamically collapsing is only {approx}1000 years.

79 ASTRONOMY AND ASTROPHYSICS↗

SECULAR GRAVITATIONAL INSTABILITY OF A DUST LAYER IN SHEAR TURBULENCE

We perform a linear stability analysis of a dust layer in a turbulent gas disk. Youdin investigated the secular gravitational instability (GI) of a dust layer using hydrodynamic equations with a turbulent diffusion term. We obtain essentially the same result independently of Youdin. In the present analysis, we restrict the area of interest to small dust particles, while investigating the secular GI in a more rigorous manner. We discuss the time evolution of the dust surface density distribution using a stochastic model and derive the advection-diffusion equation. The validity of the analysis by Youdin is confirmed in the strong drag limit. We demonstrate quantitatively that the finite thickness of a dust layer weakens the secular GI and that the density-dependent diffusion coefficient changes the growth rate. We apply the results obtained to the turbulence driven by the shear instability and find that the secular GI is faster than the radial drift when the gas density is three times as large as that in the minimum-mass disk model. If the dust particles are larger than chondrules, the secular GI grows within the lifetime of a protoplanetary disk.

79 ASTRONOMY AND ASTROPHYSICS↗

THE ROLE OF STOCHASTIC ACCELERATION IN THE PROMPT EMISSION OF GAMMA-RAY BURSTS: APPLICATION TO HADRONIC INJECTION

We study effects of particle re-acceleration (or heating) in the post-shock region via magnetohydrodynamic/plasma turbulence, in the context of a mixed hadronic-leptonic model for the prompt emission of gamma-ray bursts, using both analytical and numerical methods. We show that stochastically accelerated (or heated) leptons, which are injected via pp and p{gamma} reactions and subsequent pair cascades, are plausibly able to reproduce the Band function spectra with {alpha} {approx} 1 and {beta} {approx} 2-3 in the {approx}MeV range. An additional hard component coming from the proton-induced cascade emission is simultaneously expected, which can be compatible with observed extra power-law spectra far above the MeV range. We also discuss the specific implications of hadronic models for ongoing high-energy neutrino observations.

79 ASTRONOMY AND ASTROPHYSICS↗

REVISITING ACCELERATION OF CHARGED GRAINS IN MAGNETOHYDRODYNAMIC TURBULENCE

We study the acceleration of charged grains by magnetohydrodynamic (MHD) turbulence in the interstellar medium (ISM). We begin with revisiting gyroresonance acceleration by taking into account the fluctuations of grain guiding center along a uniform magnetic field (i.e., nonlinear theory-NLT). We calculate grain velocities due to gyroresonance by fast MHD modes using the NLT for different phases of the ISM and compare them with results obtained using quasi-linear theory (QLT). We find for the parameters applicable to the typical ISM phases that the fluctuations of the grain guiding center reduce grain velocities by less than 15%, but they can be important for more special circumstances. We confirm that large grains can be accelerated to super-Alfvenic velocities through gyroresonance. For such super-Alfvenic grains, we investigate the effect of further acceleration via transit-time damping (TTD) by fast modes. We find that due to the broadening of the resonance condition in the NLT, the TTD acceleration is not only important for the cosines of grain pitch angle relative to the magnetic field {mu} > V{sub A}/v, but also for {mu} < V{sub A}/v where v is the grain velocity and V{sub A} is the Alfven speed. We show that the TTD acceleration is dominant over the gyroresonance for large grains, and can increase substantially grain velocities induced by gyroresonance acceleration. We quantify another stochastic acceleration mechanism arising from low-frequency Alfven waves. We discuss the range of applicability of the mechanisms and their implications.

79 ASTRONOMY AND ASTROPHYSICS↗

IMPULSIVE PHASE CORONAL HARD X-RAY SOURCES IN AN X3.9 CLASS SOLAR FLARE

We present the analysis of a pair of unusually energetic coronal hard X-ray (HXR) sources detected by the Reuven Ramaty High Energy Solar Spectroscopic Imager during the impulsive phase of an X3.9 class solar flare on 2003 November 3, which simultaneously shows two intense footpoint (FP) sources. A distinct loop top (LT) coronal source is detected up to {approx}150 keV and a second (upper) coronal source up to {approx}80 keV. These photon energies, which were not fully investigated in earlier analysis of this flare, are much higher than commonly observed in coronal sources and pose grave modeling challenges. The LT source in general appears higher in altitude with increasing energy and exhibits a more limited motion compared to the expansion of the thermal loop. The high-energy LT source shows an impulsive time profile and its nonthermal power-law spectrum exhibits soft-hard-soft evolution during the impulsive phase, similar to the FP sources. The upper coronal source exhibits an opposite spatial gradient and a similar spectral slope compared to the LT source. These properties are consistent with the model of stochastic acceleration of electrons by plasma waves or turbulence. However, the LT and FP spectral index difference (varying from {approx}0 to 1) is much smaller than commonly measured and than that expected from a simple stochastic acceleration model. Additional confinement or trapping mechanisms of high-energy electrons in the corona are required. Comprehensive modeling including both kinetic effects and the macroscopic flare structure may shed light on this behavior. These results highlight the importance of imaging spectroscopic observations of the LT and FP sources up to high energies in understanding electron acceleration in solar flares. Finally, we show that the electrons producing the upper coronal HXR source may very likely be responsible for the type III radio bursts at the decimetric/metric wavelength observed during the impulsive phase of this flare.

79 ASTRONOMY AND ASTROPHYSICS↗

THE TORQUING OF CIRCUMNUCLEAR ACCRETION DISKS BY STARS AND THE EVOLUTION OF MASSIVE BLACK HOLES

An accreting massive black hole (MBH) in a galactic nucleus is surrounded by a dense stellar cluster. We analyze and simulate numerically the evolution of a thin accretion disk due to its internal viscous torques, due to the frame-dragging torques of a spinning MBH (the Bardeen-Petterson effect), and due to the orbit-averaged gravitational torques by the stars (resonant relaxation). We show that the evolution of the MBH mass accretion rate, the MBH spin growth rate, and the covering fraction of the disk relative to the central ionizing continuum source, are all strongly coupled to the stochastic fluctuations of the stellar potential via the warps that the stellar torques excite in the disk. These lead to fluctuations by factors of up to a few in these quantities over a wide range of timescales, with most of the power on timescales {approx}> (M{sub .}/M{sub d} )P(R{sub d} ), where M{sub .} and M{sub d} are the masses of the MBH and disk, and P is the orbital period at the disk's mass-weighted mean radius R{sub d}. The response of the disk is stronger the lighter it is and the more centrally concentrated the stellar cusp. As proof of concept, we simulate the evolution of the low-mass maser disk in NGC 4258 and show that its observed O(10 Degree-Sign ) warp can be driven by the stellar torques. We also show that the frame dragging of a massive active galactic nucleus disk couples the stochastic stellar torques to the MBH spin and can excite a jitter of a few degrees in its direction relative to that of the disk's outer regions.

79 ASTRONOMY AND ASTROPHYSICS↗