Radiative transfer in cylindrical threads with incident radiation. VII. Multi-thread models

Aims. Our aim is to improve on previous radiative transfer calculations in illuminated cylindrical threads in order to better understand the physical conditions in cool solar chromospheric and coronal structures commonly observed in hydrogen and helium lines. Methods. We solve the radiative transfer and statistical equilibrium equations in a two-dimensional cross-section of a cylindrical structure oriented horizontally and lying above the solar surface. The cylinder is filled with a mixture of hydrogen and helium, and is illuminated at a given altitude from the solar disc. We construct simple models made from a single thread, or from an ensemble of several threads along the line of sight. This first use of 2D multi-thread fine structure modelling combining hydrogen and helium radiative transfer allows us to compute synthetic emergent spectra from cylindrical structures and to study the effect of line-of-sight integration of an ensemble of threads under a range of physical conditions. We analyse the effects of variations in temperature distribution and in gas pressure.We consider the effect of multi-thread structures within a given field of view and the effect of peculiar velocities between the structures in a multi-thread model. These new models are compared to the single thread model, and tested with varying parameters. Results. The presence of a temperature gradient, with temperature increasing towards the edge of the cylindrical thread, reduces the relative importance of the incident radiation coming from the solar disc on the emergent intensities of most hydrogen and helium lines. We also find that when assuming randomly displaced threads in a given field of view, the integrated intensities of optically thick and thin transitions behave considerably differently. In optically thin lines, the emergent intensity increases proportionally with the number of threads, and the spatial variation of the intensity becomes increasingly homogeneous. Optically thick lines however saturate after only a few threads. As a consequence, the spatial variation of the intensity retains much similarity with that of the first few threads. The multi-thread model produces complex line profiles with significant asymmetries if randomly generated line-of-sight velocities are added for each thread. Conclusions. These new computations show for the first time the effect of integrating the radiation emitted in H and He lines by several cylindrical threads static or moving along the line of sight. They can be used to interpret high-spatial and spectral resolutions of cylindrical structures found in the solar atmosphere, such as cool coronal loops or prominence threads

Earth's Inner Core dynamics induced by the Lorentz force

Seismic studies indicate that the Earth's inner core has a complex structure and exhibits a strong elastic anisotropy with a cylindrical symmetry. Among the various models which have been proposed to explain this anisotropy, one class of models considers the effect of the Lorentz force associated with the magnetic field diffused within the inner core. In this paper we extend previous studies and use analytical calculations and numerical simulations to predict the geometry and strength of the flow induced by the poloidal component of the Lorentz force in a neutrally or stably stratified growing inner core, exploring also the effect of different types of boundary conditions at the inner core boundary (ICB). Unlike previous studies, we show that the boundary condition that is most likely to produce a significant deformation and seismic anisotropy is impermeable, with negligible radial flow through the boundary. Exact analytical solutions are found in the case of a negligible effect of buoyancy forces in the inner core (neutral stratification), while numerical simulations are used to investigate the case of stable stratification. In this situation, the flow induced by the Lorentz force is found to be localized in a shear layer below the ICB, which thickness depends on the strength of the stratification, but not on the magnetic field strength. We obtain scaling laws for the thickness of this layer, as well as for the flow velocity and strain rate in this shear layer as a function of the control parameters, which include the magnitude of the magnetic field, the strength of the density stratification, the viscosity of the inner core, and the growth rate of the inner core. We find that the resulting strain rate is probably too small to produce significant texturing unless the inner core viscosity is smaller than about $10^{12}$ Pa.s.Comment: submitted to Geophysical Journal Internationa

Solar prominence modelling and plasma diagnostics at ALMA wavelengths

