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,477 records · Page 82

Mergers in galaxy groups. I. Structure and properties of elliptical remnants

We present collisionless simulations of dry mergers in groups of 3 to 25 galaxies to test the hypothesis that elliptical galaxies form at the centers of such groups. Mock observations of the central remnants confirm their similarity to ellipticals, despite having no dissipational component. We vary the profile of the original spiral's bulge and find that ellipticals formed from spirals with exponential bulges have too low Sersic indices. Mergers of spirals with de Vaucouleurs (classical) bulges produce remnants with larger Sersic indices correlated with luminosity, as with Sloan Digital Sky Survey ellipticals. Exponential bulge mergers are better fits to faint ellipticals, whereas classical bulge mergers better match luminous ellipticals. Similarly, luminous ellipticals are better reproduced by remnants undergoing many (>5) mergers, and fainter ellipticals by those with fewer mergers. The remnants follow tight size-luminosity and velocity dispersion-luminosity (Faber-Jackson) relations (<0.12 dex scatter), demonstrating that stochastic merging can produce tight scaling relations if the merging galaxies also follow tight scaling relations. The slopes of the size-luminosity and Faber-Jackson relations are close to observations but slightly shallower in the former case. Both relations' intercepts are offset—remnants are too large but have too low dispersions at fixed luminosity. Some remnants show substantial (v/σ > 0.1) rotational support, although most are slow rotators and few are very fast rotators (v/σ > 0.5). These findings contrast with previous studies concluding that dissipation is necessary to produce ellipticals from binary mergers of spirals. Multiple, mostly minor and dry mergers can produce bright ellipticals, whereas significant dissipation could be required to produce faint, rapidly rotating ellipticals.

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic stellar cluster initial mass functions: Models and impact on integrated cluster parameter determination

Stellar clusters are regularly used to study the evolution of their host galaxy. Except for a few nearby galaxies, these studies rely on the interpretation of integrated cluster properties, especially integrated photometry observed using multiple filters (i.e., the spectral energy distribution, SED). To allow interpretation of such observations, we present a large set of GALEV cluster models using the realistic approach of adopting stochastically sampled stellar initial mass functions. We provide models for a wide range of cluster masses (10{sup 3}-2 × 10{sup 5} M {sub ☉}), metallicities (–2.3 ≤ [Fe/H] ≤ +0.18 dex), foreground extinction, and 184 regularly used filters. We analyze various sets of stochastic cluster SEDs by fitting them with non-stochastic models, which is the procedure commonly used in this field. We identify caveats and quantify the fitting uncertainties associated with this standard procedure. We show that this can yield highly unreliable fitting results, especially for low-mass clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Global hydromagnetic simulations of a planet embedded in a dead zone: Gap opening, gas accretion, and formation of a protoplanetary jet

We present global hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations with mesh refinement of accreting planets embedded in protoplanetary disks (PPDs). The magnetized disk includes Ohmic resistivity that depends on the overlying mass column, leading to turbulent surface layers and a dead zone near the midplane. The main results are: (1) the accretion flow in the Hill sphere is intrinsically three-dimensional for HD and MHD models. Net inflow toward the planet is dominated by high-latitude flows. A circumplanetary disk (CPD) forms. Its midplane flows outward in a pattern whose details differ between models. (2) The opening of a gap magnetically couples and ignites the dead zone near the planet, leading to stochastic accretion, a quasi-turbulent flow in the Hill sphere, and a CPD whose structure displays high levels of variability. (3) Advection of magnetized gas onto the rotating CPD generates helical fields that launch magnetocentrifugally driven outflows. During one specific epoch, a highly collimated, one-sided jet is observed. (4) The CPD's surface density is ∼30 g cm{sup −2}, small enough for significant ionization and turbulence to develop. (5) The accretion rate onto the planet in the MHD simulation reaches a steady value 8 × 10{sup –3} M {sub ⊕} yr{sup –1} and is similar in the viscous HD runs. Our results suggest that gas accretion onto a forming giant planet within a magnetized PPD with a dead zone allows rapid growth from Saturnian to Jovian masses. As well as being relevant for giant planet formation, these results have important implications for the formation of regular satellites around gas giant planets.

