Turbulence properties in the solar convection envelope： properties of the overshooting regions

We apply the turbulent convection model (TCM) to investigate properties of tur-bulence in the solar convective envelope, especially in overshooting regions. The results show TCM gives negative turbulent heat flux uγ′T′in overshooting regions, which is sim-ilar to other nonlocal turbulent convection theories. The turbulent temperature fluctuation T′T′shows peaks in overshooting regions. Most important, we find that the downward overshooting region below the base of the solar convection zone is a thin cellular layer filled with roll-shaped convective cells. The overshooting length for the temperature gradi-ent is much shorter than that for element mixing because turbulent heat flux of downward and upward moving convective cells counteract each other in this cellular overshooting region. Comparing the models' sound speed with observations, we find that raking the convective overshooting into account helps to improve the sound speed profile of our nonlocal solar models. Comparing the p-mode oscillation frequencies with observations,we validated that increasing the diffusion parameters and decreasing the dissipation pa-rameters of TCM make the p-mode oscillation frequencies of the solar model be in betteragreement with observations.     

Intermediate degree p-mode frequency splittings near solar maximum

We report on observations of global solar Ca K-line intensity oscillations taken in May 1991 from Mees Solar Observatory, Hawaii. We measurep-mode frequency splittings for modes of spherical harmonic degrees between 20 and 129 averaged over the radial order of the modes. Our measurement of the antisymmetric component of the splittings is comparable with previous measurements and thus indicates a decrease in the latitudinal differential rotation with depth into the convection zone and the upper radiative zone. We find evidence for a 1% variation in the rotation rate of the upper convection zone roughly in phase with the solar activity cycle. Our measurement of the symmetric component of the splittings is of the same order as was reported from the previous solar maximum and is an order of magnitude larger than has been measured near solar minimum. From the degree dependence of the symmetric component of the splittings, we find evidence for an aspherical fractional sound speed perturbation located at a depth of 0.85 ± 0.05 solar radii. This perturbation has a magnitude ofc/c +9 × 10^–4 at the equator relative to the poles. Additionally, there is evidence for a near-surface aspherical sound speed perturbation of smaller magnitudec/c +4 × 10^–4 at the equator relative to the poles. If an intense global magnetic field were the dominant source of the observed symmetric component of the splittings, instead of latitudinal gradients in the sound speed, then global fields of order 10⁵ G would be required in the convection zone.     

The structure of the solar convective overshooting zone

