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,495 records · Page 83

Gravitational waves and stalled satellites from massive galaxy mergers at z ≤ 1

We present a model for merger-driven evolution of the mass function for massive galaxies and their central supermassive black holes at late times. We discuss the current observational evidence in favor of merger-driven massive galaxy evolution during this epoch, and demonstrate that the observed evolution of the mass function can be reproduced by evolving an initial mass function under the assumption of negligible star formation. We calculate the stochastic gravitational wave signal from the resulting black hole binary mergers in the low redshift universe (z ≤ 1) implied by this model, and find that this population has a signal-to-noise ratio 2 × to 5 × larger than previous estimates for pulsar timing arrays, with a (2σ, 3σ) lower limit within this model of h{sub c}(f = 1 yr{sup –1}) = (1.1 × 10{sup –15}, 6.8 × 10{sup –16}). The strength of this signal is sufficient to make it detectable with high probability under conservative assumptions within the next several years. A principle reason that this result is larger than previous estimates is our use of a recent recalibration of the black hole-stellar mass correlation for the brightest cluster galaxies, which increases our estimate by a factor of ∼2 relative to past results. For cases where a galaxy merger fails to lead to a black hole merger, we estimate the probability for a given number of satellite black holes to remain within a massive host galaxy, and interpret the result in light of ULX observations. We find that in rare cases, wandering supermassive black holes may be bright enough to appear as ULXs.

79 ASTRONOMY AND ASTROPHYSICS↗

On the spin-axis dynamics of a Moonless Earth

The variation of a planet's obliquity is influenced by the existence of satellites with a high mass ratio. For instance, Earth's obliquity is stabilized by the Moon and would undergo chaotic variations in the Moon's absence. In turn, such variations can lead to large-scale changes in the atmospheric circulation, rendering spin-axis dynamics a central issue for understanding climate. The relevant quantity for dynamically forced climate change is the rate of chaotic diffusion. Accordingly, here we re-examine the spin-axis evolution of a Moonless Earth within the context of a simplified perturbative framework. We present analytical estimates of the characteristic Lyapunov coefficient as well as the chaotic diffusion rate and demonstrate that even in absence of the Moon, the stochastic change in Earth's obliquity is sufficiently slow to not preclude long-term habitability. Our calculations are consistent with published numerical experiments and illustrate the putative system's underlying dynamical structure in a simple and intuitive manner.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear dynamics of accretion disks with stochastic viscosity

We present a nonlinear numerical model for a geometrically thin accretion disk with the addition of stochastic nonlinear fluctuations in the viscous parameter. These numerical realizations attempt to study the stochastic effects on the disk angular momentum transport. We show that this simple model is capable of reproducing several observed phenomenologies of accretion-driven systems. The most notable of these is the observed linear rms-flux relationship in the disk luminosity. This feature is not formally captured by the linearized disk equations used in previous work. A Fourier analysis of the dissipation and mass accretion rates across disk radii show coherence for frequencies below the local viscous frequency. This is consistent with the coherence behavior observed in astrophysical sources such as Cygnus X-1.

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic microhertz gravitational radiation from stellar convection

High Reynolds-number turbulence driven by stellar convection in main-sequence stars generates stochastic gravitational radiation. We calculate the wave-strain power spectral density as a function of the zero-age main-sequence mass for an individual star and for an isotropic, universal stellar population described by the Salpeter initial mass function and redshift-dependent Hopkins-Beacom star formation rate. The spectrum is a broken power law, which peaks near the turnover frequency of the largest turbulent eddies. The signal from the Sun dominates the universal background. For the Sun, the far-zone power spectral density peaks at S(f {sub peak}) = 5.2 × 10{sup –52} Hz{sup –1} at frequency f {sub peak} = 2.3 × 10{sup –7} Hz. However, at low observing frequencies f < 3 × 10{sup –4} Hz, the Earth lies inside the Sun's near zone and the signal is amplified to S {sub near}(f {sub peak}) = 4.1 × 10{sup –27} Hz{sup –1} because the wave strain scales more steeply with distance (∝d {sup –5}) in the near zone than in the far zone (∝d {sup –1}). Hence the Solar signal may prove relevant for pulsar timing arrays. Other individual sources and the universal background fall well below the projected sensitivities of the Laser Interferometer Space Antenna and next-generation pulsar timing arrays. Stellar convection sets a fundamental noise floor for more sensitive stochastic gravitational-wave experiments in the more distant future.

