Anisotropy of turbulent transport in a convective layer

An estimate for the anisotropy of the turbulent viscositys is given in a convective layer heated from below and rotating around a vertical axis. In the case of two-dimensional convection, there is a stationary regime withs⊇2 regardless of the rotation. In the case of three-dimensional convection in a slowly rotating layer (with the Taylor number equal to 1600), nonstationary turbulent regimes take place withs⊇1.6 forR=2.5×10⁴ (R is the Rayleigh number) ands⊇1.2 forR=10⁴. The parameters plays an, important role in the theory of differential rotation of the convective solar or stellar envelopes. So far, it has been evaluated empirically or semi-empirically. Some prospects in the development of the theory of differential rotation are discussed here in terms of the moment theory of hydrodynamic fields. The relation between this strict approach and an anisotropic viscosity approximation is considered.     

Tachocline Confinement by an Oscillatory Magnetic Field

Helioseismic measurements indicate that the solar tachocline is very thin, its full thickness not exceeding 4% of the solar radius. The mechanism that inhibits differential rotation to propagate from the convective zone to deeper into the radiative zone is not known, though several propositions have been made. In this paper we demonstrate by numerical models and analytic estimates that the tachocline can be confined to its observed thickness by a poloidal magnetic field B _p of about one kilogauss, penetrating below the convective zone and oscillating with a period of 22 years, if the tachocline region is turbulent with a diffusivity of η∼10¹⁰ cm² s⁻¹ (for a turbulent magnetic Prandtl number of unity). We also show that a similar confinement may be produced for other pairs of the parameter values (B _p, η). The assumption of the dynamo field penetrating into the tachocline is consistent whenever η≳10⁹ cm² s⁻¹. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1013389631585     

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)     

On the sun's pole-equator flux differences

We study the possibility that large flux differences between the poles and the equator at the bottom of the solar convective zone are compatible with the small differences observed at the surface. The consequences of increasing the depth of the convective zone due to overshooting are explored.A Boussinesq model is used for the convective zone and we assume that the interaction of the global convection with rotation is modelled through a convective flux coefficient whose perturbed part is proportional to the local Taylor number. The numerical integration of the equations of motion and energy shows that coexistence between large pole-equator flux differences at the bottom and small ones at the surface is possible if the solar convective zone extends to a depth of 0.4R . The angular velocity distribution inside the convective zone is in agreement with the -dynamo theories of the solar cycle.     

Effect of Gravity Waves on Transport Process in Stellar Radiative Region

The aim of this work is to present a transport process which is likely to have a great importance for the internal constitution of the stars. In order to set the problem, we first give a short presentation of the physical properties of the Sun and stars, described usually under the names of `Standard Solar Model' or `Standard Stellar Models' (SSM). Next we show that an important question about SSM is that they do not explain the age dependance of lithium deficiency of stars of known age: stars of galactic clusters and the Sun. It has been suggested a long time ago to assume the presence of a macrosocpic diffusion process in the radiative zone, below the surface convective zone of solar like stars. It is then possible for the lithium present in the convective zone to be carried to the thermonuclear burning level below the convective zone. The first assumption was that differential rotation generates turbulence and therefore that a turbulent diffusion process takes place. However, this model predicts a lithium abundance which is strongly rotation dependant, contrary to the observations. Furthermore, the diffusion coefficient being large all over the radiative zone, it prevents the possibility of gravitational separation by diffusion and consequently leads to an impossibility of explaining the difference of helium abundance between the surface and the center of the Sun. The consequence is obviously that we need to take into account another physical process. Stars having a mass M < 1.3 M _⊙ have a convective zone which begins close to the stellar surface and extends down to a depth which is an appreciable fraction of stellar radius. In the convective zone, strong stochastic motions take care, at least partially, of heat transfer. These motions do not vanish at the lower boundary and generate internal waves into the radiative zone. These random internal waves are at the origin of a diffusion process which can be considered as responsible of the diffusive transport of lithium down to the lithium burning level. This is certainly not the only physical process responsible of lithium deficiency in main sequence stars, but its properties open the way to a completely consistent analysis of lithium deficiency. The model of generation of gravity waves is based on a model of heat transport in the convective zone by diving plumes. The horizontal component of the turbulent motion at the boundary of the convective zone is supposed to generate the horizontal motion of internal waves. The result is a large horizontal component of the diffusion coefficient, which produces in a short time an horizontal uniform chemical composition. It is known that gravity waves, in the absence of any dissipative process, cannot generate vertical mixing. Therefore, the vertical component of the diffusion coefficient is entirely dependant of radiative damping. It decreases quickly in the radiative zone, but is large enough to be responsible of lithium burning. Due to the radial dependance of velocity amplitude, the diffusion coeficient increases when approaching the stellar center. However, very close to the center, non-linear dissipative and radiative damping of internal waves become large and the diffusion coefficient vanishes at the very center. The development of this abstract can be found in E. Schatzman (1996, J. Fluid Mech. 322, 355). This revised version was published online in July 2006 with corrections to the Cover Date.     