Using a non-local theory of convection, we calculated the structure of the solar convection zone, paying special attention to the detailed structure of the lower overshooting zone. Our results show that an extended transition zone exists near the bottom of the convection zone, where the temperature gradient turns smoothly from adiabatic in the convection zone to radiative in solar interior. A super-radiative temperature region is found in the overshooting zone under the solar convection zone, where     ,     ,     and     . The extension of the super-radiative region (defined by     l is about 0.63  H _P (0.053 R_⊙). A careful comparison of the distribution of adiabatic sound speed and density with the local one is carried out. It is found, strikingly, that the distribution of adiabatic sound speed and density of our model is roughly consistent with the results of reversion from solar oscillation observations.     

Differential rotation and meridional flow of Arcturus

The spectroscopic variability of Arcturus hints at cyclic activity cycle and differential rotation. This could provide a test of current theoretical models of solar and stellar dynamos. To examine the applicability of current models of the flux transport dynamo to Arcturus, we compute a mean‐field model for its internal rotation, meridional flow, and convective heat transport in the convective envelope. We then compare the conditions for dynamo action with those on the Sun. We find solar‐type surface rotation with about 1/10th of the shear found on the solar surface. The rotation rate increases monotonically with depth at all latitudes throughout the whole convection zone. In the lower part of the convection zone the horizontal shear vanishes and there is a strong radial gradient. The surface meridional flow has maximum speed of 170 m/s and is directed towards the equator at high and towards the poles at low latitudes. Turbulent magnetic diffusivity is of the order 10¹⁵–10¹⁶ cm²/s. The conditions on Arcturus are not favorable for a circulation‐dominated dynamo (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Solar dynamo models with α ‐effect and turbulent pumping from local 3D convection calculations

Results from kinematic solar dynamo models employing α ‐effect and turbulent pumping from local convection calculations are presented. We estimate the magnitude of these effects to be around 2–3 m s^–1, having scaled the local quantities with the convective velocity at the bottom of the convection zone from a solar mixing‐length model. Rotation profile of the Sun as obtained from helioseismology is applied in the models; we also investigate the effects of the observed surface shear layer on the dynamo solutions. With these choices of the small‐ and large‐scale velocity fields, we obtain estimate of the ratio of the two induction effects, C _α /C _Ω ≈ 10^–3, which we keep fixed in all models. We also include a one‐cell meridional circulation pattern having a magnitude of 10–20 m s^–1 near the surface and 1–2 m s^–1 at the bottom of the convection zone. The model essentially represents a distributed turbulent dynamo, as the α ‐effect is nonzero throughout the convection zone, although it concentrates near the bottom of the convection zone obtaining a maximum around 30° of latitude. Turbulent pumping of the mean fields is predominantly down‐ and equatorward. The anisotropies in the turbulent diffusivity are neglected apart from the fact that the diffusivity is significantly reduced in the overshoot region. We find that, when all these effects are included in the model, it is possible to correctly reproduce many features of the solar activity cycle, namely the correct equatorward migration at low latitudes and the polar branch at high latitudes, and the observed negative sign of B _r B _ϕ . Although the activity clearly shifts towards the equator in comparison to previous models due to the combined action of the α ‐effect peaking at midlatitudes, meridional circulation and latitudinal pumping, most of the activity still occurs at too high latitudes (between 5° … 60°). Other problems include the relatively narrow parameter space within which the preferred solution is dipolar (A0), and the somewhat too short cycle lengths of the solar‐type solutions. The role of the surface shear layer is found to be important only in the case where the α ‐effect has an appreciable magnitude near the surface. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Testing turbulent convection theory in solar models – I. Structure of the solar convection zone

Turbulent convection models (TCMs) based on hydrodynamic moment equations are compared with the classical mixing-length theory (MLT) in solar models. The aim is to test the effects of some physical processes on the structure of the solar convection zone, such as the dissipation, diffusion and anisotropy of turbulence that have been ignored in the MLT. Free parameters introduced by the TCMs are also tested in order to find appropriate values for astrophysical applications. It is found that the TCMs usually give larger convective heat fluxes than the MLT does, and the heat transport efficiency is sensitively related to the dissipation parameters used in the TCMs. As a result of calibrating to the present solar values, our solar models usually have rather smaller values of the mixing length to local pressure scaleheight ratio than the standard solar model. The turbulent diffusion is found to have important effects on the structure of the solar convection zone. It leads to significantly lowered and expanded profiles for the Reynolds correlations, and a larger temperature gradient in the central part of the superadiabatic convection region but a smaller one near the boundaries of the convection zone. It is interesting to note that, due to a careful treatment of turbulence developing towards isotropic state, our non-local TCM results in radially dominated motion in the central part and horizontally dominated motion near the boundaries of the convection zone, just as what has been observed in many 3D numerical simulations. Our solar models with the TCMs give small but meaningful differences in the temperature and sound speed profiles compared with the standard solar model using the MLT.     

Helioseismology and the solar interior dynamics

The inversion of helioseismic modes leads to the sound velocity inside the Sun with a precision of about 0.1 per cent. Comparisoons of solar models with the “seismic sun” represent powerful tools to test the physics: depth of the convection zone, equation of state, opacities, element diffusion processes and mixing inside the radiative zone. We now have evidence that microscopic diffusion (element segregation) does occur below the convection zone, leading to a mild helium depletion in the solar outer layers. Meanwhile this process must be slowed down by some macroscopic effect, presumably rotation-induced mixing. The same mixing is also responsible for the observed lithium depletion. On the other hand, the observations of beryllium and helium 3 impose specific constraints on the depth of this mildly mixed zone. Helioseismology also gives information on the internal solar rotation: while differential rotation exists in the convection zone, solid rotation prevails in the radiative zone, and the transition layer (the so-called “tachocline”) is very small. These effects are discussed, together with the astrophysical constraints on the solar neutrino fluxes.     

The slender solar tachocline: a magnetic model

A model for the sharp transition from differential rotation in the solar convection zone to rigid rotation in the radiative interior is presented. Differential rotation in the radiative zone is shown to be quenched efficiently by an internal magnetic field. The poloidal field amplitude, B₀, is the input parameter for our model which determines the transition layer thickness and the toroidal field strength. It is illustrated analytically and confirmed numerically that a rather small field, B₀ = 10⁻⁴ Gauss, suffices to satisfy the helioseismological restrictions on the depth of the differential rotation penetration below the convection zone. The transition layer thickness decreases further with increasing B₀. The toroidal field amplitude B ≃ 200 Gauss is almost independent of B₀.     

Fluctuations of the Solar Dynamo Observed in the Solar Wind and Interplanetary Magnetic Field at 1 AU and in the Outer Heliosphere

Recent helioseismic observations have found strong fluctuations at a period of about 1.3 years in the rotation speed around the tachocline in the deep solar convection layer. Similar mid-term quasi-periodicities (MTQP; periods between 1–2 years) are known to occur in various solar atmospheric and heliospheric parameters for centuries. Since the deep convection layer is the expected location of the solar magnetic dynamo, its fluctuations could modulate magnetic flux generation and cause related MTQP fluctuations at the solar surface and beyond. Accordingly, it is likely that the heliospheric MTQP periodicities reflect similar changes in solar dynamo activity. Here we study the occurrence of the MTQP periodicities in the near and distant heliosphere in the solar wind speed and interplanetary magnetic field observed by several satellites at 1 AU and by four interplanetary probes (Pioneer 10 and 11 and Voyager 1 and 2) in the outer heliosphere. The overall structure of MTQP fluctuations in the different locations of the heliosphere is very consistent, verifying the solar (not heliospheric) origin of these periodicities. We find that the mid-term periodicities were particularly strong during solar cycle 22 and were observed at two different periods of 1.3 and 1.7 years simultaneously. These periodicities were latitudinally organized so that the 1.3-year periodicity was found in solar wind speed at low latitudes and the 1.7-year periodicity in IMF intensity at mid-latitudes. While all heliospheric results on the 1.3-year periodicity are in a good agreement with helioseismic observations, the 1.7-year periodicity has so far not been detected in helioseismic observations. This may be due to temporal changes or due to the helioseismic method where hemispherically antisymmetric fluctuations would so far have remained hidden. In fact, there is evidence that MTQP fluctuations may occur antisymmetrically in the northern and southern solar hemisphere. Moreover, we note that the MTQP pattern was quite different during solar cycles 21 and 22, implying fundamental differences in solar dynamo action between the two halves of the magnetic cycle.     

Realistic Solar Convection Simulations

We report on realistic simulations of solar surface convection that are essentially parameter-free, but include detailed physics in the equation of state and radiative energy exchange. The simulation results are compared quantitatively with observations. Excellent agreement is obtained for the distribution of the emergent continuum intensity, the profiles of weak photospheric lines, the p-mode frequencies, the asymmetrical shape of the mode velocity and intensity spectra, the p-mode excitation rate, and the depth of the convection zone. We describe how solar convection is non-local. It is driven from a thin surface thermal boundary layer where radiative cooling produces low entropy gas which forms the cores of the downdrafts in which most of the buoyancy work occurs. Turbulence and vorticity are mostly confined to the intergranular lanes and underlying downdrafts. Finally, we present some preliminary results on magneto-convection.     

Differential rotation and meridional flow in the solar convection zone and beneath

The influence of the basic rotation on anisotropic and inhomogeneous turbulence is discussed in the context of differential rotation theory. An improved representation for the original turbulence leads to a Λ‐effect which complies with the results of 3D numerical simulations. The resulting rotation law and meridional flow agree well with both the surface observations (∂Ω/∂r < 0 and meridional flow towards the poles) and with the findings of helioseismology. The computed equatorward flow at the bottom of convection zone has an amplitude of about 10 m/s and may be significant for the solar dynamo. The depth of the meridional flow penetration into the radiative zone is proportional to ν^0.5_core, where ν_core is the viscosity beneath the convection zone. The penetration is very small if the tachocline is laminar. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effects of errors in the solar radius on helioseismic inferences

Frequencies of intermediate-degree f modes of the Sun seem to indicate that the solar radius is smaller than what is normally used in constructing solar models. We investigate the possible consequences of an error in radius on results for solar structure obtained using helioseismic inversions. It is shown that solar sound speed will be overestimated if oscillation frequencies are inverted using reference models with a larger radius. Using solar models with a radius of 695.78 Mm and new data sets, the base of the solar convection zone is estimated to be at a radial distance of 0.7135 ± 0.0005 of the solar radius. The helium abundance in the convection zone as determined using models with an OPAL equation of state is 0.248 ± 0.001, where the errors reflect the estimated systematic errors in the calculation, the statistical errors being much smaller. Assuming that the OPAL opacities used in the construction of the solar models are correct, the surface Z / X is estimated to be 0.0245 ± 0.0006.     

The Contribution of Sound Waves and Instabilities to the Penetration of the Solar Differential Rotation Below the Convection Zone

The response of a layer to a horizontal shear flow at its top the surface was studied numerically as an initial value problem. The geometry was Cartesian and the conservation equations were solved with the help of the Zeus-3D code. In the initial state, the pressure, p, and density, ρ, of the layer were assumed to be related by a polytropic equation of index 1.14, which best approximates the solar values in the region of interest. The values of p and ρ at the lower boundary of the layer, namely r=R _l=0.4 R _⊙, were taken to be the solar values. The upper boundary was chosen to be the base of the solar convection zone, r=R _c=0.7 R _⊙. The shear flow at the surface, v _φ(R _c), was proportional to the solar differential rotation, and acoustical oscillations were present in the layer. It is shown that if the initial state is stable, a dynamical coupling between sound waves and the shear flow transmits the surface flow to the inner regions of the layer, even in the absence of dissipation. The shear flow in the sublayer below the one at the surface is proportional to v _φ(R _c), to the time, and to the strength of the oscillations. The constant of proportionality is calculated from the numerical integrations, performed for times of the order of 100 hr. Extrapolation of these results to longer times shows that the surface shear flow is transmitted to the inner regions in a time of the order of of 30 000 years. If the initial state is unstable to the vertical shear, the region of maximum instability depends also on the horizontal shear, and is located away from the equator (where the vertical shear is maximum). As a consequence, the longitudinal flow below the surface shows two equidistant maxima across the equator, located at intermediate latitudes.     

Convection regions and coronas of sun and stars

Photospheric models were calculated for 90 stars with effective temperatures between 2500 K and 41600 K for five logg-values ranging from 1 to 5. Molecule formation was taken into account. In order to have an idea about possible instabilities in the different stellar layers some quantities, characteristic for convection and turbulence were calculated, such as the Rayleigh-, Reynolds-, Prandtl- and Péclet-numbers. It turned out that all the investigated stars contain unstable layers, including the hottest. Nevertheless, only stars with effective temperatures of 8300 K or less contain layers where the convective energy transport is important. For all stars the convective velocities were calculated and also the generated mechanical fluxes in the convection zones were tabulated.Under the hypothesis that this mechanical energy flux is responsible for the heating of the corona, coronal models were constructed for the Sun and for some stars with effective temperatures between 5000 K and 8320 K for logg-values of 4 or 5.For Main Sequence stars the largest fluxes are generated in F-stars; stars withT _eff=7130 K and logg=4 possess also the hottest and most dense coronas with a computed temperature of 3.7·10⁶ K and logN _e =10.5.The solar corona computed in this way, on the basis of a photospheric mechanical flux of 0.14·10⁸ erg cm^–2 sec^–1, has a temperature of 1.3·10⁶ K and logN _e =9.8. This density is apparently too high, but even when including in the computations all theoretical refinements proposed in the last few years by various authors it does not appear possible to obtain a solar coronal model with a smaller density.However, when taking into account the inhomogeneous structure of the chromosphere and by associating the calculated mechanical fluxes to the coarse mottles, and lower fluxes to the undisturbed regions we find a mean coronal temperature of 1.1·10⁶ K and a mean logN _e -values of 9. The computed velocity of the solar wind at a distance of 10⁴ km above the photosphere has a value between 7 and 11 km sec^–1. These latter values are in fair agreement with the observations.     

On the influence of nonlinearities on the eigenfrequencies of five-minute oscillations of the Sun

Fitting the results of linear normal-mode analysis of the solar five-minute oscillations to the observed k - ω diagram selects a class of models of the Sun's envelope. It is a property of all the models in this class that their convection zones are too deep to permit substantial transmission of internal g modes of degree 20 or more. This is in apparent conflict with Hill and Caudell's (1979) claim to have detected such modes in the photosphere. A proposal to resolve the conflict was made by Rosenwald and Hill (1980). They pointed out that despite the impressive agreement between linearized theory and observation, nonlinear phenomena in the solar atmosphere might influence the eigenfrequencies considerably. In particular, they suggested that a correct nonlinear analysis could predict a shallow convection zone. This paper is an enquiry into whether their hypothesis is plausible. We construct k - ω diagrams assuming that the modes suffer local nonlinear distortions in the atmosphere that are insensitive to the amplitude of oscillation over the range of amplitudes that are observed. The effect of the nonlinearities on the eigenfrequencies is parameterized in a simple way. Taking a class of simple analytical models of the Sun's envelope, we compute the linear eigenfrequencies of one model and show that no other model can be found whose nonlinear eigenfrequencies agree with them. We show also that the nonlinear eigenfrequencies of a particular solar model with a shallow convective zone, computed with more realistic physics, cannot be made to agree with observation. We conclude, therefore, that the hypothesis of Rosenwald and Hill is unlikely to be correct.     

Quenching of Meridional Circulation in Flux Transport Dynamo Models

Guided by the recent observational result that the meridional circulation of the Sun becomes weaker at the time of the sunspot maximum, we have included a parametric quenching of the meridional circulation in solar dynamo models such that the meridional circulation becomes weaker when the magnetic field at the base of the convection zone is stronger. We find that a flux transport solar dynamo tends to become unstable on including this quenching of meridional circulation if the diffusivity in the convection zone is less than about 2×10¹¹ cm² s⁻¹. The quenching of α, however, has a stabilizing effect and it is possible to stabilize a dynamo with low diffusivity with sufficiently strong α-quenching. For dynamo models with high diffusivity, the quenching of meridional circulation does not produce a large effect and the dynamo remains stable. We present a solar-like solution from a dynamo model with diffusivity 2.8×10¹² cm² s⁻¹ in which the quenching of meridional circulation makes the meridional circulation vary periodically with solar cycle as observed and does not have any other significant effect on the dynamo.     

The convection zone in stars with radiation pressure in the presence of magnetic field

The effect of a uniform magnetic field on the envelope convection zone of an 8.8M star has been studied. The adiabatic exponents _i (i=1, 2, 3), adiabatic temperature gradient and specific heat of stellar matter has been computed. It is shown that the magnetic field tends to increase the values of adiabatic temperature gradient and specific heats of stellar matter in the envelope convection zone.     

Outer convection zones and internal temperatures of cool white dwarfs

The outer convection zone of the low-temperature white dwarf Van Maanen 2 has been studied for two different atmospheric models given byWeidemann (1960). A slight modification of the standard mixing length theory and the abundances derived by Weidemann have been used.The thickness of the convection zone is about 8 km for the atmospheric model withT _eff=5780 K,g=10⁸ cm sec^–2 and about 23 km forT _eff=5040,g=3.16×10⁷K. In both cases the temperature at the lower boundary of the convection zone is about 9.8×10⁵K. It is shown that this temperature corresponds approximately to the transition temperatureT _tr to the (almost) isothermal core of the white dwarf. This value is considerably lower than the values ofT _tr discussed in the literature until now.The outer convection zone consists of an upper completely non-degenerate part and a lower part with moderate degeneracy. In this lower part the degree of degeneracy is practically independent of depth.     

A New Way that Planets Can Affect the Sun

We derive a perturbation inside a rotating star that occurs when the star is accelerated by orbiting bodies. If a fluid element has rotational and orbital components of angular momentum with respect to the inertially fixed point of a planetary system that are of opposite sign, then the element may have potential energy that could be released by a suitable flow. We demonstrate the energy with a very simple model in which two fluid elements of equal mass exchange positions, calling to mind a turbulent field or natural convection. The exchange releases potential energy that, with a minor exception, is available only in the hemisphere facing the barycenter of the planetary system. We calculate its strength and spatial distribution for the strongest case (“vertical”) and for weaker horizontal cases whose motions are all perpendicular to gravity. The vertical cases can raise the kinetic energy of a few well positioned convecting elements in the Sun’s envelope by a factor ≤7. This is the first physical mechanism by which planets can have a nontrivial effect on internal solar motions. Occasional small mass exchanges near the solar center and in a recently proposed mixed shell centered at 0.16R _s would carry fresh fuel to deeper levels. This would cause stars like the Sun with appropriate planetary systems to burn somewhat more brightly and have shorter lifetimes than identical stars without planets. The helioseismic sound speed and the long record of sunspot activity offer several bits of evidence that the effect may have been active in the Sun’s core, its envelope, and in some vertically stable layers. Additional proof will require direct evidence from helioseismology or from transient waves on the solar surface.     

Modelling the differential rotation of F stars

We model stellar differential rotation based on the mean-field theory of fluid dynamics. DR is mainly driven by Reynolds stress, which is anisotropic and has a non-diffusive component because the Coriolis force affects the convection pattern. Likewise, the convective heat transport is not strictly radial but slightly tilted towards the rotation axis, causing the polar caps to be slightly warmer than the equator. This drives a flow opposite to that caused by differential rotation and so allows the system to avoid the Taylor-Proudman state. Our model reproduces the rotation pattern in the solar convection zone and allows predictions for other stars with outer convection zones. The surface shear turns out to depend mainly on the spectral type and only weakly on the rotation rate. We present results for stars of spectral type F which show signs of very strong differential rotation in some cases. Stars just below the mass limit for outer convection zones have shallow convection zones with short convective turnover times. We find solar-type rotation and meridional flow patterns at much shorter rotation periods and horizontal shear much larger than on the solar surface, in agreement with recent observations. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)