Search

Using numerical simulations, we study the formation and dynamics of post-flare loops in a local region of the solar atmosphere. The MHD equations rule the post-flare structures' dynamic evolution, including space-dependent magnetic resistivity and highly anisotropic thermal conduction on a 2.5 D slice. We use an initial magnetic configuration consisting of a vertical current sheet, which helps trigger the magnetic reconnection process. Specifically, we study two scenarios, one with only resistivity and the second with resistivity plus thermal conduction. Numerical simulations show differences in the global morphology of the post-flare substructures in both cases. In particular, localized resistivity produces more substructure on the loops related to a Ritchmyer-Meshkov Instability (RMI). On the other hand, in the scenario with resistivity plus thermal conduction, the post-flare loops are smooth, and no apparent substructures develop. Besides, in the $z-$component of the current density for the Res+TC scenario, we observe the development of multiple small magnetic islands along the current sheet., Comment: 16 pages, 14 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

This paper explores the effects of numerical algorithms on global magnetohydrodynamics (MHD) simulations of solar wind (SW) in the inner heliosphere. To do so, we use sunRunner3D, a 3-D MHD model that employs the boundary conditions generated by CORHEL and the PLUTO code to compute the plasma properties of the SW with the ideal MHD approximation up to 1.1 AU in the inner heliosphere. Mainly, we define three different combinations of numerical algorithms based on their diffusion level. This diffusion level is related to the way of solving the MHD equations using the finite volume formulation, and, therefore, we set in terms of the divergence-free condition methods, Riemann solvers, variable reconstruction schemes, limiters, and time-steeping algorithms. According to the simulation results, we demonstrate that sunRunner3D reproduces global features of Corotating Interaction Regions (CIRs) observed by Earth-based spacecraft (OMNI) for a set of Carrington rotations that cover a period that lays in the late declining phase of solar cycle 24, independently of the numerical algorithms. Moreover, statistical analyses between models and in-situ measurements show reasonable agreement with the observations, and remarkably, the high diffusive method matches better with in-situ data than low diffusive methods., Comment: 13 pages, 6 figures. Accepted for publication in Revista Mexicana de F\'isica

We perform a series of numerical simulations to recreate small-scale two-fluid jets using the JOANNA code, considering the magnetohydrodynamics of two fluids (ions + electrons and neutral particles). We first excite the jets in a uniform magnetic field by using velocity pulse perturbations located at $y_{0}=$1.3, 1.5, and 1.8 Mm, considering the base of the photosphere at $y=0$ Mm. Then, we repeat the excitation of the jets in a magnetic field that mimics a flux tube. Mainly, the jets excited at the upper chromosphere ($y\sim1.8$ Mm) reach lower heights than those excited at the lower chromosphere ($y\sim1.3$ Mm); this is due to the higher initial vertical location because of the lesser amount of plasma dragging. In both scenarios, the dynamics of the neutral particles and ions show similar behavior; however, we can still identify some differences in the velocity drift, which in our simulations is of the order of $10^{-3}$ km s$^{-1}$ at the tips of the jets once they reached their maximum heights. Also, we estimate the heat generation due to the friction between ions and neutrals ($Q^{in}_{i,n}$), which is of the order of $0.002-0.06$ W m$^{-3}$; however it is small to contribute to the heating of the surroundings of the solar corona. The jets in the two magnetic environments do not show substantial differences other than a slight variation in the maximum heights reached, particularly in the uniform magnetic field scenario. Finally, the maximum heights reached by the three different jets are in the range of some morphological parameters corresponding to macrospicules, Type I spicules, and Type II spicules., Comment: 15 pages, 6 figures. Accepted for publication in Physics

