This paper presents a three-dimensional simulation of chromospheric jets with twisted magnetic field lines. Detailed treatments of the photospheric radiative transfer and the equations of state allow ...us to model realistic thermal convection near the solar surface, which excites various MHD waves and produces chromospheric jets in the simulation. A tall chromospheric jet with a maximum height of 10-11 Mm and lifetime of 8-10 minutes is formed above a strong magnetic field concentration. The magnetic field lines are strongly entangled in the chromosphere, which helps the chromospheric jet to be driven by the Lorentz force. The jet exhibits oscillatory motion as a natural consequence of its generation mechanism. We also find that the produced chromospheric jet forms a cluster with a diameter of several Mm with finer strands. These results imply a close relationship between the simulated jet and solar spicules.
The 11-year solar magnetic cycle shows a high degree of coherence in spite of the turbulent nature of the solar convection zone. It has been found in recent high-resolution magnetohydrodynamics ...simulations that the maintenance of a large-scale coherent magnetic field is difficult with small viscosity and magnetic diffusivity (≲10¹² square centimenters per second). We reproduced previous findings that indicate a reduction of the energy in the large-scale magnetic field for lower diffusivities and demonstrate the recovery of the global-scale magnetic field using unprecedentedly high resolution. We found an efficient small-scale dynamo that suppresses small-scale flows, which mimics the properties of large diffusivity. As a result, the global-scale magnetic field is maintained even in the regime of small diffusivities—that is, large Reynolds numbers.
ABSTRACT We investigate small-scale dynamo action in the solar convection zone through a series of high-resolution MHD simulations in a local Cartesian domain with (solar radius) of horizontal extent ...and a radial extent from 0.715 to . The dependence of the solution on resolution and diffusivity is studied. For a grid spacing of less than 350 km, the rms magnetic field strength near the base of the convection zone reaches 95% of the equipartition field strength (i.e., magnetic and kinetic energy are comparable). For these solutions the Lorentz force feedback on the convection velocity is found to be significant. The velocity near the base of the convection zone is reduced to 50% of the hydrodynamic one. In spite of the significant decrease of the convection velocity, the reduction in the enthalpy flux is relatively small, since the magnetic field also suppresses the horizontal mixing of the entropy between up- and downflow regions. This effect increases the amplitude of the entropy perturbation and makes convective energy transport more efficient. We discuss potential implications of these results for solar global convection and dynamo simulations.
A new three‐dimensional high‐resolution numerical model to study equatorial plasma bubble (EPB) has been developed. The High‐Resolution Bubble (HIRB) model is developed in a magnetic dipole ...coordinate system for the equatorial and low‐latitude ionosphere with a spatial resolution of as fine as 1 km. Adopting a higher‐order numerical scheme than those used in the existing models, the HIRB model is capable of reproducing the bifurcation, pinching, and turbulent structures of EPB. From a seeding perturbation resembling large‐scale wave structure (LSWS), EPB grows nonlinearly from the crest of LSWS upwelling, bifurcates at the top of EPB, then becomes turbulent at the topside of the F region. One of the bifurcated EPB is pinched off from the primary EPB and stops growing after pinching. The narrow channel of EPB tends to have a wiggle due to the secondary instability along the wall of EPB. Because of the fringe field effect above and below the EPB, upward drifting low‐density plasma converges toward the F peak altitude, forming a narrow‐depleted channel, and diverges above the peak, forming a flattened top of the EPB. The flattened top which has a steep upward density gradient is so unstable that bifurcation can easily occur even from a very small thermal perturbation. A higher density region between the bifurcated EPB moves downward due to westward polarization electric field. The EPB is pinched off when it reaches the wall of the primary EPB. It is concluded that turbulent plume‐like irregularities can be spontaneously generated only from large‐scale perturbation at the bottomside F region.
Key Points
A new high‐resolution plasma bubble model is developed
Turbulent plasma bubble from large‐scale wave structure is reproduced
Bifurcation and pinching processes of plasma bubble are explained
It has recently been recognized that the convective velocities achieved in current solar convection simulations might be overestimated. The newly revealed effects of the prevailing small-scale ...magnetic field within the convection zone may offer possible solutions to this problem. The small-scale magnetic fields can reduce the convective amplitude of small-scale motions through the Lorentz-force feedback, which concurrently inhibits the turbulent mixing of entropy between upflows and downflows. As a result, the effective Prandtl number may exceed unity inside the solar convection zone. In this paper, we propose and numerically confirm a possible suppression mechanism of convective velocity in the effectively high-Prandtl number regime. If the effective horizontal thermal diffusivity decreases (the Prandtl number accordingly increases), the subadiabatic layer which is formed near the base of the convection zone by continuous depositions of low entropy transported by adiabatically downflowing plumes is enhanced and extended. The global convective amplitude in the high-Prandtl thermal convection is thus reduced, especially in the lower part of the convection zone via the change in the mean entropy profile, which becomes more subadiabatic near the base and less superadiabatic in the bulk.
ABSTRACT We investigate the effect of coronal temperature on the formation process of solar chromospheric jets using two-dimensional magnetohydrodynamic simulations of the region from the upper ...convection zone to the lower corona. We develop a new radiative magnetohydrodynamic code for the dynamic modeling of the solar atmosphere, employing an LTE equation of state, optically thick radiative loss in the photosphere, optically thin radiative loss in the chromosphere and the corona, and thermal conduction along the magnetic field lines. Many chromospheric jets are produced in the simulations by shock waves passing through the transition region. We find that these jets are projected farther outward when the coronal temperature is lower (similar to that in coronal holes) and shorter when the coronal temperature is higher (similar to that in active regions). When the coronal temperature is high, the deceleration of the chromospheric jets is consistent with the model in which deceleration is determined by the periodic chromospheric shock waves. However, when the coronal temperature is low, the gravitational deceleration becomes more important and the chromospheric jets approach ballistic motion.
ABSTRACT The main object of the paper is to present the condition of the nondiffusive part of the Reynolds stress for driving the double-cell structure of the solar meridional circulation, which has ...been revealed by recent helioseismic observations. By conducting a set of mean-field hydrodynamic simulations, we confirm for the first time that double-cell meridional circulation can be achieved along with the solar-like differential rotation when the Reynolds stress transports the angular momentum upward in the lower part and downward in the upper part of the convection zone. It is concluded that in a stationary state, the accumulated angular momentum via the Reynolds stress in the middle layer is advected to both the upper and lower parts of the convection zone by each of the two meridional circulation cells, respectively.
ABSTRACT Transverse magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere and may be responsible for generating the Sun's million-degree outer atmosphere. However, direct evidence of ...the dissipation process and heating from these waves remains elusive. Through advanced numerical simulations combined with appropriate forward modeling of a prominence flux tube, we provide the observational signatures of transverse MHD waves in prominence plasmas. We show that these signatures are characterized by a thread-like substructure, strong transverse dynamical coherence, an out-of-phase difference between plane-of-the-sky motions and line-of-sight velocities, and enhanced line broadening and heating around most of the flux tube. A complex combination between resonant absorption and Kelvin-Helmholtz instabilities (KHIs) takes place in which the KHI extracts the energy from the resonant layer and dissipates it through vortices and current sheets, which rapidly degenerate into turbulence. An inward enlargement of the boundary is produced in which the turbulent flows conserve the characteristic dynamics from the resonance, therefore guaranteeing detectability of the resonance imprints. We show that the features described in the accompanying paper through coordinated Hinode and Interface Region Imaging Spectrograph observations match the numerical results well.
Magnetohydrodynamic (MHD) waves permeate the solar atmosphere and constitute potential coronal heating agents. Yet, the waves detected so far may be but a small subset of the true existing wave ...power. Detection is limited by instrumental constraints but also by wave processes that localize the wave power in undetectable spatial scales. In this study, we conduct 3D MHD simulations and forward modeling of standing transverse MHD waves in coronal loops with uniform and non-uniform temperature variation in the perpendicular cross-section. The observed signatures are largely dominated by the combination of the Kelvin-Helmholtz instability (KHI), resonant absorption, and phase mixing. In the presence of a cross-loop temperature gradient, we find that emission lines sensitive to the loop core catch different signatures compared to those that are more sensitive to the loop boundary and the surrounding corona, leading to an out-of-phase intensity and Doppler velocity modulation produced by KHI mixing. In all of the considered models, common signatures include an intensity and loop width modulation at half the kink period, a fine strand-like structure, a characteristic arrow-shaped structure in the Doppler maps, and overall line broadening in time but particularly at the loop edges. For our model, most of these features can be captured with a spatial resolution of 0 33 and a spectral resolution of 25 km s−1, although we do obtain severe over-estimation of the line width. Resonant absorption leads to a significant decrease of the observed kinetic energy from Doppler motions over time, which is not recovered by a corresponding increase in the line width from phase mixing and KHI motions. We estimate this hidden wave energy to be a factor of 5-10 of the observed value.
Current analytical and numerical modeling suggest the existence of ubiquitous thin current sheets in the corona that could explain the observed heating requirements. On the other hand, new high ...resolution observations of the corona indicate that its magnetic field may tend to organize itself in fine strand-like structures of few hundred kilometers widths. The link between small structure in models and the observed widths of strand-like structure several orders of magnitude larger is still not clear. A popular theoretical scenario is the nanoflare model, in which each strand is the product of an ensemble of heating events. Here, we suggest an alternative mechanism for strand generation. Through forward modeling of three-dimensional MHD simulations we show that small amplitude transverse MHD waves can lead in a few periods time to strand-like structure in loops in EUV intensity images. Our model is based on previous numerical work showing that transverse MHD oscillations can lead to Kelvin-Helmholtz instabilities that deform the cross-sectional area of loops. While previous work has focused on large amplitude oscillations, here we show that the instability can occur even for low wave amplitudes for long and thin loops, matching those presently observed in the corona. We show that the vortices generated from the instability are velocity sheared regions with enhanced emissivity hosting current sheets. Strands result as a complex combination of the vortices and the line-of-sight angle, last for timescales of a period, and can be observed for spatial resolutions of a tenth of loop radius.