79 ASTRONOMY AND ASTROPHYSICS↗

The era of star formation in galaxy clusters

We analyze the star formation properties of 16 infrared-selected, spectroscopically confirmed galaxy clusters at 1 < z < 1.5 from the Spitzer/IRAC Shallow Cluster Survey (ISCS). We present new spectroscopic confirmation for six of these high-redshift clusters, five of which are at z > 1.35. Using infrared luminosities measured with deep Spitzer/Multiband Imaging Photometer for Spitzer observations at 24 μm, along with robust optical + IRAC photometric redshifts and spectral-energy-distribution-fitted stellar masses, we present the dust-obscured star-forming fractions, star formation rates, and specific star formation rates in these clusters as functions of redshift and projected clustercentric radius. We find that z ∼ 1.4 represents a transition redshift for the ISCS sample, with clear evidence of an unquenched era of cluster star formation at earlier times. Beyond this redshift, the fraction of star-forming cluster members increases monotonically toward the cluster centers. Indeed, the specific star formation rate in the cores of these distant clusters is consistent with field values at similar redshifts, indicating that at z > 1.4 environment-dependent quenching had not yet been established in ISCS clusters. By combining these observations with complementary studies showing a rapid increase in the active galactic nucleus (AGN) fraction, a stochastic star formation history, and a major merging episode at the same epoch in this cluster sample, we suggest that the starburst activity is likely merger-driven and that the subsequent quenching is due to feedback from merger-fueled AGNs. The totality of the evidence suggests we are witnessing the final quenching period that brings an end to the era of star formation in galaxy clusters and initiates the era of passive evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

The range of variation of the mass of the most massive star in stellar clusters derived from 35 million Monte Carlo simulations

A growing fraction of simple stellar population models, in an aim to create more realistic simulations capable of including stochastic variation in their outputs, begin their simulations with a distribution of discrete stars following a power-law function of masses. Careful attention is needed to create a correctly sampled initial mass function (IMF), and here we provide a solid mathematical method, called MASSCLEAN IMF Sampling, for doing so. We use our method to perform 10 million MASSCLEAN Monte Carlo stellar cluster simulations to determine the most massive star in a mass distribution as a function of the total mass of the cluster. We find that a maximum mass range is predicted, not a single maximum mass. This range is (1) dependent on the total mass of the cluster and (2) independent of an upper stellar mass limit, M{sub limit} , for unsaturated clusters and emerges naturally from our IMF sampling method. We then turn our analysis around, starting with our new database of 25 million simulated clusters, to constrain the highest mass star from the observed integrated colors of a sample of 40 low-mass Large Magellanic Cloud stellar clusters of known age and mass. Finally, we present an analytical description of the maximum mass range of the most massive star as a function of the cluster's total mass and present a new M{sub max} -M{sub cluster} relation.

79 ASTRONOMY AND ASTROPHYSICS↗

Ages of young star clusters, massive blue stragglers, and the upper mass limit of stars: Analyzing age-dependent stellar mass functions

