It is generally accepted that solar acoustic (p) modes are excited by near-surface turbulent motions, in particular, by downdrafts and interacting vortices in intergranular lanes. Recent analysis of Solar Dynamics Observatory data by (Zhao et al., 2015) revealed fast-moving waves around sunspots, which are consistent with magnetoacoustic waves excited approximately 5 Mm beneath the sunspot. We analyzed 3D radiative MHD simulations of solar magnetoconvection with a self-organized pore-like magnetic structure, and identified more than 600 individual acoustic events both inside and outside of this structure. By performing a case-by-case study, we found that surrounding the magnetic structure, acoustic sources are associated with downdrafts. Their depth correlates with downdraft speed and magnetic fields. The sources often can be transported into deeper layers by downdrafts. The wave front shape, in the case of a strong or inclined downdraft, can be stretched along the downdraft. Inside the magnetic structure, excitation of acoustic waves is driven by converging flows. Frequently, strong converging plasma streams hit the structure boundaries, causing compressions in its interior that excite acoustic waves. Analysis of the depth distribution of acoustic events shows the strongest concentration at 0.2 - 1 Mm beneath the surface for the outside sources and 2.5 - 3 Mm for the excitation event inside the structure.

Continued progress in observational stellar astrophysics requires a deep understanding of the underlying convection dynamics. We present results of realistic 3D radiative hydrodynamic simulations of the outer layers of a moderate mass star (1.47 Msun), including the full convection zone, the overshoot region, and the top layers of the radiative zone. The simulation results show that the surface granulation has a broad range of scales, from 2 to 12 Mm, and that large granules are organized in well-defined clusters, consisting of several granules. Comparison of the mean structure profiles from 3D simulations with the corresponding 1D standard stellar model shows an increase of the stellar radius by ~800 km, as well as significant changes in the thermodynamic structure and turbulent properties of the ionization zones. Convective downdrafts in the intergranular lanes between granulation clusters reach speeds of more than 20 km/s, penetrate through the whole convection zone, hit the radiative zone, and form a 8 Mm thick overshoot layer. Contrary to semi-empirical overshooting models, our results show that the 3D dynamic overshoot region consists of two layers: a nearly adiabatic extension of the convection zone and a deeper layer of enhanced subadiabatic stratification. This layer is formed because of heating caused by the braking of the overshooting convective plumes. This effect has to be taken into account in stellar modeling and the interpretation of asteroseismology data. In particular, we demonstrate that the 3D model can qualitatively explain deviations from the standard solar model found by helioseismology.

Realistic numerical simulations, i.e., those that make minimal use of ad hoc modeling, are essential for understanding the complex turbulent dynamics of the interiors and atmospheres of the Sun and other stars and the basic mechanisms of their magnetic activity and variability. The goal of this paper is to present a detailed description and test results of a compressible radiative MHD code, `StellarBox', specifically developed for simulating the convection zones, surface, and atmospheres of the Sun and moderate-mass stars. The code solves the three-dimensional, fully coupled compressible MHD equations using a fourth-order Pad\'e spatial differentiation scheme and a fourth-order Runge-Kutta scheme for time integration. The radiative transfer equation is solved using the Feautrier method for bi-directional ray tracing and an opacity-binning technique. A specific feature of the code is the implementation of subgrid-scale MHD turbulence models. The data structures are automatically configured, depending on the computational grid and the number of available processors, to achieve good load balancing. We present test results and illustrate the code's capabilities for simulating the granular convection on the Sun and a set of main-sequence stars. The results reveal substantial changes in the near-surface turbulent convection in these stars, which in turn affect properties of the surface magnetic fields. For example, in the solar case initially uniform vertical magnetic fields tend to self-organize into compact (pore-like) magnetic structures, while in more massive stars such structures are not formed and the magnetic field is distributed more-or-less uniformly in the intergranular lanes.

Magnetic fields are usually observed in the quiet Sun as small-scale elements that cover the entire solar surface (the `salt and pepper' patterns in line-of-sight magnetograms). By using 3D radiative MHD numerical simulations we find that these fields result from a local dynamo action in the top layers of the convection zone, where extremely weak 'seed' magnetic fields (e.g., from a $10^{-6}$ G) can locally grow above the mean equipartition field, to a stronger than 2000~G field localized in magnetic structures. Our results reveal that the magnetic flux is predominantly generated in regions of small-scale helical downflows. We find that the local dynamo action takes place mostly in a shallow, about 500~km deep, subsurface layer, from which the generated field is transported into the deeper layers by convective downdrafts. We demonstrate that the observed dominance of vertical magnetic fields at the photosphere and horizontal fields above the photosphere can be explained by small-scale magnetic loops produced by the dynamo. Such small-scale loops play an important role in the structure and dynamics of the solar atmosphere and that their detection in observations is critical for understanding the local dynamo action on the Sun.