In the present paper, we use Si IV 1393.755 \AA\ spectral line observed by the Interface Region Imaging Spectrograph (IRIS) in the quiet-Sun to determine physical nature of the solar transition region (TR) oscillations. We analyze the properties of these oscillations using wavelet tools (e.g., power, cross-power, coherence, and phase difference) along with the stringent noise model (i.e., power-law + constant). We estimate the period of the intensity and Doppler velocity oscillations at each chosen location in the quiet-Sun (QS) and quantify the distribution of the statistically significant power and associated periods in one bright and two dark regions. In the bright TR region, the mean periods in intensity and velocity are 7 min, and 8 min respectively. In the dark region, the mean periods in intensity and velocity are 7 min, and 5.4 min respectively. We also estimate the phase difference between the intensity and Doppler velocity oscillations at each location. The statistical distribution of phase difference is estimated, which peaks at -119\degree $\pm$ 13\degree, 33\degree $\pm$ 10\degree, 102\degree $\pm$ 10\degree\ in the bright region, while at -153\degree $\pm$ 13\degree, 6\degree $\pm$ 20\degree, 151\degree $\pm$ 10\degree\ in the dark region. The statistical distribution reveals that the oscillations are caused by propagating slow magnetoacoustic waves encountered with the TR. Some of these locations may also be associated with the standing slow waves. Even, in the given time domain, several locations exhibit presence of both propagating and standing oscillations at different frequencies., Comment: 16 pages, 11 figures

We study the formation and evolution of jets in the solar atmosphere using numerical simulations of partially ionized plasma. The two-fluid magnetohydrodynamic equations with ion+electron and neutral hydrogen components are used in two-dimensional (2D) Cartesian geometry. Numerical simulations show that a localized nonlinear Gaussian pulse of ion and neutral pressures initially launched from the magnetic null point of a potential arcade located below the transition region quickly develops into a shock due to the decrease of density with height. The shock propagates upwards into the solar corona and lifts the cold and dense chromospheric plasma behind in the form of a collimated jet with an inverted-Y shape. The inverted-Y shape of jets is connected with the topology of a magnetic null point. The pulse also excites a nonlinear wake in the chromosphere, which leads to quasi-periodic secondary shocks. The secondary shocks lift the chromospheric plasma upwards and create quasi-periodic jets in the lower corona. Ion and neutral fluids show generally similar behavior, but their relative velocity is higher near the upper part of jets, which leads to enhanced temperature or heating due to ion-neutral collisions. Simulations of jets with inverted-Y shape and their heating may explain the properties of some jets observed in the solar atmosphere., Comment: 14 pages, 11 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Using numerical simulations, we study the effects of thermal conduction and radiative cooling on the formation and evolution of solar jets with some macrospicules features. We initially assume that the solar atmosphere is rarely in equilibrium through energy imbalance. Therefore, we test whether the background flows resulting from an imbalance between thermal conduction and radiative cooling influence the jets' behavior. In this particular scenario, we trigger the formation of the jets by launching a vertical velocity pulse localized at the upper chromosphere for the following test cases: i) adiabatic case; ii) thermal conduction case; iii) radiative cooling case; iv) thermal conduction + radiative cooling case. According to the test results, the addition of the thermal conduction results in smaller and hotter jets than in the adiabatic case. On the other hand, the radiative cooling dissipates the jet after reaching the maximum height ($\approx$ 5.5 Mm), making it shorter and colder than in the adiabatic and thermal conduction cases. Besides, the flow generated by the radiative cooling is more substantial than the caused by thermal conduction. Despite the energy imbalance of the solar atmosphere background, the simulated jet shows morphological features of macrospicules. Furthermore, the velocity pulse steepens into a shock that propagates upward into a solar corona that maintains its initial temperature. The shocks generate the jets with a quasi-periodical behavior that follows a parabolic path on time-distance plots consistent with macrospicule jets' observed dynamics., Comment: 15 pages, 12 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Using numerical simulations, we study the effects of magnetic resistivity and thermal conductivity in the dynamics and properties of solar jets with characteristics of Type II spicules and cool coronal jets. The dynamic evolution of the jets is governed by the resistive MHD equations with thermal conduction along the magnetic field lines on a 2.5D slice. The magnetic field configuration consists of two symmetric neighboring loops with opposite polarity, used to support reconnection and followed by the plasma jet formation. In total 10 simulations were carried out with different values of resistivity and thermal conductivity, that produce jets with different morphological and thermal properties we quantify. We find that an increase in magnetic resistivity does not produce significant effects on the morphology, velocity and temperature of the jets. However, thermal conductivity affects both temperature and morphology of the jets. In particular, thermal conductivity causes jets to reach greater heights and increases the temperature of the jet-apex. Also, heat flux maps indicate the jet-apex and corona interchange energy more efficiently than the body of the jet. These results could potentially open a new avenue for plasma diagnostics in the Sun's atmosphere., Comment: 11 pages, 4 figures. Accepted for publication in the Astrophysical Journal