79 ASTRONOMY AND ASTROPHYSICS↗

Self-sustained asymmetry of lepton-number emission: a new phenomenon during the supernova shock-accretion phase in three dimensions

During the stalled-shock phase of our three-dimensional, hydrodynamical core-collapse simulations with energy-dependent, three-flavor neutrino transport, the lepton-number flux (ν {sub e} minus ν-bar {sub e}) emerges predominantly in one hemisphere. This novel, spherical-symmetry breaking neutrino-hydrodynamical instability is termed LESA for 'Lepton-number Emission Self-sustained Asymmetry'. While the individual ν {sub e} and ν-bar {sub e} fluxes show a pronounced dipole pattern, the heavy-flavor neutrino fluxes and the overall luminosity are almost spherically symmetric. Initially, LESA seems to develop stochastically from convective fluctuations. It exists for hundreds of milliseconds or more and persists during violent shock sloshing associated with the standing accretion shock instability. The ν {sub e} minus ν-bar {sub e} flux asymmetry originates predominantly below the neutrinosphere in a region of pronounced proto-neutron star (PNS) convection, which is stronger in the hemisphere of enhanced lepton-number flux. On this side of the PNS, the mass accretion rate of lepton-rich matter is larger, amplifying the lepton-emission asymmetry, because the spherical stellar infall deflects on a dipolar deformation of the stalled shock. The increased shock radius in the hemisphere of less mass accretion and minimal lepton-number flux ( ν-bar {sub e} flux maximum) is sustained by stronger convection on this side, which is boosted by stronger neutrino heating due to 〈ϵ{sub ν-bar{sub e}}〉>〈ϵ{sub ν{sub e}}〉. Asymmetric heating thus supports the global deformation despite extremely nonstationary convective overturn behind the shock. While these different elements of the LESA phenomenon form a consistent picture, a full understanding remains elusive at present. There may be important implications for neutrino-flavor oscillations, the neutron-to-proton ratio in the neutrino-heated supernova ejecta, and neutron-star kicks, which remain to be explored.

79 ASTRONOMY AND ASTROPHYSICS↗

Galactic cosmic-ray modulation near the heliopause

We investigate the modulation of galactic cosmic rays in the inner and outer heliosheaths using three-dimensional numerical simulations. The model is based on the Parker transport equation integrated using a stochastic phase-space trajectory method. Integration is performed on a plasma background obtained from a global three-dimensional magnetohydrodynamic simulations. Our results predict a negligible amount of modulation in the outer heliosheath because of weak scattering of cosmic ray ions owing to very low levels of magnetic fluctuation power at wavenumbers relevant to the transport of cosmic rays with MeV to GeV energies. This means that the heliopause may be treated as a Dirichlet-type boundary for the purpose of energetic particle modeling. We present models with and without drift velocity to facilitate comparison with papers published earlier. We also attempt to reproduce the sudden step-like increases of cosmic-ray intensity observed by Voyager 1 before its encounter with the heliopause. Our results indicate that very slow cross-field diffusion in the outer heliosheath could produce a large gradient of cosmic rays inside the heliospheric boundary. The resulting large gradient in cosmic-ray intensity near the heliopause qualitatively agrees with recent Voyager 1 observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Big Fish in Small Ponds: massive stars in the low-mass clusters of M83

We have used multi-wavelength Hubble Space Telescope WFC3 data of the starbursting spiral galaxy M83 in order to measure variations in the upper end of the stellar initial mass function (uIMF) using the production rate of ionizing photons in unresolved clusters with ages ≤ 8 Myr. As in earlier papers on M51 and NGC 4214, the uIMF in M83 is consistent with a universal IMF, and stochastic sampling of the stellar populations in the ∼<10{sup 3} M {sub ☉} clusters are responsible for any deviations in this universality. The ensemble cluster population, as well as individual clusters, also imply that the most massive star in a cluster does not depend on the cluster mass. In fact, we have found that these small clusters seem to have an over-abundance of ionizing photons when compared to an expected universal or truncated IMF. This also suggests that the presence of massive stars in these clusters does not affect the star formation in a destructive way.

