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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

STOCHASTIC PARTICLE ACCELERATION IN TURBULENCE GENERATED BY MAGNETOROTATIONAL INSTABILITY

We investigate stochastic particle acceleration in accretion flows. It is believed that magnetorotational instability (MRI) generates turbulence inside accretion flows and that cosmic rays (CRs) are accelerated by the turbulence. We calculate equations of motion for CRs in the turbulent fields generated by MRI with the shearing box approximation and without back reaction to the field. Our results show that the CRs randomly gain or lose their energy through interaction with the turbulent fields. The CRs diffuse in the configuration space anisotropically: the diffusion coefficient in the direction of the unperturbed flow is about 20 times higher than the Bohm coefficient, while those in the other directions are only a few times higher than the Bohm. The momentum distribution is isotropic and its evolution can be described by the diffusion equation in momentum space where the diffusion coefficient is a power-law function of the CR momentum. We show that the shear acceleration works efficiently for energetic particles. We also cautiously note that in the shearing box approximation, particles that cross the simulation box many times along the radial direction undergo unphysical runaway acceleration by the Lorentz transformation, which needs to be taken into account with special care.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Deep mixing in stellar variability: improved method, statistics, and applications

The preliminary results on deep-mixing manifestations in stellar variability are tested using our improved method and extended data set. We measure the timescales τ {sub m} of the stochastic change in the spectral power of rotational harmonics with numbers m ≤ 3 in the light curves of 1361 main-sequence stars from the Kepler mission archive. We find that the gradient [log(τ{sub 2})−log(τ{sub 1})]/[log(2)−log(1)] has a histogram maximum at −2/3, demonstrating agreement with Kolmogorov’s theory of turbulence and therefore confirming the manifestation of deep mixing. The squared amplitudes of the first and second rotational harmonics, corrected for integral photometry distortion, also show a quasi-Kolmogorov character with spectral index ≈−5/3. Moreover, the reduction of τ {sub 1} and τ {sub 2} to the timescales τ {sub lam1} and τ {sub lam2} of laminar convection in the deep stellar layers reveals the proximity of both τ {sub lam1} and τ {sub lam2} to the turnover time τ {sub MLT} of standard mixing length theory. Considering this result, we use the obtained stellar variability timescales instead of τ {sub MLT} in our analysis of the relation between stellar activity and the Rossby number P/τ {sub MLT}. Comparison of our diagrams with previous results and theoretical expectations shows that best-fit correspondence is achieved for τ {sub lam1}, which can therefore be used as an analog of τ {sub MLT}. This means that the laminar component (giant cells) of stellar turbulent convection indeed plays an important role in the physics of stars. Additionally, we estimate the diffusivity of magnetic elements in stellar photospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

GALACTIC COSMIC-RAY INTENSITY MODULATION BY COROTATING INTERACTION REGION STREAM INTERFACES AT 1 au

We present a new model that couples galactic cosmic-ray (GCR) propagation with magnetic turbulence transport and the MHD background evolution in the heliosphere. The model is applied to the problem of the formation of corotating interaction regions (CIRs) during the last solar minimum from the period between 2007 and 2009. The numerical model simultaneously calculates the large-scale supersonic solar wind properties and its small-scale turbulent content from 0.3 au to the termination shock. Cosmic rays are then transported through the background, and thus computed, with diffusion coefficients derived from the solar wind turbulent properties, using a stochastic Parker approach. Our results demonstrate that GCR variations depend on the ratio of diffusion coefficients in the fast and slow solar winds. Stream interfaces inside the CIRs always lead to depressions of the GCR intensity. On the other hand, heliospheric current sheet (HCS) crossings do not appreciably affect GCR intensities in the model, which is consistent with the two observations under quiet solar wind conditions. Therefore, variations in diffusion coefficients associated with CIR stream interfaces are more important for GCR propagation than the drift effects of the HCS during a negative solar minimum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SIMULATING THE TIMESCALE-DEPENDENT COLOR VARIATION IN QUASARS WITH A REVISED INHOMOGENEOUS DISK MODEL