We study the process of magnetic field annihilation and reconnection in simulations of magnetised solar photosphere and chromosphere with magnetic fields of opposite polarities and constant numerical resistivity. Exact analytical solutions for reconnective annihilations are used to interpret the features of magnetic reconnection in simulations of flux cancellation in the solar atmosphere. We use MURaM high-resolution photospheric radiative magneto-convection simulations to demonstrate the presence of magnetic field reconnection consistent with the magnetic flux pile-up models. Also, a simulated data-driven chromospheric magneto-hydrodynamic simulation is used to demonstrate magnetic field and flow structures, which are similar to the ones theoretically predicted. Both simulations demonstrate flow and magnetic field structures roughly consistent with accelerated reconnection with magnetic flux pile-up. The presence of standard Sweet-Parker type reconnection is also demonstrated in stronger photospheric magnetic fields., Comment: 7 pages, 6 figures. Resubmitted after taking into account the referee's comments

We present a new numerical model called CAFE-PakalMPI with the capability to solve the equations of classical magnetohydrodynamic (MHD) and to obtain the multispecies whose ionization states are calculated through statistical equilibrium, using the approximation of non-local thermodynamic equilibrium (NLTE) in three dimensions with the multiprocessor environment. For this, we couple the Newtonian CAFE MHD code introduced in Gonz\'alez-Avil\'es et al. (2015) with PakalMPI presented in De la Luz et al. (2010). In this model, Newtonian CAFE solves the equations of ideal MHD under the effects of magnetic resistivity and heat transfer considering a fully ionized plasma. On the other hand, PakalMPI calculates the density of ionization states using the NLTE approximation for Hydrogen, electronic densities and H-, the other species are computed by the classical local thermodynamic equilibrium (LTE) approximation. The main purpose of the model focuses on the analysis of solar phenomena within the chromospheric region. As an application of the model, we study the stability of the C7 equilibrium atmospheric model with a constant magnetic field in a 3D environment. According to the results of the test, the C7 model remains stable in the low chromosphere, while in the range $z\in$[1.5,2.1] Mm we can observe the propagation of a wave produced by the changes in the ionization rate of H., Comment: 9 pages, 9 figures