In the quiet Sun, magnetic fields are usually observed as small-scale magnetic elements, `salt and pepper', covering the entire solar surface. By using 3D radiative MHD numerical simulations we demonstrate that these fields are a result of local dynamo action in the top layers of the convection zone, where extremely weak `seed' magnetic fields can locally grow above the mean equipartition field (e.g., from a $10^{-6}$ G `seed' field to more than 1000 G magnetic structures). We find that the local dynamo action takes place only in a shallow, about 500 km deep, subsurface layer, from which the generated field is transported into deeper layers by convection downdrafts. We demonstrate that the observed dominance of vertical magnetic fields at the photosphere and the horizontal fields above the photosphere can be explained by multi-scale magnetic loops produced by the dynamo.

The solar surface is covered by high-speed jets transporting mass and energy into the solar corona and feeding the solar wind. The most prominent of these jets have been known as spicules. However, the mechanism initiating these eruptions events is still unknown. Using realistic numerical simulations we find that small-scale eruptions are produced by ubiquitous magnetized vortex tubes generated by the Sun's turbulent convection in subsurface layers. The swirling vortex tubes (resembling tornadoes) penetrate into the solar atmosphere, capture and stretch background magnetic field, and push surrounding material up, generating quasiperiodic shocks. Our simulations reveal a complicated high-speed flow patterns, and thermodynamic and magnetic structure in the erupting vortex tubes. We found that the eruptions are initiated in the subsurface layers and are driven by the high-pressure gradients in the subphotosphere and photosphere, and by the Lorentz force in the higher atmosphere layers.