79 ASTRONOMY AND ASTROPHYSICS↗

Implications of the delayed 2013 outburst of ESO 243-49 HLX-1

After showing four quasi-periodic outbursts spaced by ∼1 yr from 2009 to 2012, the hyper luminous X-ray source ESO 243-49 HLX-1, currently the best intermediate mass black hole (IMBH) candidate, showed an outburst in 2013 delayed by more than a month. In Lasota et al., we proposed that the X-ray light curve is the result of enhanced mass transfer episodes at periapsis from a donor star orbiting the IMBH in a highly eccentric orbit. In this scenario, the delay can be explained only if the orbital parameters can change suddenly from orbit to orbit. To investigate this, we ran Newtonian smooth particle hydrodynamical simulations starting with an incoming donor approaching an IMBH on a parabolic orbit. We survey a large parameter space by varying the star-to-black hole mass ratio (10{sup –5}-10{sup –3}) and the periapsis separation r{sub p} from 2.2 to 2.7r{sub t} with r{sub t} , the tidal radius. To model the donor, we choose several polytropes (Γ = 5/2, n = 3/2, Γ = 3/2, n = 2, Γ = 5/3, n = 2, and Γ = 5/3, n = 3). Once the system is formed, the orbital period decreases until reaching a minimum that may be shallow. Then, the period tends to increase over several periapsis passages due to tidal effects and increasing mass transfer, leading ultimately to the ejection of the donor. We show that the development of stochastic fluctuations inside the donor by adding or removing orbital energy from the system could lead to sudden changes in the orbital period from orbit to orbit with the appropriate order of magnitude to that which has been observed for HLX-1. We also show that given the constraints on the black hole (BH) mass (M {sub BH} > 10{sup 4} M {sub ☉}) and assuming that the HLX-1 system is currently near a minimum in period of ∼1 yr, the donor has to be a white dwarf or a stripped giant core. We predict that if HLX-1 is indeed emerging from a minimum in orbital period, then the period would generally increase with each passage, although substantial stochastic fluctuations can be superposed on this trend.

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic acceleration of electrons by fast magnetosonic waves in solar flares: the effects of anisotropy in velocity and wavenumber space

We develop a model for stochastic acceleration of electrons in solar flares. As in several previous models, the electrons are accelerated by turbulent fast magnetosonic waves ({sup f}ast waves{sup )} via transit-time-damping (TTD) interactions. (In TTD interactions, fast waves act like moving magnetic mirrors that push the electrons parallel or anti-parallel to the magnetic field). We also include the effects of Coulomb collisions and the waves' parallel electric fields. Unlike previous models, our model is two-dimensional in both momentum space and wavenumber space and takes into account the anisotropy of the wave power spectrum F{sub k} and electron distribution function f {sub e}. We use weak turbulence theory and quasilinear theory to obtain a set of equations that describes the coupled evolution of F{sub k} and f {sub e}. We solve these equations numerically and find that the electron distribution function develops a power-law-like non-thermal tail within a restricted range of energies E in (E {sub nt}, E {sub max}). We obtain approximate analytic expressions for E {sub nt} and E {sub max}, which describe how these minimum and maximum energies depend upon parameters such as the electron number density and the rate at which fast-wave energy is injected into the acceleration region at large scales. We contrast our results with previous studies that assume that F{sub k} and f {sub e} are isotropic, and we compare one of our numerical calculations with the time-dependent hard-X-ray spectrum observed during the 1980 June 27 flare. In our numerical calculations, the electron energy spectra are softer (steeper) than in models with isotropic F{sub k} and f {sub e} and closer to the values inferred from observations of solar flares.

79 ASTRONOMY AND ASTROPHYSICS↗

The role of turbulent pressure as a coherent pulsational driving mechanism: the case of the δ Scuti star HD 187547

