This paper presents the results of a set of radiative hydrodynamic (RHD) simulations of convection in the near-surface regions of a rapidly rotating star. The simulations use microphysics consistent with stellar models, and include the effects of realistic convection and radiative transfer. We find that the overall effect of rotation is to reduce the strength of turbulence. The combination of rotation and radiative cooling creates a zonal velocity profile in which the motion of fluid parcels near the surface is independent of rotation. Their motion is controlled by the strong up and down flows generated by radiative cooling. The fluid parcels in the deeper layers, on the other hand, are controlled by rotation.

We extend the analysis of Penev et al. (2007) to calculate effective viscosities for the surface convective zones of three main sequence stars of 0.775Msun, 0.85Msun and the present day Sun. In addition we also pay careful attention to all normalization factors and assumptions in order to derive actual numerical prescriptions for the effective viscosity as a function of the period and direction of the external shear. Our results are applicable for periods that are too long to correspond to eddies that fall within the inertial subrange of Kolmogorov scaling, but no larger than the convective turnover time, when the assumptions of the calculation break down. We find linear scaling of effective viscosity with period and magnitudes at least three times larger than the Zahn (1966, 1989) prescription.

The development of 2D and 3D simulations of solar convection has lead to a picture of convection quite unlike the usually assumed Kolmogorov spectrum turbulent flow. We investigate the impact of this changed structure on the dissipation properties of the convection zone, parametrized by an effective viscosity coefficient. We use an expansion treatment developed by Goodman & Oh 1997, applied to a numerical model of solar convection (Robinson et al. 2003) to calculate an effective viscosity as a function of frequency and compare this to currently existing prescriptions based on the assumption of Kolmogorov turbulence (Zahn 1966, Goldreich & Keeley 1977). The results match quite closely a linear scaling with period, even though this same formalism applied to a Kolmogorov spectrum of eddies gives a scaling with power-law index of 5/3.

The space telescope MOST is now providing us with extremely accurate low frequency p-mode oscillation data for the star Eta Boo. We demonstrate in this paper that these data, when combined with ground based measurements of the high frequency p-mode spectrum, can be reproduced with stellar models that include the effects of turbulence in their outer layers. Without turbulence, the l=0 modes of our models deviate from either the ground based or the space data by about 1.5-4.0 micro Hz. This discrepancy can be completely removed by including turbulence in the models and we can exactly match 12 out of 13 MOST frequencies that we identified as l=0 modes in addition to 13 out of 21 ground based frequencies within their observational 2 sigma tolerances. The better agreement between model frequencies and observed ones depends for the most part on the turbulent kinetic energy which was taken from a 3D convection simulation for the Sun.

Compared to the Sun, the atmospheric structure and convective flow in Procyon A exhibit the following characteristics: (1) the highly superadiabatic transition layer (SAL) is located at much shallower optical depth; it is in a dynamically active region, and its outer region is located part of the time in the optically thin atmosphere; (2) the outer region of the SAL moves from an optically thin region to thick region and back again over a time of 20-30 minutes. This motion, which is driven by the granulation, takes place in a time approximately half the turnover time of the largest granules; The main reason for the radically different radiative-convective behaviour in Procyon A compared to the Sun is the role played by turbulent eddies in determining the overall flow/thermal structure. The turbulent pressure and turbulent kinetic energy can exceed 50 % of the local gas pressure (compared to about 10-20 % in the Sun). The Procyon A simulation thus reveals two distinct timescales - the autocorrelation time of the vertical velocity and the characteristic timescale of the SAL which is tied to granulation. Just below the surface the autocorrelation decay time is about 5 minutes in Procyon A, and the SAL motion timescale is 20-30 mins. When the SAL penetrates the optically thin region there are efficient radiative losses and the peak of the SAL is low. We speculate that these losses damp out the relative amplitudes in luminosity (temperature fluctuations) compared to velocity (Doppler). Although this will not affect the frequencies of the peaks in the power spectrum, it will probably lower the average amplitude of the peaks relative to the noise background.