Seismology of the solar convection zone

An attempt is made to infer the structure of the solar convection zone from observedp-mode frequencies of solar oscillations. The differential asymptotic inversion technique is used to find the sound speed in the solar envelope. It is found that envelope models which use the Canuto-Mazzitelli (CM) formulation for calculating the convective flux give significantly better agreement with observations than models constructed using the mixing length formalism. This inference can be drawn from both the scaled frequency differences and the sound speed difference. The sound speed in the CM envelope model is within 0.2% of that in the Sun except in the region withr > 0.99R _⊙. The envelope models are extended below the convection zone, to find some evidence for the gravitational settling of helium beneath the base of the convection zone. It turns out that for models with a steep composition gradient below the convection zone, the convection zone depth has to be increased by about 6 Mm in order to get agreement with helioseismic observations.     

Transport Effects in the Evolution of the Global Solar Magnetic Field

The axisymmetric component of the large-scale solar magnetic fields has a pronounced poleward branch at higher latitudes. In order to clarify the origin of this branch we construct an axisymmetric model of the passive transport of the mean poloidal magnetic field in the convective zone, including meridional circulation, anisotropic diffusivity, turbulent pumping and density pumping. For realistic values of the transport coefficients we find that diffusivity is prevalent, and the latitudinal distribution of the field at the surface simply reflects the conditions at the bottom of the convective zone. Pumping effects concentrate the field to the bottom of the convective zone; a significant part of this pumping occurs in a shallow subsurface layer, normally not resolved in dynamo models. The phase delay of the surface poloidal field relative to the bottom poloidal field is found to be small. These results support the double dynamo wave models, may be compatible with some form of a mixed transport scenario, and exclude the passive transport theory for the origin of the polar branch.     

The magnetic field transfer in the solar convective zone

The large-scale azimuth magnetic field is pumping to the bottom of the solar convective zone due to the diamagnetic action of turbulent conductive fluids. When the field at the bottom is of about 10³ G, an equilibrium is established between diamagnetic pumping and buoyancy.If, in addition to the density gradient, an additional anisotropy exists (for instance, due to rotation), another mechanism of the magnetic field transfer appears, the efficiency of which greatly depends on the magnitude of the anistropy parameter.     

Smooth models of overshooting at the base of the solar convective zone

A set of smoothed temperature gradient profiles around overshooting layers at the solar convective zone bottom is considered. In classical local theories of convection the one point defined according to the Schwarzschild criterion is enough to describe a convective boundary. To get a sophisticated picture of the overshooting we use four points to compute the transition overshooting functions. Analyzing the transition gradient profiles we found that the overshooting convective flux may be either positive or negative. A negative overshooting flux appears in nonlocal convective theories and causes a steep temperature gradient profile. But we propose an evenly smoothed gradient which corresponds to a convective flux positive everywhere. To outline the effect of the temperature gradient on the solar oscillations the squared Brunt–Väisälä frequency N ² is calculated. In local convective theories the N ² profile shows the discontinuity of the first derivative at the convective boundary, while all smoothed profiles eliminate the break.     