The UV–optical variability of active galactic nuclei and quasars is useful for understanding the physics of the accretion disk and is gradually being attributed to stochastic fluctuations over the accretion disk. Quasars generally appear bluer when they brighten in the UV–optical bands; the nature of this phenomenon remains controversial. Recently, Sun et al. discovered that the color variation of quasars is timescale-dependent, in the way that faster variations are even bluer than longer term ones. While this discovery can directly rule out models that simply attribute the color variation to contamination from the host galaxies, or to changes in the global accretion rates, it favors the stochastic disk fluctuation model as fluctuations in the inner-most hotter disk could dominate the short-term variations. In this work, we show that a revised inhomogeneous disk model, where the characteristic timescales of thermal fluctuations in the disk are radius-dependent (i.e., τ ∼ r ; based on that originally proposed by Dexter and Agol), can reproduce well a timescale-dependent color variation pattern, similar to the observed one and unaffected by the uneven sampling and photometric error. This demonstrates that one may statistically use variation emission at different timescales to spatially resolve the accretion disk in quasars, thus opening a new window with which to probe and test the accretion disk physics in the era of time domain astronomy. Caveats of the current model, which ought to be addressed in future simulations, are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A NUMERICAL SIMULATION OF COSMIC RAY MODULATION NEAR THE HELIOPAUSE. II. SOME PHYSICAL INSIGHTS

Cosmic ray (CR) transport near the heliopause (HP) is studied using a hybrid transport model, with the parameters constrained by observations from the Voyager 1 spacecraft. We simulate the CR radial flux along different directions in the heliosphere. There is no well-defined thin layer between the solar wind region and the interstellar region along the tail and polar directions of the heliosphere. By analyzing the radial flux curve along the direction of Voyager 2 , together with its trajectory information, the crossing time of the HP by Voyager 2 is predicted to be in 2017.14. We simulate the CR radial flux for different energy values along the direction of Voyager 1 . We find that there is only a modest modulation region of about 10 au wide beyond the HP, so that Voyager 1 observing the Local Interstellar Spectra is justified in numerical modeling. We analyze the heliospheric exit information of pseudo-particles in our stochastic numerical (time-backward) method, conjecturing that they represent the behavior of CR particles, and we find that pseudo-particles that have been traced from the nose region exit in the tail region. This implies that many CR particles diffuse directly from the heliospheric tail region to the nose region near the HP. In addition, when pseudo-particles were traced from the Local Interstellar Medium (LISM), it is found that their exit location (entrance for real particles) from the simulation domain is along the prescribed Interstellar Magnetic Field direction. This indicates that parallel diffusion dominates CR particle transport in the LISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

COMBINING DIFFUSIVE SHOCK ACCELERATION WITH ACCELERATION BY CONTRACTING AND RECONNECTING SMALL-SCALE FLUX ROPES AT HELIOSPHERIC SHOCKS

Computational and observational evidence is accruing that heliospheric shocks, as emitters of vorticity, can produce downstream magnetic flux ropes and filaments. This led Zank et al. to investigate a new paradigm whereby energetic particle acceleration near shocks is a combination of diffusive shock acceleration (DSA) with downstream acceleration by many small-scale contracting and reconnecting (merging) flux ropes. Using a model where flux-rope acceleration involves a first-order Fermi mechanism due to the mean compression of numerous contracting flux ropes, Zank et al. provide theoretical support for observations that power-law spectra of energetic particles downstream of heliospheric shocks can be harder than predicted by DSA theory and that energetic particle intensities should peak behind shocks instead of at shocks as predicted by DSA theory. In this paper, a more extended formalism of kinetic transport theory developed by le Roux et al. is used to further explore this paradigm. We describe how second-order Fermi acceleration, related to the variance in the electromagnetic fields produced by downstream small-scale flux-rope dynamics, modifies the standard DSA model. The results show that (i) this approach can qualitatively reproduce observations of particle intensities peaking behind the shock, thus providing further support for the new paradigm, and (ii) stochastic acceleration by compressible flux ropes tends to be more efficient than incompressible flux ropes behind shocks in modifying the DSA spectrum of energetic particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THERMAL NON-EQUILIBRIUM REVEALED BY PERIODIC PULSES OF RANDOM AMPLITUDES IN SOLAR CORONAL LOOPS