Our aim is to test potential solar prominence plasma diagnostics as obtained with the new solar capability of the Atacama Large Millimeter / submillimeter Array (ALMA). We investigate the thermal and plasma diagnostic potential of ALMA for solar prominences through the computation of brightness temperatures at ALMA wavelengths. The brightness temperature, for a chosen line of sight, is calculated using densities of hydrogen and helium obtained from a radiative transfer code under non local thermodynamic equilibrium (NLTE) conditions, as well as the input internal parameters of the prominence model in consideration. Two distinct sets of prominence models were used: isothermal-isobaric fine-structure threads, and large-scale structures with radially increasing temperature distributions representing the prominence-to-corona transition region. We compute brightness temperatures over the range of wavelengths in which ALMA is capable of observing (0.32 - 9.6mm), however we particularly focus on the bands available to solar observers in ALMA cycles 4 and 5, namely 2.6 - 3.6mm (Band 3) and 1.1 - 1.4mm (Band 6). We show how the computed brightness temperatures and optical thicknesses in our models vary with the plasma parameters (temperature and pressure) and the wavelength of observation. We then study how ALMA observables such as the ratio of brightness temperatures at two frequencies can be used to estimate the optical thickness and the emission measure for isothermal and non-isothermal prominences. From this study we conclude that, for both sets of models, ALMA presents a strong thermal diagnostic capability, provided that the interpretation of observations is supported by the use of non-LTE simulation results.Comment: Submitted to Solar Physic

Geophysically consistent values of the perovskite to post-perovskite transition Clapeyron slope

International audienceThe doublecrossing hypothesis posits that post-perovskite bearing rock in Earth's D 00 layer exists as a layer above the core-mantle boundary bounded above and below by intersections between a curved thermal boundary layer geotherm and a relatively steep phase boundary. Increasing seismic evidence for the existence of pairs of discontinuities predicted to occur at the top and bottom of this layer motivates an examination of the consistency of this model with mineral physics constraints for the Clapeyron slope of this phase transition. Using independent constraints for a lower bound on temperature in Earth's deep mantle and the temperature of Earth's inner core boundary, we show that a post-perovskite doublecrossing is inconsistent with plausible core temperatures for a Clapeyron slope less than about 7 MPa/K, with the higher range of experimental values yielding better agreement with recent estimates of the melting temperature of Earth's core

Modeling of the hydrogen Lyman lines in solar flares

The hydrogen Lyman lines (91.2 nm < λ < 121.6 nm) are significant contributors to the radiative losses of the solar chromosphere, and they are enhanced during flares. We have shown previously that the Lyman lines observed by the Extreme Ultraviolet Variability instrument onboard the Solar Dynamics Observatory exhibit Doppler motions equivalent to speeds on the order of 30 km s−1. However, contrary to expectations, both redshifts and blueshifts were present and no dominant flow direction was observed. To understand the formation of the Lyman lines, particularly their Doppler motions, we have used the radiative hydrodynamic code, RADYN, along with the radiative transfer code, RH, to simulate the evolution of the flaring chromosphere and the response of the Lyman lines during solar flares. We find that upflows in the simulated atmospheres lead to blueshifts in the line cores, which exhibit central reversals. We then model the effects of the instrument on the profiles, using the Extreme Ultraviolet Variability Experiment (EVE) instrument's properties. What may be interpreted as downflows (redshifted emission) in the lines, after they have been convolved with the instrumental line profile, may not necessarily correspond to actual downflows. Dynamic features in the atmosphere can introduce complex features in the line profiles that will not be detected by instruments with the spectral resolution of EVE, but which leave more of a signature at the resolution of the Spectral Investigation of the Coronal Environment instrument onboard the Solar Orbiter

Toward detailed prominence seismology - I. Computing accurate 2.5D magnetohydrodynamic equilibria