Within the framework of resistive MHD, implementing the C7 equilibrium atmosphere model and a 3D potential magnetic field realistic configuration, we simulate the formation of a plasm jet with the morphology, upward velocity up to 130 km/s and timescale formation between 60 and 90 s after beginning of simulation, similar to those expected for Type II spicules. Initial results of this simulation were published in Paper (e.g., Gonz\'alez-Avil\'es et al. 2018) and present paper is devoted to the analysis of transverse displacements and rotational type motion of the jet. Our results suggest that 3D magnetic reconnection may be responsible for the formation of the jet in Paper (Gonz\'alez-Avil\'es et al. 2018). In this paper, by calculating times series of the velocity components $v_x$ and $v_y$ in different points near to the jet for various heights we find transverse oscillations in agreement with spicule observations. We also obtain a time-distance plot of the temperature in a cross-cut at the plane $x=$0.1 Mm and find significant transverse displacements of the jet. By analyzing temperature isosurfaces of $10^{4}$ K with the distribution of $v_x$, we find that if the line-of-sight (LOS) is approximately perpendicular to the jet axis then there is both motion towards and away from the observer across the width of the jet. This red-blue shift pattern of the jet is caused by rotational motion, initially clockwise and anti-clockwise afterwards, which could be interpreted as torsional motion. From a nearly vertical perspective of the jet the LOS velocity component shows a central blue-shift region surrounded by red-shifted plasma., Comment: 8 pages, 7 figures

Using simulated data-driven three-dimensional resistive MHD simulations of the solar atmosphere, we show that magnetic reconnection can be responsible of the formation of jets with characteristic of Type II spicules. For this, we numerically model the photosphere-corona region using the C7 equilibrium atmosphere model. The initial magnetic configuration is a 3D potential magnetic field, extrapolated up to the solar corona region from a dynamic realistic simulation of solar photospheric magnetoconvection model which is mimicking quiet-Sun. In this case we consider a uniform and constant value of the magnetic resistivity of $12.56 ~\Omega~{\rm m}$. We have found that formation of the jets depends on the Lorentz force, which helps to accelerate the plasma upwards. Analyzing various properties of the jet dynamics, we found that the jet structure shows Doppler shift near to regions with high vorticity. The morphology, upward velocity, covering a range up to 100 $\rm km$ $\rm s^{-1}$, and life-time of the estructure, bigger than 100 ${\rm s}$, are similar to those expected for Type II spicules., Comment: 11 pages, 12 figures; Accepted for publication in the Astrophysical Journal

Using numerical simulations, we show that jets with features of type II spicules and cold coronal jets corresponding to temperatures $10^{4}$ K can be formed due to magnetic reconnection in a scenario in presence of magnetic resistivity. For this we model the low chromosphere-corona region using the C7 equilibrium solar atmosphere model and assuming Resistive MHD rules the dynamics of the plasma. The magnetic filed configurations we analyze correspond to two neighboring loops with opposite polarity. The separation of the loops' feet determines the thickness of a current sheet that triggers a magnetic reconnection process, and the further formation of a high speed and sharp structure. We analyze the cases where the magnetic filed strength of the two loops is equal and different. In the first case, with a symmetric configuration the spicules raise vertically whereas in an asymmetric configuration the structure shows an inclination. With a number of simulations carried out under a 2.5D approach, we explore various properties of excited jets, namely, the morphology, inclination and velocity. The parameter space involves magnetic field strength between 20 and 40 G, and the resistivity is assumed to be uniform with a constant value of the order $10^{-2}\Omega\cdot m$, Comment: Submitted to ApJ

We present a new code designed to solve the equations of classical ideal magneto-hydrodynamics (MHD) in three dimensions, submitted to a constant gravitational field. The purpose of the code centers on the analysis of solar phenomena within the photosphere-corona region. We present 1D and 2D standard tests to demonstrate the quality of the numerical results obtained with our code. As solar tests we present the transverse oscillations of Alfvenic pulses in coronal loops using a 2.5D model, and as 3D tests we present the propagation of impulsively generated MHD-gravity waves and vortices in the solar atmosphere. The code is based on high-resolution shock-capturing methods, uses the HLLE flux formula combined with Minmod, MC and WENO5 reconstructors. The divergence free magnetic field constraint is controlled using the Flux Constrained Transport method., Comment: 17 pages, 48 png fiugres,3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We solve numerically the ideal MHD equations with an external gravitational field in 2D in order to study the effects of impulsively generated linear and non-linear Alfv\'en waves into isolated solar arcades and coronal funnels. We analyze the region containing the interface between the photosphere and the corona. The main interest is to study the possibility that Alfv\'en waves triggers the energy flux transfer toward the quiet solar corona and heat it, including the case that two consecutive waves can occur. We find that in the case of arcades, short or large, the transferred fluxes by Alfv\'en waves are sufficient to heat the quiet corona only during a small lapse of time and in a certain region. In the case of funnels the threshold is achieved only when the wave is faster than 10 km/s, which is extremely high. We conclude from our analysis, that Alfv\'en waves, even in the optimistic scenario of having two consecutive Alfv\'en wave pulses, cannot transport enough energy as to heat the quiet corona., Comment: 13 pages, 48 png figures. Accepted for publication in MNRAS

We present the dynamics of rotating Bose Condensate galactic dark matter halos, made of an ultralight spinless boson. We restrict to the case of adding axisymmetric rigid rotation to initially spherically symmetric structures and show there are three regimes: i) small angular momentum, that basically retains the drawbacks of spherically symmetric halos related to compactness and failure at explaining galactic RCs, ii) an intermediate range of values of angular momentum that allow the existence of long-lived structures with acceptable RC profiles, and iii) high angular momentum, in which the structure is dispersed away by rotation. We also present in detail the new code used to solve the Gross-Pitaevskii Poisson system of equations in three dimensions., Comment: 10 revtex pages, 30 eps figures

In this paper we show that spherically symmetric BEC dark matter halos, with the $\sin r/r$ density profile, that accurately fit galactic rotation curves and represent a potential solution to the cusp-core problem are unstable. We do this by introducing back the density profiles into the fully time-dependent Gross-Pitaevskii-Poisson system of equations. Using numerical methods to track the evolution of the system, we found that these galactic halos lose mass at an approximate rate of half of its mass in a time scale of dozens of Myr. We consider this time scale is enough as to consider these halos are unstable and unlikely to be formed. We provide some arguments to show that this behavior is general and discuss some other drawbacks of the model that restrict its viability., Comment: 11 pages, 8 eps figures. Accepted for publication in JCAP