No Arabic abstract
This paper presents the results of the analysis of 3D simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation-runs are compared: one including both ambipolar diffusion and Hall effect; one including only ambipolar diffusion; and one without any of these two effects. The magnetic field is amplified from initial field to saturation level by the action of turbulent local dynamo. In each of these cases, we study 2 hours of simulated solar time after the local dynamo reaches the saturation regime. We analyze the power spectra of vorticity, of magnetic field fluctuations and of the different components of the magnetic Poynting flux responsible for the transport of vertical or horizontal perturbations. Our preliminary results show that the ambipolar diffusion produces a strong reduction of vorticity in the upper chromospheric layers and that it dissipates the vortical perturbations converting them into thermal energy. The Hall effect acts in the opposite way, strongly enhancing the vorticity. When the Hall effect is included, the magnetic field in the simulations becomes, on average, more vertical and long-lived flux tube-like structures are produced. We trace a single magnetic structure to study its evolution pattern and the magnetic field intensification, and their possible relation to the Hall effect.
The partial ionization of the solar plasma causes several nonideal effects such as the ambipolar diffusion, the Hall effect, and the Biermann battery effect. Here we report on the first three-dimensional realistic simulations of solar local dynamo where all three effects were taken into account. The simulations started with a snapshot of already saturated battery-seeded dynamo, where two new series were developed: one with solely ambipolar diffusion and another one also taking into account the Hall term in the generalized Ohms law. The simulations were then run for about 4 hours of solar time to reach the stationary regime and improve the statistics. In parallel, a purely MHD dynamo simulation was also run for the same amount of time. The simulations are compared in a statistical way. The results show that, with the inclusion of the ambipolar diffusion, the amplitudes of the incompressible perturbations related to Alfven waves are reduced, and the Poynting flux is absorbed, with a frequency dependence. The Hall effect causes the opposite action: significant excess of incompressible perturbations is generated and an excess of the Poynting flux is observed in the chromospheric layers. The model with ambipolar diffusion shows, on average, sharper current sheets and slightly more abundant fast magneto-acoustic shocks in the chromosphere. The model with the Hall effect has higher temperatures at the lower chromosphere and stronger and more vertical magnetic field concentrations all over the chromosphere. The study of high-frequency waves reveals that significant power of incompressible perturbations is associated with areas with intense and more vertical magnetic fields and larger temperatures. We find a positive correlation between the magnitude of the ambipolar heating and the temperature increase at the same location after a characteristic time of 10^2 sec.
Over the last decades, realistic 3D radiative-MHD simulations have become the dominant theoretical tool for understanding the complex interactions between the plasma and the magnetic field on the Sun. Most of such simulations are based on approximations of magnetohydrodynamics, without directly considering the consequences of the very low degree of ionization of the solar plasma in the photosphere and bottom chromosphere. The presence of large amount of neutrals leads to a partial decoupling of the plasma and the magnetic field. As a consequence of that, a series of non-ideal effects (ambipolar diffusion, Hall effect and battery effect) arises. The ambipolar effect is the dominant one in the solar chromosphere. Here we report on the first three-dimensional realistic simulations of magneto-convection including ambipolar diffusion and battery effects. The simulations are done using the newly developed Mancha3D code. Our results reveal that ambipolar diffusion causes measurable effects on the amplitudes of waves excited by convection in the simulations, on the absorption of Poynting flux and heating and on the formation of chromospheric structures. We provide a low limit on the chromospheric temperature increase due to the ambipolar effect using the simulations with battery-excited dynamo fields.
Simulating the long-term evolution of temperature and magnetic fields in neutron stars is a major effort in astrophysics, having significant impact in several topics. A detailed evolutionary model requires, at the same time, the numerical solution of the heat diffusion equation, the use of appropriate numerical methods to control non-linear terms in the induction equation, and the local calculation of realistic microphysics coefficients. Here we present the latest extension of the magneto-thermal 2D code in which we have coupled the crustal evolution to the core evolution, including ambipolar diffusion. It has also gained in modularity, accuracy, and efficiency. We revise the most suitable numerical methods to accurately simulate magnetar-like magnetic fields, reproducing the Hall-driven magnetic discontinuities. From the point of view of computational performance, most of the load falls on the calculation of microphysics coefficients. To a lesser extent, the thermal evolution part is also computationally expensive because it requires large matrix
The chromosphere is a partially ionized layer of the solar atmosphere, the transition between the photosphere where the gas motion is determined by the gas pressure and the corona dominated by the magnetic field. We study the effect of partial ionization for 2D wave propagation in a gravitationally stratified, magnetized atmosphere with properties similar to the solar chromosphere. We adopt an oblique uniform magnetic field in the plane of propagation with strength suitable for a quiet sun region. The theoretical model used is a single fluid magnetohydrodynamic approximation, where ion-neutral interaction is modeled by the ambipolar diffusion term. Magnetic energy can be converted into internal energy through the dissipation of the electric current produced by the drift between ions and neutrals. We use numerical simulations where we continuously drive fast waves at the bottom of the atmosphere. The collisional coupling between ions and neutrals decreases with the decrease of the density and the ambipolar effect becomes important. Fast waves excited at the base of the atmosphere reach the equipartition layer and reflect or transmit as slow waves. While the waves propagate through the atmosphere and the density drops, the waves steepen into shocks. The main effect of ambipolar diffusion is damping of the waves. We find that for the parameters chosen in this work, the ambipolar diffusion affects the fast wave before it is reflected, with damping being more pronounced for waves which are launched in a direction perpendicular to the magnetic field. Slow waves are less affected by ambipolar effects. The damping increases for shorter periods and larger magnetic field strengths. Small scales produced by the nonlinear effects and the superposition of different types of waves created at the equipartition height are efficiently damped by ambipolar diffusion.
We present an investigation of the velocity fields in early to late M-type star hydrodynamic models, and we simulate their influence on FeH molecular line shapes. The M star model parameters range between log g of 3.0 - 5.0 and Teff of 2500 K and 4000 K. Our aim is to characterize the Teff- and log g -dependence of the velocity fields and express them in terms of micro- and macro-turbulent velocities in the one dimensional sense. We present also a direct comparison between 3D hydrodynamical velocity fields and 1D turbulent velocities. The velocity fields strongly affect the line shapes of FeH, and it is our goal to give a rough estimate for the log g and Teff parameter range in which 3D spectral synthesis is necessary and where 1D synthesis suffices. In order to calculate M-star structure models we employ the 3D radiative-hydrodynamics (RHD) code CO5BOLD. The spectral synthesis on these models is performed with the line synthesis code LINFOR3D. We describe the 3D velocity fields in terms of a Gaussian standard deviation and project them onto the line of sight to include geometrical and limb-darkening effects. The micro- and macro-turbulent velocities are determined with the Curve of Growth method and convolution with a Gaussian velocity profile, respectively. To characterize the log g and Teff dependence of FeH lines, the equivalent width, line width, and line depth are regarded. The velocity fields in M-stars strongly depend on log g and Teff. They become stronger with decreasing log g and increasing Teff.