We recently detected variations in extreme ultraviolet intensity in coronal loops repeating with periods of several hours. Models of loops including stratified and quasi-steady heating predict the development of a state of thermal non-equilibrium (TNE): cycles of evaporative upflows at the footpoints followed by falling condensations at the apex. Based on Fourier and wavelet analysis, we demonstrate that the observed periodic signals are indeed not signatures of vibrational modes. Instead, superimposed on the power law expected from the stochastic background emission, the power spectra of the time series exhibit the discrete harmonics and continua expected from periodic trains of pulses of random amplitudes. These characteristics reinforce our earlier interpretation of these pulsations as being aborted TNE cycles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Constraining warm dark matter mass with cosmic reionization and gravitational waves

We constrain the warm dark matter (WDM) particle mass with observations of cosmic reionization and CMB optical depth. We suggest that the gravitational waves (GWs) from stellar-mass black holes (BHs) could give a further constraint on WDM particle mass for future observations. The star formation rates (SFRs) of Population I/II (Pop I/II) and Population III (Pop III) stars are also derived. If the metallicity of the universe is enriched beyond the critical value of Z{sub crit}=10{sup −3.5} Z{sub ⊙}, the star formation shifts from Pop III to Pop I/II stars. Our results show that the SFRs are quite dependent on the WDM particle mass, especially at high redshifts. Combined with the reionization history and CMB optical depth derived from the recent Planck mission, we find that the current data require the WDM particle mass to be in a narrow range of 1 keV≲m{sub x}≲3 keV. Furthermore, we suggest that the stochastic gravitational wave background (SGWB) produced by stellar BHs could give a further constraint on the WDM particle mass for future observations. For m{sub x}=3 keV, with Salpeter (Chabrier) initial mass function (IMF), the SGWB from Pop I/II BHs has a peak amplitude of Ω{sub GW}≈2.8×10{sup −9} (5.0×10{sup −9}) at f=316Hz, while the GW radiation at f<10 Hz is seriously suppressed. For m{sub x}=1 keV, the SGWB peak amplitude is the same as that for m{sub x}=1 keV, but a little lower at low frequencies. Therefore, it is hard to constrain the WDM particle mass by the SGWB from Pop I/II BHs. To assess the detectability of the GW signal, we also calculate the signal-to-noise ratios (S/N), which are S/N=37.7 (66.5) and 27 (47.7) for m{sub x}=3 keV and m{sub x}=1 keV for the Einstein Telescope with Salpeter (Chabrier) IMF, respectively. The SGWB from Pop III BHs is very dependent on the WDM particle mass, the GW strength could be an order of magnitude different, and the frequency band could be two times different for m{sub x}=1 keV and m{sub x}=3 keV. Moreover, the SGWB from Pop III BHs with m{sub x}=1 keV could be detected by the Laser Interferometer Space Antenna for one year of observation, but it cannot be detected for those with m{sub x}=3 keV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

EXPONENTIAL DISKS FROM STELLAR SCATTERING. III. STOCHASTIC MODELS

Stellar scattering off irregularities in a galaxy disk has been shown to make an exponential radial profile, but no fundamental reason for this has been suggested. Here, we show that exponentials are mathematically expected from random scattering in a disk when there is a slight inward bias in the scattering probability. Such a bias was present in our previous scattering experiments that formed exponential profiles. Double exponentials can arise when the bias varies with radius. This is a fundamental property of scattering and may explain why piece-wise exponential profiles are ubiquitous in galaxies, even after minor mergers and other disruptive events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE IMPLICIT CONTRIBUTION OF SLAB MODES TO THE PERPENDICULAR DIFFUSION COEFFICIENT OF PARTICLES INTERACTING WITH TWO-COMPONENT TURBULENCE