Massive stars rapidly change their masses through strong stellar winds and mass transfer in binary systems. The latter aspect is important for populations of massive stars as more than 70% of all O stars are expected to interact with a binary companion during their lifetime. We show that such mass changes leave characteristic signatures in stellar mass functions of young star clusters that can be used to infer their ages and to identify products of binary evolution. We model the observed present-day mass functions of the young Galactic Arches and Quintuplet star clusters using our rapid binary evolution code. We find that the shaping of the mass function by stellar wind mass loss allows us to determine the cluster ages as 3.5 ± 0.7 Myr and 4.8 ± 1.1 Myr, respectively. Exploiting the effects of binary mass exchange on the cluster mass function, we find that the most massive stars in both clusters are rejuvenated products of binary mass transfer, i.e., the massive counterpart of classical blue straggler stars. This resolves the problem of an apparent age spread among the most luminous stars exceeding the expected duration of star formation in these clusters. We perform Monte Carlo simulations to probe stochastic sampling, which support the idea of the most massive stars being rejuvenated binary products. We find that the most massive star is expected to be a binary product after 1.0 ± 0.7 Myr in Arches and after 1.7 ± 1.0 Myr in Quintuplet. Today, the most massive 9 ± 3 stars in Arches and 8 ± 3 in Quintuplet are expected to be such objects. Our findings have strong implications for the stellar upper mass limit and solve the discrepancy between the claimed 150 M {sub ☉} limit and observations of four stars with initial masses of 165-320 M {sub ☉} in R136 and of supernova 2007bi, which is thought to be a pair-instability supernova from an initial 250 M {sub ☉} star. Using the stellar population of R136, we revise the upper mass limit to values in the range 200-500 M {sub ☉}.

79 ASTRONOMY AND ASTROPHYSICS↗

Stellar granulation as the source of high-frequency flicker in Kepler light curves

A large fraction of cool, low-mass stars exhibit brightness fluctuations that arise from a combination of convective granulation, acoustic oscillations, magnetic activity, and stellar rotation. Much of the short-timescale variability takes the form of stochastic noise, whose presence may limit the progress of extrasolar planet detection and characterization. In order to lay the groundwork for extracting useful information from these quasi-random signals, we focus on the origin of the granulation-driven component of the variability. We apply existing theoretical scaling relations to predict the star-integrated variability amplitudes for 508 stars with photometric light curves measured by the Kepler mission. We also derive an empirical correction factor that aims to account for the suppression of convection in F-dwarf stars with magnetic activity and shallow convection zones. So that we can make predictions of specific observational quantities, we performed Monte Carlo simulations of granulation light curves using a Lorentzian power spectrum. These simulations allowed us to reproduce the so-called flicker floor (i.e., a lower bound in the relationship between the full light-curve range and power in short-timescale fluctuations) that was found in the Kepler data. The Monte Carlo model also enabled us to convert the modeled fluctuation variance into a flicker amplitude directly comparable with observations. When the magnetic suppression factor described above is applied, the model reproduces the observed correlation between stellar surface gravity and flicker amplitude. Observationally validated models like these provide new and complementary evidence for a possible impact of magnetic activity on the properties of near-surface convection.

79 ASTRONOMY AND ASTROPHYSICS↗

The statistical nature of the brightest group galaxies

We examine the statistical properties of the brightest group galaxies (BGGs) by using a complete spectroscopic sample of groups/clusters of galaxies selected from the Data Release 7 of the Sloan Digital Sky Survey. We test whether BGGs and other bright members of groups are consistent with an ordered population among the total population of group galaxies. We find that the luminosity distributions of BGGs do not follow the predictions from the order statistics (OS). The average luminosities of BGGs are systematically brighter than OS predictions. On the other hand, by properly taking into account the brightening effect of the BGGs, the luminosity distributions of the second brightest galaxies are in excellent agreement with the expectations of OS. The brightening of BGGs relative to the OS expectation is consistent with a scenario that the BGGs on average have overgrown about 20% masses relative to the other member galaxies. The growth (ΔM) is not stochastic but correlated with the magnitude gap (G {sub 1,} {sub 2}) between the brightest and the second brightest galaxy. The growth (ΔM) is larger for the groups having more prominent BGGs (larger G {sub 1,} {sub 2}) and averagely contributes about 30% of the final G {sub 1,} {sub 2} of the groups of galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

A 33 yr constancy of the X-ray coronae of AR Lac and eclipse diagnosis of scale height

