The emergence of minimally structured (uniform and horizontal) magnetic field from a depth of 20 Mm has been simulated. The field emerges first in a mixed polarity pepper and salt pattern, but then collects into separate, unipolar concentrations and produces pores. The field strength was then artificially increased to produce spot-like structures. The field strength at continuum optical depth unity peaks at 1 kG, with a maximum of 4 kG. Where the vertical field is strong, the spots persist (at present an hour of solar time has been simulated). Where the field is weak, the spot gets filled in and disappears. Stokes profiles have been calculated and processed with the Hinode annular mtf, the slit diffraction and frequency smoothing. These data are available at steinr.pa.msu.edu/~bob/stokes.

We simulate the rise through the upper convection zone and emergence through the solar surface of initially uniform, untwisted, horizontal magnetic flux with the same entropy as the non-magnetic plasma that is advected into a domain 48 Mm wide from from 20 Mm deep. The magnetic field is advected upward by the diverging upflows and pulled down in the downdrafts, which produces a hierarchy of loop like structures of increasingly smaller scale as the surface is approached. There are significant differences between the behavior of fields of 10 kG and 20 or 40 kG strength at 20 Mm depth. The 10 kG fields have little effect on the convective flows and show little magnetic buoyancy effects, reaching the surface in the typical fluid rise time from 20 Mm depth of 32 hours. 20 and 40 kG fields significantly modify the convective flows, leading to long thin cells of ascending fluid aligned with the magnetic field and their magnetic buoyancy makes them rise to the surface faster than the fluid rise time. The 20 kG field produces a large scale magnetic loop that as it emerges through the surface leads to the formation of a bipolar pore-like structure.

Results of realistic simulations of solar surface convection on the scale of supergranules (96 Mm wide by 20 Mm deep) are presented. The simulations cover only 10% of the geometric depth of the solar convection zone, but half its pressure scale heights. They include the hydrogen, first and most of the second helium ionization zones. The horizontal velocity spectrum is a power law and the horizontal size of the dominant convective cells increases with increasing depth. Convection is driven by buoyancy work which is largest close to the surface, but significant over the entire domain. Close to the surface buoyancy driving is balanced by the divergence of the kinetic energy flux, but deeper down it is balanced by dissipation. The damping length of the turbulent kinetic energy is 4 pressure scale heights. The mass mixing length is 1.8 scale heights. Two thirds of the area is upflowing fluid except very close to the surface. The internal (ionization) energy flux is the largest contributor to the convective flux for temperatures less than 40,000 K and the thermal energy flux is the largest contributor at higher temperatures. This data set is useful for validating local helioseismic inversion methods. Sixteen hours of data are available as four hour averages, with two hour cadence, at steinr.msu.edu/~bob/96averages, as idl save files. The variables stored are the density, temperature, sound speed, and three velocity components. In addition, the three velocity components at 200 km above mean continuum optical depth unity are available at 30 sec. cadence.

We apply time-distance helioseismology, local correlation tracking and Fourier spatial-temporal filtering methods to realistic supergranule scale simulations of solar convection and compare the results with high-resolution observations from the SOHO Michelson Doppler Imager (MDI). Our objective is to investigate the surface and sub-surface convective structures and test helioseismic measurements. The size and grid of the computational domain are sufficient to resolve various convective scales from granulation to supergranulation. The spatial velocity spectrum is approximately a power law for scales larger than granules, with a continuous decrease in velocity amplitude with increasing size. Aside from granulation no special scales exist, although a small enhancement in power at supergranulation scales can be seen. We calculate the time-distance diagram for f- and p-modes and show that it is consistent with the SOHO/MDI observations. From the simulation data we calculate travel time maps for surface gravity waves (f-mode). We also apply correlation tracking to the simulated vertical velocity in the photosphere to calculate the corresponding horizontal flows. We compare both of these to the actual large-scale (filtered) simulation velocities. All three methods reveal similar large scale convective patterns and provide an initial test of time-distance methods.

Recent progress in realistic simulations of solar convection have given us an unprecedented opportunity to evaluate the robustness of solar interior structures and dynamics obtained by methods of local helioseismology. We present results of testing the time-distance method using realistic simulations. By computing acoustic wave propagation time and distance relations for different depths of the simulated data, we confirm that acoustic waves propagate into the interior and then turn back to the photosphere. This demonstrates that in the numerical simulations properties of acoustic waves (p-modes) are similar to the solar conditions, and that these properties can be analyzed by the time-distance technique. For the surface gravity waves (f-mode), we calculate perturbations of their travel times, caused by localized downdrafts, and demonstrate that the spatial pattern of these perturbations (representing so-called sensitivity kernels) is similar to the patterns obtained from the real Sun, displaying characteristic hyperbolic structures. We then test the time-distance measurements and inversions by calculating acoustic travel times from a sequence of vertical velocities at the photosphere of the simulated data, and inferring a mean 3D flow fields by performing inversion based on the ray approximation. The inverted horizontal flow fields agree very well with the simulated data in subsurface areas up to 3 Mm deep, but differ in deeper areas. Due to the cross-talk effects between the horizontal divergence and downward flows, the inverted vertical velocities are significantly different from the mean convection velocities of the simulation dataset.