We explore the transport of energetic particles in two-component turbulence in which the stochastic magnetic field is assumed to be a superposition of slab and two-dimensional modes. It is known that in magnetostatic slab turbulence, the motion of particles across the mean magnetic field is subdiffusive. If a two-dimensional component is added, diffusion is recovered. It was also shown before that in two-component turbulence, the slab modes do not explicitly contribute to the perpendicular diffusion coefficient. In the current paper, the implicit contribution of slab modes is explored and it is shown that this contribution leads to a reduction of the perpendicular diffusion coefficient. This effect improves the agreement between simulations and analytical theory. Furthermore, the obtained results are relevant for investigations of diffusive shock acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A PURE HYDRODYNAMIC INSTABILITY IN SHEAR FLOWS AND ITS APPLICATION TO ASTROPHYSICAL ACCRETION DISKS

We provide a possible resolution for the century-old problem of hydrodynamic shear flows, which are apparently stable in linear analysis but shown to be turbulent in astrophysically observed data and experiments. This mismatch is noticed in a variety of systems, from laboratory to astrophysical flows. There are so many uncountable attempts made so far to resolve this mismatch, beginning with the early work of Kelvin, Rayleigh, and Reynolds toward the end of the nineteenth century. Here we show that the presence of stochastic noise, whose inevitable presence should not be neglected in the stability analysis of shear flows, leads to pure hydrodynamic linear instability therein. This explains the origin of turbulence, which has been observed/interpreted in astrophysical accretion disks, laboratory experiments, and direct numerical simulations. This is, to the best of our knowledge, the first solution to the long-standing problem of hydrodynamic instability of Rayleigh-stable flows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A UNIFIED MODEL OF GRAIN ALIGNMENT: RADIATIVE ALIGNMENT OF INTERSTELLAR GRAINS WITH MAGNETIC INCLUSIONS

The radiative torque (RAT) alignment of interstellar grains with ordinary paramagnetic susceptibilities has been supported by earlier studies. The alignment of such grains depends on the so-called RAT parameter q {sup max}, which is determined by the grain shape. In this paper, we elaborate on our model of RAT alignment for grains with enhanced magnetic susceptibility due to iron inclusions, such that RAT alignment is magnetically enhanced, which we term the MRAT mechanism. Such grains can be aligned with high angular momentum at the so-called high- J attractor points, achieving a high degree of alignment. Using our analytical model of RATs, we derive the critical value of the magnetic relaxation parameter δ {sub m} to produce high- J attractor points as functions of q {sup max} and the anisotropic radiation angle relative to the magnetic field ψ . We find that if about 10% of the total iron abundance present in silicate grains is forming iron clusters, this is sufficient to produce high- J attractor points for all reasonable values of q {sup max}. To calculate the degree of grain alignment, we carry out numerical simulations of MRAT alignment by including stochastic excitations from gas collisions and magnetic fluctuations. We show that large grains can achieve perfect alignment when the high- J attractor point is present, regardless of the values of q {sup max}. Our obtained results pave the way for the physical modeling of polarized thermal dust emission as well as magnetic dipole emission. We also find that millimeter-sized grains in accretion disks may be aligned with the magnetic field if they are incorporated with iron nanoparticles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A common origin for globular clusters and ultra-faint dwarfs in simulations of the first galaxies

