Search

In sunspot umbrae, the core of some chromospheric lines exhibits periodic brightness enhancements known as umbral flashes. The consensus is that they are produced by the upward propagation of shock waves. This view has recently been challenged by the detection of downflowing umbral flashes and the confirmation of the existence of a resonant cavity above sunspots. We aim to determine waves' propagating or standing nature in the low umbral chromosphere and confirm or refute the existence of downflowing umbral flashes. Spectroscopic temporal series of Ca II 8542 \AA, Ca II H, and Halpha in a sunspot were acquired with the Swedish Solar Telescope. The Halpha velocity was inferred using bisectors. Simultaneous inversions of the Ca II 8542 \AA\ line and the Ca II H core were performed using the NICOLE code. The nature of the oscillations and insights into the resonant oscillatory pattern were determined by analyzing the phase shift between the velocity signals and examining the temporal evolution. Propagating waves in the low chromosphere are more common in regions with frequent umbral flashes, where the transition region is shifted upward, making resonant cavity signatures less noticeable. In contrast, areas with fewer umbral flashes show velocity fluctuations that align with standing oscillations. Evidence suggests dynamic changes in the location of velocity resonant nodes due to variations in transition region height. Downflowing profiles appear at the onset of some umbral flashes, but upflowing motion dominates during most of the flash. These downflowing flashes are more common in standing umbral flashes. We confirm the existence of a chromospheric resonant cavity above sunspot umbrae produced by wave reflections at the transition region. The oscillatory pattern depends on the transition region height, which exhibits spatial and temporal variations due to the impact of the waves., Comment: Accepted for publication in A&A

Solar partially ionized plasma is frequently modeled using single-fluid (1F) or two-fluid (2F) approaches. In the 1F case, charge-neutral interactions are often described through ambipolar diffusion, while the 2F model fully considers charge-neutral drifts. Here, we expand the definition of the ambipolar diffusion coefficient to include inelastic collisions (ion/rec) in two cases: a VAL3C 1D model and a 2F simulations of the Rayleigh-Taylor instability (RTI) in a solar prominence thread based on \cite{PopLukKho2021aa, PopLukKho2021ab}. On one side, we evaluate the relative importance of the inelastic contribution, compared to elastic and charge-exchange collisions. On the other side, we compare the contributions of ion/rec, thermal pressure, viscosity, and magnetic forces to the charge-neutral drift velocity of the turbulent flow of the RTI. Our analysis reveals that the contribution of inelastic collisions to the ambipolar diffusion coefficient is negligible across the chromosphere, allowing the classical definition of this coefficient to be safely used in 1F modeling. However, in the transition region, the contribution of inelastic collisions can become as significant as that of elastic collisions. Furthermore, we ascertain that the thermal pressure force predominantly influences the charge-neutral drifts in the RTI model, surpassing the impact of the magnetic force., Comment: Accepted

The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic processes in solar/stellar atmospheres. The code includes non-ideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high quality realistic simulations of magneto-convection, as well as idealized simulations of particular processes, such as wave propagation, instabilities or energetic events. The paper summarizes the equations and methods used in the Mancha3D code. It also describes its numerical stability and parallel performance and efficiency. The code is based on a finite difference discretization and memory-saving Runge-Kutta (RK) scheme. It handles non-ideal effects through super-time stepping and Hall diffusion schemes, and takes into account thermal conduction by solving an additional hyperbolic equation for the heat flux. The code is easily configurable to perform different kinds of simulations. Several examples of the code usage are given. It is demonstrated that splitting variables into equilibrium and perturbation parts is essential for simulations of wave propagation in a static background. A perfectly matched layer (PML) boundary condition built into the code greatly facilitates a non-reflective open boundary implementation. Spatial filtering is an important numerical remedy to eliminate grid-size perturbations enhancing the code stability. Parallel performance analysis reveals that the code is strongly memory bound, which is a natural consequence of the numerical techniques used, such as split variables and PML boundary conditions. Both strong and weak scalings show adequate performance up till several thousands of CPUs.

Several generalizations of the well-known fluid model of Braginskii (Rev. of Plasma Phys., 1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multi-species plasmas with arbitrary masses and temperatures, where all the fluid moments are described by their evolution equations. The 21-moment model contains two "heat flux vectors" (3rd and 5th-order moments) and two "viscosity-tensors" (2nd and 4th-order moments). The Braginskii model is then obtained as a particular case of a one ion-electron plasma with similar temperatures, with de-coupled heat fluxes and viscosity-tensors expressed in a quasi-static approximation. We provide all the numerical values of the Braginskii model in a fully analytic form (together with the 4th and 5th-order moments). For multi-species plasmas, the model makes calculation of transport coefficients straightforward. Formulation in fluid moments (instead of Hermite moments) is also suitable for implementation into existing numerical codes. It is emphasized that it is the quasi-static approximation which makes some Braginskii coefficients divergent in a weakly-collisional regime. Importantly, we show that the heat fluxes and viscosity-tensors are coupled even in the linear approximation, and that the fully contracted (scalar) perturbations of the 4th-order moment, which are accounted for in the 22-moment model, modify the energy exchange rates. We also provide several Appendices, which can be useful as a guide for deriving the Braginskii model with the moment method of Grad., Comment: Updated version. Title change from "Braginskii-type" to "Braginskii Type" to match the published version (accepted to ApJ Supplements, 29 Jan 2022). Bookmarks included, TOC is now at the beginning, numbering in Appendix I is fixed, few typos corrected