Meridional flow in the solar interior plays an important role in redistributing angular momentum and transporting magnetic flux inside the Sun. Although it has long been recognized that the meridional flow is predominantly poleward at the Sun's surface and in its shallow interior, the location of the equatorward return flow and the meridional flow profile in the deeper interior remain unclear. Using the first 2 yr of continuous helioseismology observations from the Solar Dynamics Observatory/Helioseismic Magnetic Imager, we analyze travel times of acoustic waves that propagate through different depths of the solar interior carrying information about the solar interior dynamics. After removing a systematic center-to-limb effect in the helioseismic measurements and performing inversions for flow speed, we find that the poleward meridional flow of a speed of 15 m s −1 extends in depth from the photosphere to about 0.91 R . An equatorward flow of a speed of 10 m s −1 is found between 0.82 and 0.91 R in the middle of the convection zone. Our analysis also shows evidence of that the meridional flow turns poleward again below 0.82 R , indicating an existence of a second meridional circulation cell below the shallower one. This double-cell meridional circulation profile with an equatorward flow shallower than previously thought suggests a rethinking of how magnetic field is generated and redistributed inside the Sun.

Direct numerical simulations of Rayleigh-Bénard convection in a plane layer with periodic boundary conditions at Rayleigh numbers up to 10(7) show that flow structures can be objectively classified as large or small scale structures because of a gap in spatial spectra. The typical size of the large scale structures does not always vary monotonically as a function of the Rayleigh number but broadly increases with increasing Rayleigh number. A mean flow (whose average over horizontal planes differs from zero) is also excited but is weak in comparison with the large scale structures. The large scale circulation observed in experiments should therefore be a manifestation of the large scale structures identified here.

Numerical simulations of Rayleigh–Bénard convection in a fluid layer heated from below between two rigid horizontal boundaries have been performed for Rayleigh numbers Ra up to 107, Prandtl numbers in the range between 0.7 and 60 and for aspect ratios $\Gamma$ up to 20. Periodic boundary conditions in the horizontal plane have been used. To a considerable extent, the evolution towards turbulent convection at high values of Ra is governed by processes exhibited by instabilities of steady or time periodic forms of convection at lower Rayleigh numbers. With increasing Ra, the properties of convection are increasingly determined by the thermal boundary layers. The role of mean flows which may be symmetric or antisymmetric with respect to the mid-plane of the layer is emphasized.

We show that under certain conditions, subsurface structures in the solar interior can alter the average acoustic power observed at the photosphere above them. By using numerical simulations of wave propagation, we show that this effect is large enough for it to be potentially used for detecting emerging active regions before they appear on the surface. In our simulations, simplified subsurface structures are modeled as regions with enhanced or reduced acoustic wave speed. We investigate the dependence of the acoustic power above a subsurface region on the sign, depth, and strength of the wave-speed perturbation. Observations from the Solar and Heliospheric Observatory/Michelson Doppler Imager (SOHO/MDI) prior and during the emergence of NOAA active region 10488 are used to test the use of acoustic power as a potential precursor of the emergence of magnetic flux.

Turbulent flows preferentially concentrate inertial particles depending on their stopping time or Stokes number, which can lead to significant spatial variations in the particle concentration. Cascade models are one way to describe this process in statistical terms. Here, we use a direct numerical simulation (DNS) dataset of homogeneous, isotropic turbulence to determine probability distribution functions (PDFs) for cascade multipliers, which determine the ratio by which a property is partitioned into sub-volumes as an eddy is envisioned to decay into smaller eddies. We present a technique for correcting effects of small particle numbers in the statistics. We determine multiplier PDFs for particle number, flow dissipation, and enstrophy, all of which are shown to be scale dependent. However, the particle multiplier PDFs collapse when scaled with an appropriately defined local Stokes number. As anticipated from earlier works, dissipation and enstrophy multiplier PDFs reach an asymptote for sufficiently small spatial scales. From the DNS measurements, we derive a cascade model that is used it to make predictions for the radial distribution function (RDF) for arbitrarily high Reynolds numbers, Re, finding good agreement with the asymptotic, infinite Re inertial range theory of Zaichik & Alipchenkov [New Journal of Physics 11, 103018 (2009)]. We discuss implications of these results for the statistical modeling of the turbulent clustering process in the inertial range for high Reynolds numbers inaccessible to numerical simulations.