HD 187547 was the first candidate that led to the suggestion that solar-like oscillations are present in δ Scuti stars. Longer observations, however, show that the modes interpreted as solar-like oscillations have either very long mode lifetimes, longer than 960 days, or are coherent. These results are incompatible with the nature of 'pure' stochastic excitation as observed in solar-like stars. Nonetheless, one point is certain: the opacity mechanism alone cannot explain the oscillation spectrum of HD 187547. Here we present new theoretical investigations showing that convection dynamics can intrinsically excite coherent pulsations in the chemically peculiar δ Scuti star HD 187547. More precisely, it is the perturbations of the mean Reynold stresses (turbulent pressure) that drives the pulsations and the excitation takes place predominantly in the hydrogen ionization zone.

79 ASTRONOMY AND ASTROPHYSICS↗

The formation of kappa-distribution accelerated electron populations in solar flares

Driven by recent RHESSI observations of confined loop-top hard X-ray sources in solar flares, we consider stochastic acceleration of electrons in the presence of Coulomb collisions. If electron escape from the acceleration region can be neglected, the electron distribution function is determined by a balance between diffusive acceleration and collisions. Such a scenario admits a stationary solution for the electron distribution function that takes the form of a kappa distribution. We show that the evolution toward this kappa distribution involves a 'wave front' propagating forward in velocity space, so that electrons of higher energy are accelerated later; the acceleration timescales with energy according to τ{sub acc} ∼ E {sup 3/2}. At sufficiently high energies escape from the finite-length acceleration region will eventually dominate. For such energies, the electron velocity distribution function is obtained by solving a time-dependent Fokker-Planck equation in the 'leaky-box' approximation. Solutions are obtained in the limit of a small escape rate from an acceleration region that can effectively be considered a thick target.

79 ASTRONOMY AND ASTROPHYSICS↗

FORMATION OF MULTIPLE-SATELLITE SYSTEMS FROM LOW-MASS CIRCUMPLANETARY PARTICLE DISKS

Circumplanetary particle disks would be created in the late stage of planetary formation either by impacts of planetary bodies or disruption of satellites or passing bodies, and satellites can be formed by accretion of disk particles spreading across the Roche limit. Previous N-body simulation of lunar accretion focused on the formation of single-satellite systems from disks with large disk-to-planet mass ratios, while recent models of the formation of multiple-satellite systems from disks with smaller mass ratios do not take account of gravitational interaction between formed satellites. In the present work, we investigate satellite accretion from particle disks with various masses, using N-body simulation. In the case of accretion from somewhat less massive disks than the case of lunar accretion, formed satellites are not massive enough to clear out the disk, but can become massive enough to gravitationally shepherd the disk outer edge and start outward migration due to gravitational interaction with the disk. When the radial location of the 2:1 mean motion resonance of the satellite reaches outside the Roche limit, the second satellite can be formed near the disk outer edge, and then the two satellites continue outward migration while being locked in the resonance. Co-orbital satellites are found to be occasionally formed on the orbit of the first satellite. Our simulations also show that stochastic nature involved in gravitational interaction and collision between aggregates in the tidal environment can lead to diversity in the final mass and orbital architecture, which would be expected in satellite systems of exoplanets.

79 ASTRONOMY AND ASTROPHYSICS↗

THE BLACK HOLE FORMATION PROBABILITY

A longstanding question in stellar evolution is which massive stars produce black holes (BHs) rather than neutron stars (NSs) upon death. It has been common practice to assume that a given zero-age main sequence (ZAMS) mass star (and perhaps a given metallicity) simply produces either an NS or a BH, but this fails to account for a myriad of other variables that may effect this outcome, such as spin, binarity, or even stochastic differences in the stellar structure near core collapse. We argue that instead a probabilistic description of NS versus BH formation may be better suited to account for the current uncertainties in understanding how massive stars die. We present an initial exploration of the probability that a star will make a BH as a function of its ZAMS mass, P {sub BH}(M {sub ZAMS}). Although we find that it is difficult to derive a unique P {sub BH}(M {sub ZAMS}) using current measurements of both the BH mass distribution and the degree of chemical enrichment by massive stars, we demonstrate how P {sub BH}(M {sub ZAMS}) changes with these various observational and theoretical uncertainties. We anticipate that future studies of Galactic BHs and theoretical studies of core collapse will refine P {sub BH}(M {sub ZAMS}) and argue that this framework is an important new step toward better understanding BH formation. A probabilistic description of BH formation will be useful as input for future population synthesis studies that are interested in the formation of X-ray binaries, the nature and event rate of gravitational wave sources, and answering questions about chemical enrichment.