This paper describes a series of 3D simulations of shallow inefficient convection in the outer layers of the Sun. The computational domain is a closed box containing the convection-radiation transition layer, located at the top of the solar convection zone. The most salient features of the simulations are that: i)The position of the lower boundary can have a major effect on the characteristics of solar surface convection (thermal structure, kinetic energy and turbulent pressure). ii)The width of the box has only a minor effect on the thermal structure, but a more significant effect on the dynamics (rms velocities). iii)Between the surface and a depth of 1 Mm, even though the density and pressure increase by an order of magnitude, the vertical correlation length of vertical velocity is always close to 600 km. iv) In this region the vertical velocity cannot be scaled by the pressure or the density scale height. This casts doubt on the applicability of the mixing length theory, not only in the superadiabatic layer, but also in the adjacent underlying layers. v) The final statistically steady state is not strictly dependent on the initial atmospheric stratification.

The general consensus is that in order to reproduce the observed solar p-mode oscillation frequencies, turbulence should be included in solar models. However, until now there has not been any well-tested efficient method to incorporate turbulence into solar modeling. We present here two methods to include turbulence in solar modeling within the framework of the mixing length theory, using the turbulent velocity obtained from numerical simulations of the highly superadiabatic layer of the sun at three stages of its evolution. The first approach is to include the turbulent pressure alone, and the second is to include both the turbulent pressure and the turbulent kinetic energy. The latter is achieved by introducing two variables: the turbulent kinetic energy per unit mass, and the effective ratio of specific heats due to the turbulent perturbation. These are treated as additions to the standard thermodynamic coordinates (e.g. pressure and temperature). We investigate the effects of both treatments of turbulence on the structure variables, the adiabatic sound speed, the structure of the highly superadiabatic layer, and the p-mode frequencies. We find that the second method reproduces the SAL structure obtained in 3D simulations, and produces a p-mode frequency correction an order of magnitude better than the first method.

We present results of two simulations of the convection zone, obtained by solving the full hydrodynamic equations in a section of a spherical shell. The first simulation has cylindrical rotation contours (parallel to the rotation axis) and a strong meridional circulation, which traverses the entire depth. The second simulation has isorotation contours about mid-way between cylinders and cones, and a weak meridional circulation, concentrated in the uppermost part of the shell. We show that the solar differential rotation is directly related to a latitudinal entropy gradient, which pervades into the deep layers of the convection zone. We also offer an explanation of the angular velocity shear found at low latitudes near the top. A non-zero correlation between radial and zonal velocity fluctuations produces a significant Reynolds stress in that region. This constitutes a net transport of angular momentum inwards, which causes a slight modification of the overall structure of the differential rotation near the top. In essence, the {\it thermodynamics controls the dynamics through the Taylor-Proudman momentum balance}. The Reynolds stresses only become significant in the surface layers, where they generate a weak meridional circulation and an angular velocity `bump'.

We have constructed models for the sun at three stages of its evolution: a zero-age main sequence model, the present sun, and a subgiant model. For each model, the turbulent pressure and turbulent kinetic energy were calculated from 3-d radiative hydrodynamical simulations (described in the poster by Robinson et al.), and inserted into the 1-d stellar models. We note that in these simulations, the turbulent pressure is not a free parameter, but can be computed from the resulting velocity field. We show the calculated p-mode frequencies for the model of the present sun, with and without turbulent pressure, and compare them to the observed solar frequencies. When the turbulent pressure is included in the models, the calculated frequencies are brought closer to the observed frequencies in the sun by up to two $\mu Hz$, strictly from structural effects. The effect of including turbulent pressure on p-mode frequencies is also shown for the zero-age main sequence model. Our models also suggest that the importance of turbulent pressure increases as the star evolves into the subgiant region. We discuss the importance of also including realistic turbulence as well as radiation in the non-adiabatic calculation of oscillation frequencies.

We present results of three dimensional simulations of the uppermost part of the sun, at 3 stages of its evolution. Each model includes physically realistic radiative-hydrodynamics (the Eddington approximation is used in the optically thin region), varying opacities and a realistic equation of state (full treatment of the ionization of H and He). In each evolution model, we investigate a domain, which starts at the top of the photosphere and ends just inside the convection zone (about 2400 km in the sun model). This includes all of the super-adiabatic layer (SAL). Due to the different positions of the three models in the $log (g) $ vs $log T_{eff}$ plane, the more evolved models have lower density atmospheres. The reduction in density causes the amount of overshoot into the radiation layer, to be greater in the more evolved models.