Context. Prominence seismology exploits our knowledge of the linear eigenoscillations for representative magnetohydro- dynamic models of filaments. To date, highly idealized models for prominences have been used, especially with respect to the overall magnetic configurations. Aims. We initiate a more systematic survey of filament wave modes, where we consider full multi-dimensional models with twisted magnetic fields representative of the surrounding magnetic flux rope. This requires the ability to compute accurate 2.5 dimensional magnetohydrodynamic equilibria that balance Lorentz forces, gravity, and pressure gradients, while containing density enhancements (static or in motion). Methods. The governing extended Grad-Shafranov equation is discussed, along with an analytic prediction for circular flux ropes for the Shafranov shift of the central magnetic axis due to gravity. Numerical equilibria are computed with a finite element-based code, demonstrating fourth order accuracy on an explicitly known, non-trivial test case. Results. The code is then used to construct more realistic prominence equilibria, for all three possible choices of a free flux-function. We quantify the influence of gravity, and generate cool condensations in hot cavities, as well as multi- layered prominences. Conclusions. The internal flux rope equilibria computed here have the prerequisite numerical accuracy to allow a yet more advanced analysis of the complete spectrum of linear magnetohydrodynamic perturbations, as will be demonstrated in the companion paper.Comment: Accepted by Astronomy & Astrophysics, 15 pages, 15 figure

Automatic detection of limb prominences in 304 A EUV images

A new algorithm for automatic detection of prominences on the solar limb in 304 A EUV images is presented, and results of its application to SOHO/EIT data discussed. The detection is based on the method of moments combined with a classifier analysis aimed at discriminating between limb prominences, active regions, and the quiet corona. This classifier analysis is based on a Support Vector Machine (SVM). Using a set of 12 moments of the radial intensity profiles, the algorithm performs well in discriminating between the above three categories of limb structures, with a misclassification rate of 7%. Pixels detected as belonging to a prominence are then used as starting point to reconstruct the whole prominence by morphological image processing techniques. It is planned that a catalogue of limb prominences identified in SOHO and STEREO data using this method will be made publicly available to the scientific community

On acceleration of Krylov-subspace-based Newton and Arnoldi iterations for incompressible CFD: replacing time steppers and generation of initial guess

We propose two techniques aimed at improving the convergence rate of steady state and eigenvalue solvers preconditioned by the inverse Stokes operator and realized via time-stepping. First, we suggest a generalization of the Stokes operator so that the resulting preconditioner operator depends on several parameters and whose action preserves zero divergence and boundary conditions. The parameters can be tuned for each problem to speed up the convergence of a Krylov-subspace-based linear algebra solver. This operator can be inverted by the Uzawa-like algorithm, and does not need a time-stepping. Second, we propose to generate an initial guess of steady flow, leading eigenvalue and eigenvector using orthogonal projection on a divergence-free basis satisfying all boundary conditions. The approach, including the two proposed techniques, is illustrated on the solution of the linear stability problem for laterally heated square and cubic cavities

Kinematics and helicity evolution of a loop-like eruptive prominence

We aim at investigating the morphology, kinematic and helicity evolution of a loop-like prominence during its eruption. We use multi-instrument observations from AIA/SDO, EUVI/STEREO and LASCO/SoHO. The kinematic, morphological, geometrical, and helicity evolution of a loop-like eruptive prominence are studied in the context of the magnetic flux rope model of solar prominences. The prominence eruption evolved as a height expanding twisted loop with both legs anchored in the chromosphere of a plage area. The eruption process consists of a prominence activation, acceleration, and a phase of constant velocity. The prominence body was composed of left-hand (counter-clockwise) twisted threads around the main prominence axis. The twist during the eruption was estimated at 6pi (3 turns). The prominence reached a maximum height of 526 Mm before contracting to its primary location and partially reformed in the same place two days after the eruption. This ejection, however, triggered a CME seen in LASCO C2. The prominence was located in the northern periphery of the CME magnetic field configuration and, therefore, the background magnetic field was asymmetric with respect to the filament position. The physical conditions of the falling plasma blobs were analysed with respect to the prominence kinematics. The same sign of the prominence body twist and writhe, as well as the amount of twisting above the critical value of 2pi after the activation phase indicate that possibly conditions for kink instability were present. No signature of magnetic reconnection was observed anywhere in the prominence body and its surroundings. The filament/prominence descent following the eruption and its partial reformation at the same place two days later suggest a confined type of eruption. The asymmetric background magnetic field possibly played an important role in the failed eruption.Comment: 9 pages, 8 figures, in press in A&