Present problems of the solar interior

The standard model of solar evolution is reviewed and a number of problems highlighted. A fundamental question is whether there is any mixing of matter in the central regions, since such mixing could radically alter the model of the present Sun and modify our understanding of the evolution of other stars. Standard models of solar evolution become unstable to ³He driven global oscillations at an age of 3 × 10⁸ years and this may drive some mixing, even if this is not the case the finite amptitude limit of these oscillations is likely to produce modifications in the standard model. Convective overshooting at the bottom of the outer convective zone leads to an increased depth of this zone and small changes in the interior. It is pointed out that the young Sun had a ¹²C driven convective core whose extent and duration depends on the extent of overshooting. Such a core is likely to produce a magnetic field which will affect the internal dynamics. The internal rotation of the Sun remains an enigma and absence of knowledge of any internal magnetic field makes it difficult to study the problem. Rotationally driven instabilities are ineffective in the central chemically inhomogenous regions but may contribute to the inward diffusion of lithium from the convective zone. These and other problems are considered, but few solutions are proposed.     

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.     

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)     

太阳大气锂丰度的衰减

太阳大气锂的丰度７Ｌｉ／Ｈ＝１０－１１（按原子数计）。或［７Ｌｉ］＝ｌｏｇ（７Ｌｉ／Ｈ）＋１２＝１０．它比太阳系原始星云和银河系星际介质钾的丰度要低约两个数量级．因此太阳在它形成之后，其大气锂必定经受了严重衰减．然而年轻的银河疏散星团（如昂星团和英仙ａ星团）中有效温度高于５５００ Ｋ的主序星，其锂丰度都基本是正常的，井末呈现明显的衰减．这充分说明，太阳型恒星锂的衰减主要发生在主序阶段，而非在主序前的演化阶段． 在恒星中，７Ｌｉ是通过核反应７Ｌｉ（ｐ，ａ）４Ｈｅ，而毁坏．上述反应在Ｔ≥ ２．５×１０…     

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)     

Inhomogeneous convection and the equatorial acceleration of the sun

The interaction of rotation and turbulent convection is assumed to give rise to an inhomogeneous, but isotropic, latitude dependent turbulent energy transport, which is described by a convective conduction coefficient _c which varies with latitude. Energy balance in the convective zone is then possible only with a slow meridian circulation in the outer convective zone of the sun. The angular momentum transported by this circulation is balanced in a steady state by turbulent viscous transport down an angular velocity gradient. A detailed model is constructed allowing for the transition from convective transport to radiative transport at the boundaries of the convective zone, by using a perturbation analysis in which the latitude variation of _c is small. The solution for a thin compressible shell gives equatorial acceleration and a hotter equator than pole, assuming that the convection is preferentially stabilised at the equator. For agreement with the sun's equatorial acceleration the model predicts an equatorial temperature excess of 70 K and a surface meridional velocity of 350 cm/sec from pole to equator.     

Three-dimensional model for generation of the mean solar magnetic field

A three-dimensional (non-axisymmetric) model for the solar mean magnetic field generation is studied. The sources of generation are the differential rotation and mean helicity in the convective shell. The system is described by two equations of the first order in time and the fourth order in space coordinates. The solution is sought for in the form of expansion over the spherical function Y_n^m. The modes of different m are separated. A finite-difference scheme similar to the Peaceman-Rachford scheme is constructed so to find coefficients of the expansion depending on the time and radial coordinates. It is shown that a mode with a smaller azimuthal number m is primarily excited. The axisymmetric mode m = o describes the 22 year solar cycle oscillations. The modes of m o have no such periodicity, the oscillate with a period of rotation of the low boundary of the solar convective shell, The solutions which are symmetric relative to the equator plane are excited more easily compared with the antisymmetrical ones. The results obtained are confronted to the observational picture of the non-axisymmetric large-scale solar magnetic fields.     

