Nonlinear anelastic modal theory for solar convection

Solar envelope models are developed using single-mode anelastic equations as a description of turbulent convection which provide estimates for the variation with depth of the largest convective cellular flows, with horizontal sizes comparable to the total depth of the convection zone. These models can be used to describe compressible motions occurring over many density scale heights. Single-mode anelastic solutions are obtained for a solar envelope whose mean stratification is nearly adiabatic over most of its vertical extent because of the enthalpy flux explicitly carried by the big cell, while a subgrid scale representation of turbulent heat transport is incorporated into the treatment near the surface. It is shown that the single-mode equations allow two solutions for the same horizontal wavelength which are distinguished by the sense of the vertical velocity at the center of the three-dimensional cell. It is found that the upward directed flow experiences large pressure effects which can modify the density fluctuations so that the sense of the buoyancy force is changed, with buoyancy braking actually achieved near the top of the convection zone. It is suggested that such dynamical processes may explain why the amplitudes of flows related to the largest scales of convection are so weak in the solar atmosphere.