In this paper, the first in a series on galaxy formation before reionization, we focus on understanding what determines the size and morphology of stellar objects in the first low-mass galaxies, using parsec-scale cosmological simulations performed with an adaptive mesh hydrodynamics code. Although the dense gas in which stars are formed tends to have a disk structure, stars are found in spheroids with little rotation. Halos with masses between 10{sup 6} M{sub ⊙} and 5×10{sup 8} M{sub ⊙} form stars stochastically, with stellar masses in the range 10{sup 4} M{sub ⊙} to 2×10{sup 6} M{sub ⊙}. We observe, nearly independent of stellar mass, a large range of half-light radii for the stars, from a few parsecs to a few hundred parsecs and surface brightnesses and mass-to-light ratios ranging from those typical of globular clusters to ultra-faint dwarfs. In our simulations, stars form in dense stellar clusters with high gas-to-star conversion efficiencies and rather uniform metallicities. A fraction of these clusters remain bound after the gas is removed by feedback, but others are destroyed, and their stars, which typically have velocity dispersions of 20–40 km s{sup −1}, expand until they become bound by the dark matter halo. We thus speculate that the stars in ultra-faint dwarf galaxies may show kinematic and chemical signatures consistent with their origin in a few distinct stellar clusters. On the other hand, some globular clusters may form at the center of primordial dwarf galaxies and may contain dark matter, perhaps detectable in the outer parts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Turning around along the cosmic web

A bound violation designates a case in which the turnaround radius of a bound object exceeds the upper limit imposed by the spherical collapse model based on the standard ΛCDM paradigm. Given that the turnaround radius of a bound object is a stochastic quantity and that the spherical model overly simplifies the true gravitational collapse, which actually proceeds anisotropically along the cosmic web, the rarity of the occurrence of a bound violation may depend on the web environment. Assuming a Planck cosmology, we numerically construct the bound-zone peculiar velocity profiles along the cosmic web (filaments and sheets) around the isolated groups with virial mass M{sub v}⩾3×10{sup 13} h{sup −1} M{sub ⊙} identified in the Small MultiDark Planck simulations and determine the radial distances at which their peculiar velocities equal the Hubble expansion speed as the turnaround radii of the groups. It is found that although the average turnaround radii of the isolated groups are well below the spherical bound limit on all mass scales, the bound violations are not forbidden for individual groups, and the cosmic web has an effect of reducing the rarity of the occurrence of a bound violation. Explaining that the spherical bound limit on the turnaround radius in fact represents the threshold distance up to which the intervention of the external gravitational field in the bound-zone peculiar velocity profiles around the nonisolated groups stays negligible, we discuss the possibility of using the threshold distance scale to constrain locally the equation of state of dark energy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

ANALYZING METEOROID FLIGHTS USING PARTICLE FILTERS

Fireball observations from camera networks provide position and time information along the trajectory of a meteoroid that is transiting our atmosphere. The complete dynamical state of the meteoroid at each measured time can be estimated using Bayesian filtering techniques. A particle filter is a novel approach to modeling the uncertainty in meteoroid trajectories and incorporates errors in initial parameters, the dynamical model used, and observed position measurements. Unlike other stochastic approaches, a particle filter does not require predefined values for initial conditions or unobservable trajectory parameters. The Bunburra Rockhole fireball, observed by the Australian Desert Fireball Network (DFN) in 2007, is used to determine the effectiveness of a particle filter for use in fireball trajectory modeling. The final mass is determined to be 2.16±1.33 kg with a final velocity of 6030±216 m s{sup −1}, similar to previously calculated values. The full automatability of this approach will allow an unbiased evaluation of all events observed by the DFN and lead to a better understanding of the dynamical state and size frequency distribution of asteroid and cometary debris in the inner solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

An Updated 2017 Astrometric Solution for Betelgeuse