Extensive X-ray and EUV photometric observations of the eclipsing RS CVn system AR Lac were obtained over the years 1997-2013 with the Chandra X-Ray Observatory Extreme-Ultraviolet Explorer (EUVE). During primary eclipse, High Resolution Camera count rates decrease by ∼40%. A similar minimum is seen during one primary eclipse observed by EUVE but not in others owing to intrinsic source variability. Little evidence for secondary eclipses is present in either the X-ray or EUV data, reminiscent of earlier X-ray and EUV observations. Primary eclipses allow us to estimate the extent of a spherically symmetric corona on the primary G star of about 1.3 R {sub ☉}, or 0.86 R {sub *}, and indicate that the G star is likely brighter than the K component by a factor of 2-5. Brightness changes not attributable to eclipses appear to be dominated by stochastic variability and are generally non-repeating. X-ray and EUV light curves cannot therefore be reliably used to reconstruct the spatial distribution of emission assuming that only eclipses and rotational modulation are at work. Moderate flaring is observed, where count rates increase by up to a factor of three above quiescence. Combined with older ASCA, Einstein, EXOSAT, ROSAT, and BeppoSAX observations, the data show that the level of quiescent coronal emission at X-ray wavelengths has remained remarkably constant over 33 yr, with no sign of variation due to magnetic cycles. Variations in base level X-ray emission seen by Chandra over 13 yr are only ∼10%, while variations back to pioneering Einstein observations in 1980 amount to a maximum of 45% and more typically about 15%.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetohydrodynamic simulations of global accretion disks with vertical magnetic fields

We report results of three-dimensional magnetohydrodynamical (MHD) simulations of global accretion disks threaded with weak vertical magnetic fields. We perform the simulations in the spherical coordinates with different temperature profiles and accordingly different rotation profiles. In the cases with a spatially constant temperature, because the rotation frequency is vertically constant in the equilibrium condition, general properties of the turbulence excited by magnetorotational instability are quantitatively similar to those obtained in local shearing box simulations. On the other hand, in the cases with a radially variable temperature profile, the vertical differential rotation, which is inevitable in the equilibrium condition, winds up the magnetic field lines in addition to the usual radial differential rotation. As a result, the coherent wound magnetic fields contribute to the Maxwell stress in the surface regions. We obtain nondimensional density and velocity fluctuations ∼0.1-0.2 at the midplane. The azimuthal power spectra of the magnetic fields show shallower slopes, ∼m {sup 0} – m {sup –1}, than those of velocity and density. The Poynting flux associated with the MHD turbulence drives intermittent and structured disk winds as well as sound-like waves toward the midplane. The mass accretion mainly occurs near the surfaces, and the gas near the midplane slowly moves outward in the time domain of the present simulations. The vertical magnetic fields are also dragged inward in the surface regions, while they stochastically move outward and inward around the midplane. We also discuss an observational implication of induced spiral structure in the simulated turbulent disks.

79 ASTRONOMY AND ASTROPHYSICS↗

An upper bound from helioseismology on the stochastic background of gravitational waves

The universe is expected to be permeated by a stochastic background of gravitational radiation of astrophysical and cosmological origin. This background is capable of exciting oscillations in solar-like stars. Here we show that solar-like oscillators can be employed as giant hydrodynamical detectors for such a background in the μHz to mHz frequency range, which has remained essentially unexplored until today. We demonstrate this approach by using high-precision radial velocity data for the Sun to constrain the normalized energy density of the stochastic gravitational-wave background around 0.11 mHz. These results open up the possibility for asteroseismic missions like CoRoT and Kepler to probe fundamental physics.

79 ASTRONOMY AND ASTROPHYSICS↗

Verification of the helioseismology travel-time measurement technique and the inversion procedure for sound speed using artificial data