79 ASTRONOMY AND ASTROPHYSICS↗

A KINETIC TRANSPORT THEORY FOR PARTICLE ACCELERATION AND TRANSPORT IN REGIONS OF MULTIPLE CONTRACTING AND RECONNECTING INERTIAL-SCALE FLUX ROPES

Simulations of particle acceleration in turbulent plasma regions with multiple contracting and merging (reconnecting) magnetic islands emphasize the key role of temporary particle trapping in island structures for the efficient acceleration of particles to form hard power-law spectra. Statistical kinetic transport theories have been developed that capture the essential physics of particle acceleration in multi-island regions. The transport theory of Zank et al. is further developed by considering the acceleration effects of both the mean and the variance of the electric fields induced by the dynamics of multiple inertial-scale flux ropes. A focused transport equation is derived that includes new Fokker-Planck terms for particle scattering and stochastic acceleration due to the variance in multiple flux-rope magnetic fields, plasma flows, and reconnection electric fields. A Parker transport equation is also derived in which a new expression for momentum diffusion appears, combining stochastic acceleration by particle scattering in the mean multi-flux-rope electric fields with acceleration by the variance in these electric fields. Test particle acceleration is modeled analytically considering drift acceleration by the variance in the induced electric fields of flux ropes in the slow supersonic, radially expanding solar wind. Hard power-law spectra occur for sufficiently strong inertial-scale flux ropes with an index modified by adiabatic cooling, solar wind advection, and diffusive escape from flux ropes. Flux ropes might be sufficiently strong behind interplanetary shocks where the index of suprathermal ion power-law spectra observed in the supersonic solar wind can be reproduced.

79 ASTRONOMY AND ASTROPHYSICS↗

SDSS J14584479+3720215: A BENCHMARK JHK{sub S} BLAZAR LIGHT CURVE FROM THE 2MASS CALIBRATION SCANS

Active galactic nuclei (AGNs) are well-known to exhibit flux variability across a wide range of wavelength regimes, but the precise origin of the variability at different wavelengths remains unclear. To investigate the relatively unexplored near-IR (NIR) variability of the most luminous AGNs, we conduct a search for variability using well sampled JHK{sub s}-band light curves from the Two Micron All Sky Survey (2MASS) calibration fields. Our sample includes 27 known quasars with an average of 924 epochs of observation over three years, as well as one spectroscopically confirmed blazar (SDSS J14584479+3720215) with 1972 epochs of data. This is the best-sampled NIR photometric blazar light curve to date, and it exhibits correlated, stochastic variability that we characterize with continuous auto-regressive moving average (CARMA) models. None of the other 26 known quasars had detectable variability in the 2MASS bands above the photometric uncertainty. A blind search of the 2MASS calibration field light curves for active galactic nucleus (AGN) candidates based on fitting CARMA(1,0) models (damped-random walk) uncovered only seven candidates. All seven were young stellar objects within the ρ Ophiuchus star forming region, five with previous X-ray detections. A significant γ-ray detection (5σ) for the known blazar using 4.5 yr of Fermi photon data is also found. We suggest that strong NIR variability of blazars, such as seen for SDSS J14584479+3720215, can be used as an efficient method of identifying previously unidentified γ-ray blazars, with low contamination from other AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗

COUPLED SPIN AND SHAPE EVOLUTION OF SMALL RUBBLE-PILE ASTEROIDS: SELF-LIMITATION OF THE YORP EFFECT