We extend semi-analytical computations of excitation rates for solar oscillation modes to those of other solar-like oscillating stars to compare them with recent observations. Numerical 3D simulations of surface convective zones of several solar-type oscillating stars are used to characterize the turbulent spectra as well as to constrain the convective velocities and turbulent entropy fluctuations in the uppermost part of the convective zone of such stars. These constraints, coupled with a theoretical model for stochastic excitation, provide the rate 'P' at which energy is injected into the p-modes by turbulent convection. These energy rates are compared with those derived directly from the 3D simulations. The excitation rates obtained from the 3D simulations are systematically lower than those computed from the semi-analytical excitation model. We find that Pmax, the excitation rate maximum, scales as (L/M)^s where s is the slope of the power law and L and M are the mass and luminosity of the 1D stellar model built consistently with the associated 3D simulation. The slope is found to depend significantly on the adopted form of the eddy time-correlation ; using a Lorentzian form results in s=2.6, whereas a Gaussian one gives s=3.1. Finally, values of Vmax, the maximum in the mode velocity, are estimated from the computed power laws for Pmax and we find that Vmax increases as (L/M)^sv. Comparisons with the currently available ground-based observations show that the computations assuming a Lorentzian eddy time-correlation yield a slope, sv, closer to the observed one than the slope obtained when assuming a Gaussian. We show that the spatial resolution of the 3D simulations must be high enough to obtain accurate computed energy rates.

We compute the rates P at which energy is injected into the p modes for a set of 3D simulations of outer layers of stars. We found that Pmax - the maximum in P - scales as (L/M)^s where s is the slope of the power law, L and M are the luminosity and the mass of the 1D stellar models associated with the simulations. The slope is found to depend significantly on the adopted representation for the turbulent eddy-time correlation function, chi_k. According to the expected performances of COROT, it will likely be possible to measure Pmax as a function of L/M and to constrain the properties of stellar turbulence as the turbulent eddy time-correlation.

Hydrodynamical, 3D simulations of the outer layers of the Sun and Alpha Cen A are used to obtain constraints on the properties of turbulent convection in such stars. These constraints enable us to compute - on the base of a theoretical model of stochastic excitation - the rate P at which p modes are excited by turbulent convection in those two stars. Results are then compared with solar seismic observations and recent observations of Alpha Cen A. For the Sun, a good agreement between observations and computed P is obtained. For Alpha Cen A a large discrepancy is obtained which origin cannot be yet identified: it can either be caused by the present data quality which is not sufficient for our purpose or by the way the intrinsic amplitudes and the life-times of the modes are determined or finally attributed to our present modelling. Nevertheless, data with higher quality or/and more adapted data reductions will likely provide constraints on the p-mode excitation mechanism in Alpha Cen A.

The solar acoustic p-mode line profiles are asymmetric. Velocity spectra have more power on the low-frequency sides, whereas intensity profiles show the opposite sense of asymmetry. Numerical simulations of the upper convection zone have resonant p-modes with the same asymmetries and asymmetry reversal as the observed modes. The temperature and velocity power spectra at optical depth $\tau_{\rm cont} = 1$ have the opposite asymmetry as is observed for the intensity and velocity spectra. At a fixed geometrical depth, corresponding to $<\tau_{\rm cont}>=1$, however, the temperature and velocity spectra have the same asymmetry. This indicates that the asymmetry reversal is produced by radiative transfer effects and not by correlated noise.

We use the three-dimensional hydrodynamic code of Stein and Nordlund to realistically simulate the upper layers of the solar convection zone in order to study physical characteristics of solar oscillations. Our first result is that the properties of oscillation modes in the simulation closely match the observed properties. Recent observations from SOHO/MDI and GONG have confirmed the asymmetry of solar oscillation line profiles, initially discovered by Duvall et al. In this paper we compare the line profiles in the power spectra of the Doppler velocity and continuum intensity oscillations from the SOHO/MDI observations with the simulation. We also compare the phase differences between the velocity and intensity data. We have found that the simulated line profiles are asymmetric and have the same asymmetry reversal between velocity and intensity as observed. The phase difference between the velocity and intensity signals is negative at low frequencies and jumps in the vicinity of modes as is also observed. Thus, our numerical model reproduces the basic observed properties of solar oscillations, and allows us to study the physical properties which are not observed.