We performed three-dimensional numerical simulations of the solar surface acoustic wave field for the quiet Sun and for three models with different localized sound-speed perturbations in the interior with deep, shallow, and two-layer structures. We used the simulated data generated by two solar acoustics codes that employ the same standard solar model as a background model, but utilize different integration techniques and different models of stochastic wave excitation. Acoustic travel times were measured using a time-distance helioseismology technique, and compared with predictions from ray theory frequently used for helioseismic travel-time inversions. It is found that the measured travel-time shifts agree well with the helioseismic theory for sound-speed perturbations, and for the measurement procedure with and without phase-speed filtering of the oscillation signals. This testing verifies the whole measuring-filtering-inversion procedure for static sound-speed anomalies with small amplitude inside the Sun outside regions of strong magnetic field. It is shown that the phase-speed filtering, frequently used to extract specific wave packets and improve the signal-to-noise ratio, does not introduce significant systematic errors. Results of the sound-speed inversion procedure show good agreement with the perturbation models in all cases. Due to its smoothing nature, the inversion procedure may overestimate sound-speed variations in regions with sharp gradients of the sound-speed profile.

79 ASTRONOMY AND ASTROPHYSICS↗

The panchromatic Hubble Andromeda Treasury. V. Ages and masses of the year 1 stellar clusters

We present ages and masses for 601 star clusters in M31 from the analysis of the six filter integrated light measurements from near-ultraviolet to near-infrared wavelengths, made as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We derive the ages and masses using a probabilistic technique, which accounts for the effects of stochastic sampling of the stellar initial mass function. Tests on synthetic data show that this method, in conjunction with the exquisite sensitivity of the PHAT observations and their broad wavelength baseline, provides robust age and mass recovery for clusters ranging from ∼10{sup 2} to 2 × 10{sup 6} M {sub ☉}. We find that the cluster age distribution is consistent with being uniform over the past 100 Myr, which suggests a weak effect of cluster disruption within M31. The age distribution of older (>100 Myr) clusters falls toward old ages, consistent with a power-law decline of index –1, likely from a combination of fading and disruption of the clusters. We find that the mass distribution of the whole sample can be well described by a single power law with a spectral index of –1.9 ± 0.1 over the range of 10{sup 3}-3 × 10{sup 5} M {sub ☉}. However, if we subdivide the sample by galactocentric radius, we find that the age distributions remain unchanged. However, the mass spectral index varies significantly, showing best-fit values between –2.2 and –1.8, with the shallower slope in the highest star formation intensity regions. We explore the robustness of our study to potential systematics and conclude that the cluster mass function may vary with respect to environment.

79 ASTRONOMY AND ASTROPHYSICS↗

Disk-loss and disk-renewal phases in classical Be stars. II. Contrasting with stable and variable disks

Recent observational and theoretical studies of classical Be stars have established the utility of polarization color diagrams (PCDs) in helping to constrain the time-dependent mass decretion rates of these systems. We expand on our pilot observational study of this phenomenon, and report the detailed analysis of a long-term (1989-2004) spectropolarimetric survey of nine additional classical Be stars, including systems exhibiting evidence of partial disk-loss/disk-growth episodes as well as systems exhibiting long-term stable disks. After carefully characterizing and removing the interstellar polarization along the line of sight to each of these targets, we analyze their intrinsic polarization behavior. We find that many steady-state Be disks pause at the top of the PCD, as predicted by theory. We also observe sharp declines in the Balmer jump polarization for later spectral type, near edge-on steady-state disks, again as recently predicted by theory, likely caused when the base density of the disk is very high, and the outer region of the edge-on disk starts to self absorb a significant number of Balmer jump photons. The intrinsic V-band polarization and polarization position angle of γ Cas exhibits variations that seem to phase with the orbital period of a known one-armed density structure in this disk, similar to the theoretical predictions of Halonen and Jones. We also observe stochastic jumps in the intrinsic polarization across the Balmer jump of several known Be+sdO systems, and speculate that the thermal inflation of part of the outer region of these disks could be responsible for producing this observational phenomenon. Finally, we estimate the base densities of this sample of stars to be between ≈8 × 10{sup –11} and ≈4 × 10{sup –12} g cm{sup –3} during quasi steady state periods given there maximum observed polarization.