We present the first self-consistent simulations of the coupled spin-shape evolution of small gravitational aggregates under the influence of the YORP effect. Because of YORP’s sensitivity to surface topography, even small centrifugally driven reconfigurations of aggregates can alter the YORP torque dramatically, resulting in spin evolution that can differ qualitatively from the rigid-body prediction. One-third of our simulations follow a simple evolution described as a modified YORP cycle. Two-thirds exhibit one or more of three distinct behaviors—stochastic YORP, self-governed YORP, and stagnating YORP—which together result in YORP self-limitation. Self-limitation confines rotation rates of evolving aggregates to far narrower ranges than those expected in the classical YORP cycle, greatly prolonging the times over which objects can preserve their sense of rotation. Simulated objects are initially randomly packed, disordered aggregates of identical spheres in rotating equilibrium, with low internal angles of friction. Their shape evolution is characterized by rearrangement of the entire body, including the deep interior. They do not evolve to axisymmetric top shapes with equatorial ridges. Mass loss occurs in one-third of the simulations, typically in small amounts from the ends of a prolate-triaxial body. We conjecture that YORP self-limitation may inhibit formation of top-shapes, binaries, or both, by restricting the amount of angular momentum that can be imparted to a deformable body. Stochastic YORP, in particular, will affect the evolution of collisional families whose orbits drift apart under the influence of Yarkovsky forces, in observable ways.

79 ASTRONOMY AND ASTROPHYSICS↗

SPECTRAL EVOLUTION OF ANOMALOUS COSMIC RAYS AT VOYAGER 1 BEYOND THE TERMINATION SHOCK

When the Voyager 1 spacecraft crossed the termination shock (TS) on 2004 December 16, the energy spectra of anomalous cosmic rays (ACRs) could not have been produced by steady-state diffusive shock acceleration. However, over the next few years, in the declining phase of the solar cycle, the spectra began to evolve into the expected power-law profile. Observations at the shock led to a broad range of alternative theories for ACR acceleration. In spite of that, in this work we show that the observations could be explained by assuming ACRs are accelerated at the TS. In this paper, we propose that the solar cycle had an important effect on the unrolling of the spectra in the heliosheath. To investigate the spectral evolution of ACRs, a magnetohydrodynamic background model with stationary solar-wind inner boundary conditions was used to model the transport of helium and oxygen ions. We used a backward-in-time stochastic integration technique where phase-space trajectories are integrated until the so-called “injection energy” is reached. Our simulation results were compared with Voyager 1 observations using three different diffusion models. It is shown that the spectral evolution of ACRs in the heliosheath at Voyager 1 could be explained by an increase in the source strength and an enhancement in diffusion as a result of a decrease of the turbulent correlation length in the declining phase of the solar cycle. At the same time, drift effects seem to have had a smaller effect on the evolution of the spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Metal-poor, strongly star-forming galaxies in the DEEP2 Survey: The relationship between stellar mass, temperature-based metallicity, and star formation rate

We report on the discovery of 28 z≈0.8 metal-poor galaxies in DEEP2. These galaxies were selected for their detection of the weak [O iii] λ4363 emission line, which provides a “direct” measure of the gas-phase metallicity. A primary goal for identifying these rare galaxies is to examine whether the fundamental metallicity relation (FMR) between stellar mass, gas metallicity, and star formation rate (SFR) holds for low stellar mass and high SFR galaxies. The FMR suggests that higher SFR galaxies have lower metallicity (at fixed stellar mass). To test this trend, we combine spectroscopic measurements of metallicity and dust-corrected SFR with stellar mass estimates from modeling the optical photometry. We find that these galaxies are 1.05 ± 0.61 dex above the z∼1 stellar mass–SFR relation and 0.23 ± 0.23 dex below the local mass–metallicity relation. Relative to the FMR, the latter offset is reduced to 0.01 dex, but significant dispersion remains (0.29 dex with 0.16 dex due to measurement uncertainties). This dispersion suggests that gas accretion, star formation, and chemical enrichment have not reached equilibrium in these galaxies. This is evident by their short stellar mass doubling timescale of ≈100{sub −75}{sup +310} Myr, which suggests stochastic star formation. Combining our sample with other z∼1 metal-poor galaxies, we find a weak positive SFR–metallicity dependence (at fixed stellar mass) that is significant at 94.4% confidence. We interpret this positive correlation as recent star formation that has enriched the gas but has not had time to drive the metal-enriched gas out with feedback mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