Quick matching technique to study the relationship between solar radius and luminosity variations

A simple matching technique is developed which allows us to compute the response of the solar envelope to perturbations which occur within the solar convective region, and in timescales of importance to climate. This technique is applied to perturbation of the convective efficiency (-mechanism), and of the non-gas component of the pressure in different regions of the convection zone (-mechanism). The results indicate that while either perturbation affects the solar luminosity, the -mechanism has almost no effect on the solar outer radius, regardless of the affected region, whereas the -mechanism produces radius changes which may be large if the location of perturbation is deep enough.     

Differential rotation caused by anisotropic turbulent viscosity

The law of rotation as well as the corresponding meridional circulations in the hydrogen convection zone (HCZ) are investigated by solving numerically the time independent Navier-Stokes equations. The HCZ is assumed to be a spherical layer of fluid with constant density and viscosity. It is assumed further that the viscosity is caused by unisotropic turbulent motions.The results show differential rotation together with circulations. The detailed behaviour depends on a parameters characterizing the nonisotropic friction and on the kinematic viscosityv. If the friction is larger in radial direction than in lateral directions (0 s < 1) the poles rotate faster than the equator and the circulation rises at the equator and falls at the poles; if friction is smaller in radial direction (s > 1) the equator rotates faster and the sense of the circulation is reversed. The differential rotation observed at the solar surface is obtained for the values = 1.2.For small values ofv the angular velocity is constant on cylindrical surfaces, for large values ofv it is constant on spherical surfaces. The solar law of rotation turns out to be very close to the first case.Based on the author's Thesis in Göttingen.     

Solar abundances of light nuclei and mixing of the Sun

Radial profiles of the light nuclei (A 15) are calculated in the non-mixing Sun, taking into account the changes of solar structure with time. The results are discussed in relation to models of solar mixing and compared with abundance determinations at the solar surface or in the solar wind. B cannot be depleted in the outer convective zone without producing a large increase in the He³/He⁴ ratio. A decrease in He³/He⁴ would be accompanied by changes in C¹³/C¹² and N¹⁵/N¹⁴ of a magnitude which is not observed.It is shown that boron could be depleted in the pre-main sequence period of the Sun, if mixing was on a time-scale of 10⁶ yr. The simultaneous small increase in He³/He⁴ does not contradict observation. However, Be would be depleted more strongly than B.A He³/He⁴ decrease is always accompanied by large changes in N¹⁵/N¹⁴ and C¹³/C¹². Since such changes are not observed, it is concluded that the He³/He⁴ ratio in the outer convective zone is a reliable upper limit for (He³ + D)/He⁴ in the solar nebula. Thus the D/H ratio in the protosolar material was much lower than it is in sea water or in carbonaceous chondrites.     

Properties of r-modes in the Sun

Global oscillations of the Sun (r-modes) with very long periods 1 month are reviewed and studied. Such modes would be trapped in an acoustic cavity formed either by most of the convective envelope or by most of the radiative interior. A turning point frequency giving cavity boundaries is defined and the run of eigenvalues for angular harmonics l 3 are plotted for a conventional solar convection zone. The r-modes show equipartition of oscillatory energy among shells which each contain one antinode in the radial dimension. Toroidal motion is dominant to at least the 14^th radial harmonic mode. Viscosity from convective turbulence is strong and would damp any mode in just a few solar rotations if it were the only significant nonadiabatic effect. Radial fine splitting which lifts the degeneracy in n is very small (20 nHz or less) for all n 14 trapped in the envelope. But, if splitting could be detected, we would have a valuable new constraint on solar convection theories.