Turbulent properties of the quiet Sun represent the basic state of surface conditions, and a background for various processes of solar activity. Therefore understanding of properties and dynamics of this `basic' state is important for investigation of more complex phenomena, formation and development of observed phenomena in the photosphere and atmosphere. For characterization of the turbulent properties we compare kinetic energy spectra on granular and sub-granular scales obtained from infrared TiO observations with the New Solar Telescope (Big Bear Solar Observatory) and from 3D radiative MHD numerical simulations ('SolarBox' code). We find that the numerical simulations require a high spatial resolution with 10 - 25 km grid-step in order to reproduce the inertial (Kolmogorov) turbulence range. The observational data require an averaging procedure to remove noise and potential instrumental artifacts. The resulting kinetic energy spectra show a good agreement between the simulations and observations, opening new perspectives for detailed joint analysis of more complex turbulent phenomena on the Sun, and possibly on other stars. In addition, using the simulations and observations we investigate effects of background magnetic field, which is concentrated in self-organized complicated structures in intergranular lanes, and find an increase of the small-scale turbulence energy and its decrease at larger scales due to magnetic field effects.

We use 3D radiative MHD simulations to investigate the formation and dynamics of small-scale (less than 0.5 Mm in diameter) vortex tubes spontaneously generated by turbulent convection in quiet-Sun regions with initially weak mean magnetic fields. The results show that the vortex tubes penetrate into the chromosphere and substantially affect the structure and dynamics of the solar atmosphere. The vortex tubes are mostly concentrated in intergranular lanes and are characterized by strong (near sonic) downflows and swirling motions that capture and twist magnetic field lines, forming magnetic flux tubes that expand with height and which attain magnetic field strengths ranging from 200 G in the chromosphere to more than 1 kG in the photosphere. We investigate in detail the physical properties of these vortex tubes, including thermodynamic properties, flow dynamics, and kinetic and current helicities, and conclude that magnetized vortex tubes provide an important path for energy and momentum transfer from the convection zone into the chromosphere.

Investigation of the turbulent properties of solar convection is extremely important for understanding the multi-scale dynamics observed on the solar surface. In particular, recent high-resolution observations have revealed ubiquitous vortical structures, and numerical simulations have demonstrated links between vortex tube dynamics and magnetic field organization and have shown the importance of vortex tube interactions in the mechanisms of acoustic wave excitation on the Sun. In this paper we investigate the mechanisms of the formation of vortex tubes in highly-turbulent convective flows near the solar surface by using realistic radiative hydrodynamic LES simulations. Analysis of data from the simulations indicates two basic processes of vortex tube formation: 1) development of small-scale convective instability inside convective granules, and 2) a Kelvin-Helmholtz type instability of shearing flows in intergranular lanes. Our analysis shows that vortex stretching during these processes is a primary source of generation of small-scale vorticity on the Sun.

In this study we investigate the effects of turbulent convection on formation of large-scale inhomogeneous magnetic structures by means of Large-Eddy Simulation (LES) for convection in solar-type stars. The main idea of this study is the implementation of a new subgrid-scale model for the effective Lorentz force in a three-dimensional nonlinear radiative magnetohydrodynamics (MHD) code developed for simulating the upper solar convection zone and lower atmosphere. To this end we derived the energy budget equations, which include the effects of the subgrid-scale turbulence on the Lorentz-force, and implemented the new subgrid-scale turbulence model (TELF-Model) in a three-dimensional nonlinear MHD LES code. Using imposed initial vertical and horizontal uniform magnetic fields in LES with the TELF-Model, we have shown that the magnetic flux tubes formation is started when the initial mean magnetic field is larger than a threshold value (about 100 G). This is in agreement with the theoretical studies by Rogachevskii and Kleeorin (2007). We have determined the vertical profiles of the velocity and magnetic fluctuations, total MHD energy and anisotropy of turbulent magneto-convection, kinetic and current and cross helicities.

The sunspot penumbra is a transition zone between the strong vertical magnetic field area (sunspot umbra) and the quiet Sun. The penumbra has a fine filamentary structure that is characterized by magnetic field lines inclined toward the surface. Numerical simulations of solar convection in inclined magnetic field regions have provided an explanation of the filamentary structure and the Evershed outflow in the penumbra. In this paper, we use radiative MHD simulations to investigate the influence of the magnetic field inclination on the power spectrum of vertical velocity oscillations. The results reveal a strong shift of the resonance mode peaks to higher frequencies in the case of a highly inclined magnetic field. The frequency shift for the inclined field is significantly greater than that in vertical field regions of similar strength. This is consistent with the behavior of fast MHD waves.

Five-minutes oscillations is one of the basic properties of solar convection. Observations show mixture of a large number of acoustic wave fronts propagating from their sources. We investigate the process of acoustic waves excitation from the point of view of individual events, by using realistic 3D radiative hydrodynamic simulation of the quiet Sun. The results show that the excitation events are related to dynamics vortex tubes (or swirls) in the intergranular lanes. These whirlpool-like flows are characterized by very strong horizontal velocities (7 - 11 km/s) and downflows (~ 7 km/s), and are accompanied by strong decreases of the temperature, density and pressure at the surface and in a ~ 0.5-1 Mm deep layer below the surface. High-speed whirlpool flows can attract and capture other vortices. According to our simulation results, the processes of the vortex interaction, such as vortex annihilation, can cause the excitation of acoustic waves.

Accurate subgrid-scale turbulence models are needed to perform realistic numerical magnetohydrodynamic (MHD) simulations of the subsurface flows of the Sun. To perform large-eddy simulations (LES) of turbulent MHD flows, three unknown terms have to be modeled. As a first step, this work proposes to use a priori tests to measure the accuracy of various models proposed to predict the SGS term appearing in the transport equation of the filtered magnetic field. It is proposed to evaluate the SGS model accuracy in term of "structural" and "functional" performance, i.e. the model capacity to locally approximate the unknown term and to reproduce its energetic action, respectively. From our tests, it appears that a mixed model based on the scale-similarity model has better performance.

One of the puzzling features of solar magnetism is formation of long-living compact magnetic structures; such as sunspots and pores, in the highly turbulent upper layer of the solar convective zone. We use realistic radiative 3D MHD simulations to investigate the interaction between magnetic field and turbulent convection. In the simulations, a weak vertical uniform magnetic field is imposed in a region of fully developed granular convection; and the total magnetic flux through the top and bottom boundaries is kept constant. The simulation results reveal a process of spontaneous formation of stable magnetic structures, which may be a key to understanding of the magnetic self-organization on the Sun and formation of pores and sunspots. This process consists of two basic steps: 1) formation of small-scale filamentary magnetic structures associated with concentrations of vorticity and whirlpool-type motions, and 2) merging of these structures due to the vortex attraction, caused by converging downdrafts around magnetic concentration below the surface. In the resulting large-scale structure maintained by the converging plasma motions, the magnetic field strength reaches ~1.5 kG at the surface and ~6 kG in the interior; and the surface structure resembles solar pores. The magnetic structure remains stable for the whole simulation run of several hours with no sign of decay.

Recent spectro-polarimetric observations of a sunspot showed the formation of bipolar magnetic patches in the mid penumbra and their propagation toward the outer penumbral boundary. The observations were interpreted as being caused by sea-serpent magnetic fields near the solar surface (Sainz Dalda & Bellot Rubio 2008). In this Letter, we develop a 3D radiative MHD numerical model to explain the sea-serpent structure and the wave-like behavior of the penumbral magnetic field lines. The simulations reproduce the observed behavior, suggesting that the sea-serpent phenomenon is a consequence of magnetoconvection in a strongly inclined magnetic field. It involves several physical processes: filamentary structurization, high-speed overturning convective motions in strong, almost horizontal magnetic fields with partially frozen field lines, and traveling convective waves. The results demonstrate a correlation of the bipolar magnetic patches with high-speed Evershed downflows in the penumbra. This is the first time that a 3D numerical model of the penumbra results in downward directed magnetic fields, an essential ingredient of sunspot penumbrae that has eluded explanation until now.

Discovered in 1909 the Evershed effect represents strong mass outflows in sunspot penumbra, where the magnetic field of sunspots is filamentary and almost horizontal. These flows play important role in sunspots and have been studied in detail using large ground-based and space telescopes, but the basic understanding of its mechanism is still missing. We present results of realistic numerical simulations of the Sun's subsurface dynamics, and argue that the key mechanism of this effect is in non-linear magnetoconvection that has properties of traveling waves in the presence of strong, highly inclined magnetic field. The simulations reproduce many observed features of the Evershed effect, including the high-speed "Evershed clouds", the filamentary structure of the flows, and the non-stationary quasi-periodic behavior. The results provide a synergy of previous theoretical models and lead to an interesting prediction of a large-scale organization of the outflows.

We have developed 3D, compressible, non-linear radiative MHD simulations to study the influence of the magnetic field of various strength and geometry on the turbulent convective cells and on the excitation mechanisms of the acoustic oscillations. The results reveal substantial changes of the granulation structure with increased magnetic field, and a frequency-dependent reduction in the oscillation power. These simulation results reproduce the enhanced high-frequency acoustic emission observed at the boundaries of active region ("acoustic halo" phenomenon). In the presence of inclined magnetic field the solar convection develops filamentary structure with flows concentrated along magnetic filaments, and also exhibits behavior of running magnetoconvective waves, resembling recent observations of the sunspot penumbra dynamics from Hinode/SOT.

The goal of this research is to investigate how magnetic field affects the dynamics of granular convection and excitation of solar oscillations by means of realistic numerical simulations. We have used a 3D, compressible, non-linear radiative magnetohydrodynamics code developed at the NASA Ames Research Center. This code takes into account several physical phenomena: compressible fluid flow in a highly stratified medium, sub-grid scale turbulence models, radiative energy transfer between the fluid elements, and a real-gas equation of state. We have studied the influence of the magnetic field of various strength on the convective cells and on the excitation mechanisms of the acoustic oscillations by calculating spectral properties of the convective motions and oscillations. The results reveal substantial changes of the granulation structure with increased magnetic field, and a frequency-dependent reduction in the oscillation power in a good agreement with solar observations. These simulations suggest that the enhanced high-frequency acoustic emission at the boundaries of active region ("acoustic halo" phenomenon) is caused by the changes of the spatial-temporal spectrum of the turbulent convection in magnetic field, resulting in turbulent motions of smaller scales and higher frequencies than in quiet Sun regions.

The goal of this research is to investigate how well various turbulence models can describe physical properties of the upper convective boundary layer of the Sun. An accurate modeling of the turbulence motions is necessary for understanding the excitation mechanisms of solar oscillation modes. We have carried out realistic numerical simulations using several different physical Large Eddy Simulation (LES) models (Hyperviscosity approach, Smagorinsky, and dynamic models) to investigate how the differences in turbulence modeling affect the damping and excitation of the oscillations and their spectral properties and compare with observations. We have first calculated the oscillation power spectra of radial and non-radial modes supported by the computational box with the different turbulence models. Then we have calculated the work integral input to the modes to estimate the influence of the turbulence model on the depth and strength of the oscillation sources. We have compared these results with previous studies and with the observed properties of solar oscillations. We find that the dynamic turbulence model provides the best agreement with the helioseismic observations.