79 ASTRONOMY AND ASTROPHYSICS↗

A comparison of two- and three-dimensional neutrino-hydrodynamics simulations of core-collapse supernovae

We present numerical results on two- (2D) and three-dimensional (3D) hydrodynamic core-collapse simulations of an 11.2 M {sub ☉} star. By changing numerical resolutions and seed perturbations systematically, we study how the postbounce dynamics are different in 2D and 3D. The calculations were performed with an energy-dependent treatment of the neutrino transport based on the isotropic diffusion source approximation scheme, which we have updated to achieve a very high computational efficiency. All of the computed models in this work, including nine 3D models and fifteen 2D models, exhibit the revival of the stalled bounce shock, leading to the possibility of explosion. All of them are driven by the neutrino-heating mechanism, which is fostered by neutrino-driven convection and the standing-accretion-shock instability. Reflecting the stochastic nature of multi-dimensional (multi-D) neutrino-driven explosions, the blast morphology changes from model to model. However, we find that the final fate of the multi-D models, whether an explosion is obtained or not, is little affected by the explosion stochasticity. In agreement with some previous studies, higher numerical resolutions lead to slower onset of the shock revival in both 2D and 3D. Based on the self-consistent supernova models leading to the possibility of explosions, our results systematically show that the revived shock expands more energetically in 2D than in 3D.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmic ray transport in heliospheric magnetic structures. I. Modeling background solar wind using the CRONOS magnetohydrodynamic code

The transport of energetic particles such as cosmic rays is governed by the properties of the plasma being traversed. While these properties are rather poorly known for galactic and interstellar plasmas due to the lack of in situ measurements, the heliospheric plasma environment has been probed by spacecraft for decades and provides a unique opportunity for testing transport theories. Of particular interest for the three-dimensional (3D) heliospheric transport of energetic particles are structures such as corotating interaction regions, which, due to strongly enhanced magnetic field strengths, turbulence, and associated shocks, can act as diffusion barriers on the one hand, but also as accelerators of low energy CRs on the other hand as well. In a two-fold series of papers, we investigate these effects by modeling inner-heliospheric solar wind conditions with a numerical magnetohydrodynamic (MHD) setup (this paper), which will serve as an input to a transport code employing a stochastic differential equation approach (second paper). In this first paper, we present results from 3D MHD simulations with our code CRONOS: for validation purposes we use analytic boundary conditions and compare with similar work by Pizzo. For a more realistic modeling of solar wind conditions, boundary conditions derived from synoptic magnetograms via the Wang-Sheeley-Arge (WSA) model are utilized, where the potential field modeling is performed with a finite-difference approach in contrast to the traditional spherical harmonics expansion often utilized in the WSA model. Our results are validated by comparing with multi-spacecraft data for ecliptical (STEREO-A/B) and out-of-ecliptic (Ulysses) regions.

79 ASTRONOMY AND ASTROPHYSICS↗

A stochastically forced time delay solar dynamo model: Self-consistent recovery from a maunder-like grand minimum necessitates a mean-field alpha effect

Fluctuations in the Sun's magnetic activity, including episodes of grand minima such as the Maunder minimum have important consequences for space and planetary environments. However, the underlying dynamics of such extreme fluctuations remain ill-understood. Here, we use a novel mathematical model based on stochastically forced, non-linear delay differential equations to study solar cycle fluctuations in which time delays capture the physics of magnetic flux transport between spatially segregated dynamo source regions in the solar interior. Using this model, we explicitly demonstrate that the Babcock-Leighton poloidal field source based on dispersal of tilted bipolar sunspot flux, alone, cannot recover the sunspot cycle from a grand minimum. We find that an additional poloidal field source effective on weak fields—e.g., the mean-field α effect driven by helical turbulence—is necessary for self-consistent recovery of the sunspot cycle from grand minima episodes.

79 ASTRONOMY AND ASTROPHYSICS↗