We provide an update for the astrometric solution for the Type II supernova progenitor Betelgeuse using the revised Hipparcos Intermediate Astrometric Data (HIAD) of van Leeuwen, combined with existing VLA and new e-MERLIN and ALMA positions. The 2007 Hipparcos refined abscissa measurements required the addition of so-called Cosmic Noise of 2.4 mas to find an acceptable 5-parameter stochastic solution. We find that a measure of radio Cosmic Noise should also be included for the radio positions because surface inhomogeneities exist at a level significant enough to introduce additional intensity centroid uncertainty. Combining the 2007 HIAD with the proper motions based solely on the radio positions leads to a parallax of π=5.27±0.78 mas (190{sub −25}{sup +33} pc), smaller than the Hipparcos 2007 value of 6.56 ± 0.83 mas (152{sub −17}{sup +22} pc). Furthermore, combining the VLA and new e-MERLIN and ALMA radio positions with the 2007 HIAD, and including radio Cosmic Noise of 2.4 mas, leads to a nominal parallax solution of 4.51 ± 0.80 mas (222{sub −34}{sup +48} pc), which, while only 0.7σ different from the 2008 solution of Harper et al., is 2.6σ different from the solution of van Leeuwen. An accurate and precise parallax for Betelgeuse is always going to be difficult to obtain because it is small compared to the stellar angular diameter (θ=44 mas). We outline an observing strategy utilizing future mm and sub-mm high-spatial resolution interferometry that must be used if substantial improvements in the precision and accuracy of the parallax and distance are to be achieved.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Young Substellar Companion ROXs 12 B: Near-infrared Spectrum, System Architecture, and Spin–Orbit Misalignment

ROXs 12 (2MASS J16262803–2526477) is a young star hosting a directly imaged companion near the deuterium-burning limit. We present a suite of spectroscopic, imaging, and time-series observations to characterize the physical and environmental properties of this system. Moderate-resolution near-infrared spectroscopy of ROXs 12 B from Gemini-North/NIFS and Keck/OSIRIS reveals signatures of low surface gravity including weak alkali absorption lines and a triangular H-band pseudocontinuum shape. No signs of Paβ emission are evident. As a population, however, we find that about half (46% ± 14%) of young (≲15 Myr) companions with masses ≲20 M {sub Jup} possess actively accreting subdisks detected via Paβ line emission, which represents a lower limit on the prevalence of circumplanetary disks in general, as some are expected to be in a quiescent phase of accretion. The bolometric luminosity of the companion and age of the host star (6{sub −2}{sup +4} Myr) imply a mass of 17.5 ± 1.5 M {sub Jup} for ROXs 12 B based on hot-start evolutionary models. We identify a wide (5100 au) tertiary companion to this system, 2MASS J16262774–2527247, that is heavily accreting and exhibits stochastic variability in its K2 light curve. By combining v sin i {sub *} measurements with rotation periods from K2, we constrain the line-of-sight inclinations of ROXs 12 A and 2MASS J16262774–2527247 and find that they are misaligned by 60{sub −11}{sup +7∘}. In addition, the orbital axis of ROXs 12 B is likely misaligned from the spin axis of its host star, ROXs 12 A, suggesting that ROXs 12 B formed akin to fragmenting binary stars or in an equatorial disk that was torqued by the wide stellar tertiary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Chaotic Excitation and Tidal Damping in the GJ 876 System

The M-dwarf GJ 876 is the closest known star to harbor a multi-planetary system. With three outer planets locked in a chaotic Laplace-type resonance and an appreciably eccentric short-period super-Earth, this system represents a unique exposition of extrasolar planetary dynamics. A key question that concerns the long-term evolution of this system, and the fate of close-in planets in general, is how the significant eccentricity of the inner-most planet is maintained against tidal circularization on timescales comparable to the age of the universe. Here, we employ stochastic secular perturbation theory and N-body simulations to show that the orbit of the inner-most planet is shaped by a delicate balance between extrinsic chaotic forcing and tidal dissipation. As such, the planet’s orbital eccentricity represents an indirect measure of its tidal quality factor. Based on the system’s present-day architecture, we estimate that the extrasolar super-Earth GJ 876 d has a tidal Q ∼ 10{sup 4}–10{sup 5}, a value characteristic of solar system gas giants.

79 ASTRONOMY AND ASTROPHYSICS↗