The 'standard solar' flare model collects all physical ingredients identified by multi-wavelength observations of our Sun: magnetic reconnection, fast particle acceleration and the resulting emission at various wavelengths, especially in soft to hard X-ray channels. Its cartoon representation is found throughout textbooks on solar and plasma astrophysics, and guides interpretations of unresolved energetic flaring events on other stars, accretion disks and jets. To date, a fully self-consistent model that reproduces the standard scenario in all its facets is lacking, since this requires the combination of a large scale, multi-dimensional magnetohydrodynamic (MHD) plasma description with a realistic fast electron treatment. Here, we demonstrate such a novel combination, where MHD combines with an analytic fast electron model, adjusted to handle time-evolving, reconnecting magnetic fields and particle trapping. This allows to study (1) the role of fast electron deposition in the triggering of chromospheric evaporation flows; (2) the physical mechanisms that generate various hard X-ray sources at chromospheric footpoints or looptops; and (3) the relationship between soft X-ray and hard X-ray fluxes throughout the entire flare loop evolution. For the first time, this self-consistent solar flare model demonstrates the observationally suggested relationship between flux swept out by the hard X-ray footpoint regions, and the actual reconnection rate at the X-point, which is a major unknown in flaring scenarios. We also demonstrate that a looptop hard X-ray source can result from fast electron trapping.

We report on the latest additions to our open-source, block-grid adaptive framework MPI-AMRVAC, which is a general toolkit for especially hyperbolic/parabolic partial differential equations (PDEs). Applications traditionally focused on shock-dominated, magnetized plasma dynamics described by either Newtonian or special relativistic (magneto)hydrodynamics, but its versatile design easily extends to different PDE systems. Here, we demonstrate applications covering any-dimensional scalar to system PDEs, with e.g. Korteweg-de Vries solutions generalizing early findings on soliton behaviour, shallow water applications in round or square pools, hydrodynamic convergence tests as well as challenging computational fluid and plasma dynamics applications. The recent addition of a parallel multigrid solver opens up new avenues where also elliptic constraints or stiff source terms play a central role. This is illustrated here by solving several multi-dimensional reaction-diffusion-type equations. We document the minimal requirements for adding a new physics module governed by any nonlinear PDE system, such that it can directly benefit from the code flexibility in combining various temporal and spatial discretisation schemes. Distributed through GitHub, MPI-AMRVAC can be used to perform 1D, 1.5D, 2D, 2.5D or 3D simulations in Cartesian, cylindrical or spherical coordinate systems, using parallel domain-decomposition, or exploiting fully dynamic block quadtree-octree grids.

Prominences show a surprising amount of fine structure and it is widely believed that their threads, as seen in H$\alpha$ observations, provide indirect information concerning magnetic field topology. We investigate the spontaneous emergence and evolution of fine structure in high-density condensations formed through the process of thermal instability under typical solar coronal conditions. Our study reveals intricate multidimensional processes that occur through in situ condensations in a representative coronal volume in a low plasma beta regime. We performed 2D and 3D numerical simulations of interacting slow magnetohydrodynamic (MHD) wave modes when all relevant non-adiabatic effects are included. We show that the interaction of multiple slow MHD wave modes in a regime unstable to the thermal mode leads to thermal instability. This initially forms pancake-like structures almost orthogonal to the local magnetic field, while low-pressure induced inflows of matter generate rebound shocks. This is succeeded by the rapid disruption of these pancake sheets through thin-shell instabilities evolving naturally from minute ram pressure imbalances. This eventually creates high-density blobs accompanied by thread-like features from shear flow effects. The further evolution of the blobs follows the magnetic field lines, such that a dynamical realignment with the background magnetic field appears. However, the emerging thread-like features are not at all field-aligned, implying only a very weak link between fine structure orientation and magnetic field topology which has far-reaching implications for field topology interpretations based on H$\alpha$ observations.

The formation of shocks within the solar atmosphere remains one of the few observable signatures of energy dissipation arising from the plethora of magnetohydrodynamic waves generated close to the solar surface. Active region observations offer exceptional views of wave behavior and its impact on the surrounding atmosphere. The stratified plasma gradients present in the lower solar atmosphere allow for the potential formation of many theorized shock phenomena. In this study, using chromospheric Ca II 854.2nm spectropolarimetric data of a large sunspot, we examine fluctuations in the plasma parameters in the aftermath of powerful shock events that demonstrate polarimetric reversals during their evolution. Modern inversion techniques are employed to uncover perturbations in the temperatures, line-of-sight velocities, and vector magnetic fields occurring across a range of optical depths synonymous with the shock formation. Classification of these non-linear signatures is carried out by comparing the observationally-derived slow, fast, and Alfv\'en shock solutions to the theoretical Rankine-Hugoniot relations. Employing over 200,000 independent measurements, we reveal that the Alfv\'en (intermediate) shock solution provides the closest match between theory and observations at optical depths of log(tau) = -4, consistent with a geometric height at the boundary between the upper photosphere and lower chromosphere. This work uncovers first-time evidence of the manifestation of chromospheric intermediate shocks in sunspot umbrae, providing a new method for the potential thermalization of wave energy in a range of magnetic structures, including pores, magnetic flux ropes, and magnetic bright points.

The late stellar evolutionary phases of low and intermediate-mass stars are strongly constrained by their mass-loss rates. The wind surrounding cool evolved stars frequently shows non-spherical features, thought to be due to an unseen companion orbiting the donor star. We study the morphology of the circumbinary envelope, in particular around oxygen-rich asymptotic giant branch (AGB) stars. We run a grid of 70 3D hydrodynamics simulations of a progressively accelerating wind propagating in the Roche potential formed by a mass-loosing evolved star in orbit with a main sequence companion. We resolve the flow structure both in the immediate vicinity of the secondary, where bow shocks, outflows and wind-captured disks form, and up to 40 orbital separations, where spiral arms, arcs and equatorial density enhancements develop. When the companion is deeply engulfed in the wind, the lower terminal wind speeds and more progressive wind acceleration around oxygen-rich AGB stars make them more prone than carbon-rich AGB stars to display more disturbed outflows, a disk-like structure around the companion and a wind concentrated in the orbital plane. In these configurations, a large fraction of the wind is captured by the companion which leads to a significant shrinking of the orbit over the mass-loss timescale, if the donor star is at least a few times more massive than its companion. Provided the companion has a mass of at least a tenth of the mass of the donor star, it can compress the wind in the orbital plane up to large distances. Our grid of models covers a wide scope of configurations function of the dust chemical content, the terminal wind speed relative to the orbital speed, the extension of the dust condensation region around the cool evolved star and the mass ratio. It provides a frame of reference to interpret high-resolution maps of the outflows surrounding cool evolved stars.

We review and discuss insights on ideal magnetohydrodynamic (MHD) instabilities that can play a role in destabilizing solar coronal flux rope structures. For single flux ropes, failed or actual eruptions may result from internal or external kink evolutions, or from torus unstable configurations. We highlight recent findings from 3D magnetic field reconstructions and simulations where kink and torus instabilities play a prominent role. For interacting current systems, we critically discuss different routes to coronal dynamics and global eruptions, due to current channel coalescence or to tilt-kink scenarios. These scenarios involve the presence of two nearby current channels and are clearly distinct from the popular kink or torus instability. Since the solar corona is pervaded with myriads of magnetic loops -- creating interacting flux ropes typified by parallel or antiparallel current channels as exemplified in various recent observational studies -- coalescence or tilt-kink evolutions must be very common for destabilizing adjacent flux rope systems. Since they also evolve on ideal MHD timescales, they may well drive many sympathetic eruptions witnessed in the solar corona. Moreover, they necessarily lead to thin current sheets that are liable to reconnection. We review findings from 2D and 3D MHD simulations for tilt and coalescence evolutions, as well as on particle acceleration aspects derived from computed charged particle motions in tilt-kink disruptions and coalescing flux ropes. The latter were recently studied in two-way coupled kinetic-fluid models, where the complete phase-space information of reconnection is incorporated.

The Kelvin-Helmholtz (KH) instability is commonly found in many astrophysical, laboratory, and space plasmas. It could mix plasma components of different properties and convert dynamic fluid energy from large scale structure to smaller ones. In this study, we combined the ground-based New Vacuum Solar Telescope (NVST) and the Solar Dynamic Observatories (SDO) / Atmospheric Imaging Assembly (AIA) to observe the plasma dynamics associated with active region 12673 on 09 September 2017. In this multi-temperature view, we identified three adjacent layers of plasma flowing at different speeds, and detected KH instabilities at their interfaces. We could unambiguously track a typical KH vortex and measure its motion. We found that the speed of this vortex suddenly tripled at a certain stage. This acceleration was synchronized with the enhancements in emission measure and average intensity of the 193 \AA{} data. We interpret this as evidence that KH instability triggers plasma heating. The intriguing feature in this event is that the KH instability observed in the NVST channel was nearly complementary to that in the AIA 193 \AA{}. Such a multi-thermal energy exchange process is easily overlooked in previous studies, as the cold plasma component is usually not visible in the extreme ultraviolet channels that are only sensitive to high temperature plasma emissions. Our finding indicates that embedded cold layers could interact with hot plasma as invisible matters. We speculate that this process could occur at a variety of length scales and could contribute to plasma heating.

The Kelvin-Helmholtz instability (KHI) is a known mechanism for penetration of solar wind matter into the magnetosphere. Using three-dimensional, resistive magnetohydrodynamic simulations, the double mid-latitude reconnection (DMLR) process was shown to efficiently exchange solar wind matter into the magnetosphere, through mixing and reconnection. Here, we compute test particle orbits through DMLR configurations. In the instantaneous electromagnetic fields, charged particle trajectories are integrated using the guiding centre approximation. The mechanisms involved in the electron particle orbits and their kinetic energy evolutions are studied in detail, to identify specific signatures of the DMLR through particle characteristics. The charged particle orbits are influenced mainly by magnetic curvature drifts. We identify complex, temporarily trapped, trajectories where the combined electric field and (reconnected) magnetic field variations realize local cavities where particles gain energy before escaping. By comparing the orbits in strongly deformed fields due to the KHI development, with the textbook mirror-drift orbits resulting from our initial configuration, we identify effects due to current sheets formed in the DMLR process. We do this in various representative stages during the DMLR development.

Recent advances in black hole astrophysics, particularly the first visual evidence of a supermassive black hole at the center of the galaxy M87 by the Event Horizon Telescope (EHT), and the detection of an orbiting "hot spot" nearby the event horizon of Sgr A* in the Galactic center by the Gravity Collaboration, require the development of novel numerical methods to understand the underlying plasma microphysics. Non-thermal emission related to such hot spots is conjectured to originate from plasmoids that form due to magnetic reconnection in thin current layers in the innermost accretion zone. Resistivity plays a crucial role in current sheet formation, magnetic reconnection, and plasmoid growth in black hole accretion disks and jets. We included resistivity in the three-dimensional general-relativistic magnetohydrodynamics (GRMHD) code BHAC and present the implementation of an Implicit-Explicit scheme to treat the stiff resistive source terms of the GRMHD equations. The algorithm is tested in combination with adaptive mesh refinement to resolve the resistive scales and a constrained transport method to keep the magnetic field solenoidal. Several novel methods for primitive variable recovery, a key part in relativistic magnetohydrodynamics codes, are presented and compared for accuracy, robustness, and efficiency. We propose a new inversion strategy that allows for resistive-GRMHD simulations of low gas-to-magnetic pressure ratio and highly magnetized regimes as applicable for black hole accretion disks, jets, and neutron star magnetospheres. We apply the new scheme to study the effect of resistivity on accreting black holes, accounting for dissipative effects as reconnection.

We conduct forward-modeling analysis based on our 2.5 dimensional magnetohydrodynamics (MHD) simulation of magnetic flux rope (MFR) formation and eruption driven by photospheric converging motion. The current sheet (CS) evolution during the MFR formation and eruption process in our MHD simulation can be divided into four stages. The current sheet (CS) evolution during the MFR formation and eruption process in our MHD simulation can be divided into four stages. The first stage shows the CS forming and gradually lengthening. Resistive instabilities that disrupt the CS mark the beginning of the second stage. Magnetic islands disappear in the third stage and reappear in the fourth stage. Synthetic images and light curves of the seven Solar Dynamics Observatory/Atmospheric Imaging Assembly (AIA) channels, i.e., 9.4 nm, 13.1 nm, 17.1 nm, 19.3 nm, 21.1 nm, 30.4 nm, and 33.5 nm, and the 3-25 keV thermal X-ray are obtained with forward-modeling analysis. The loop-top source and the coronal sources of the soft X-ray are reproduced in forward modeling. The light curves of the seven SDO/AIA channels start to rise once resistive instabilities develop. The light curve of the 3-25 keV thermal X-ray starts to go up when the reconnection rate reaches one of its peaks. Quasiperiodic pulsations (QPPs) appear twice in the SDO/AIA 17.1 nm, 21.1 nm, and 30.4 nm channels, corresponding to the period of chaotic (re)appearance and CS-guided displacements of the magnetic islands. QPPs appear once in the SDO/AIA 9.4 nm and 33.5 nm channels after the disruption of the CS by resistive instabilities and in the 19.3 nm channel when the chaotic motion of the magnetic islands reappears.

The Black Hole Accretion Code (BHAC) has recently been extended with the ability to evolve charged test particles according to the Lorentz force within resistive relativistic magnetohydrodynamics simulations. We apply this method to evolve particles in a reconnecting current sheet that forms due to the coalescence of two magnetic flux tubes in 2D Minkowski spacetime. This is the first analysis of charged test particle evolution in resistive relativistic magnetohydrodynamics simulations. The energy distributions of an ensemble of 100.000 electrons are analyzed, as well as the acceleration of particles in the plasmoids that form in the reconnection layer. The effect of the Lundquist number, magnetization, and plasma-$\beta$ on the particle energy distribution is explored for a range of astrophysically relevant parameters. We find that electrons accelerate to non-thermal energies in the thin current sheets in all cases. We find two separate acceleration regimes: An exponential increase of the Lorentz factor during the island coalescence where the acceleration depends linearly on the resistivity and a nonlinear phase with high variability. These results are relevant for determining energy distributions and acceleration sites obtaining radiation maps in large-scale magnetohydrodynamics simulations of black hole accretion disks and jets.

We apply the general relativistic resistive magnetohydrodynamics code {\tt BHAC} to perform a 2D study of the formation and evolution of a reconnection layer in between two merging magnetic flux tubes in Minkowski spacetime. Small-scale effects in the regime of low resistivity most relevant for dilute astrophysical plasmas are resolved with very high accuracy due to the extreme resolutions obtained with adaptive mesh refinement. Numerical convergence in the highly nonlinear plasmoid-dominated regime is confirmed for a sweep of resolutions. We employ both uniform resistivity and non-uniform resistivity based on the local, instantaneous current density. For uniform resistivity we find Sweet-Parker reconnection, from $\eta = 10^{-2}$ down to $\eta = 10^{-4}$, for a reference case of magnetisation $\sigma = 3.33$ and plasma-$\beta = 0.1$. {For uniform resistivity $\eta=5\times10^{-5}$ the tearing mode is recovered, resulting in the formation of secondary plasmoids. The plasmoid instability enhances the reconnection rate to $v_{\rm rec} \sim 0.03c$ compared to $v_{\rm rec} \sim 0.01c$ for $\eta=10^{-4}$.} For non-uniform resistivity with a base level $\eta_0 = 10^{-4}$ and an enhanced current-dependent resistivity in the current sheet, we find an increased reconnection rate of $v_{\rm rec} \sim 0.1c$. The influence of the magnetisation $\sigma$ and the plasma-$\beta$ is analysed for cases with uniform resistivity $\eta=5\times10^{-5}$ and $\eta=10^{-4}$ in a range $0.5 \leq \sigma \leq 10$ and $0.01 \leq \beta \leq 1$ in regimes that are applicable for black hole accretion disks and jets. The plasmoid instability is triggered for Lundquist numbers larger than a critical value of $S_{\rm c} \approx 8000$.

We compare relativistic particle integrators commonly used in plasma physics showing several test cases relevant for astrophysics. Three explicit particle pushers are considered, namely the Boris, Vay, and Higuera-Cary schemes. We also present a new relativistic fully implicit particle integrator that is energy conserving. Furthermore, a method based on the relativistic guiding center approximation is included. The algorithms are described such that they can be readily implemented in magnetohydrodynamics codes or Particle-in-Cell codes. Our comparison focuses on the strengths and key features of the particle integrators. We test the conservation of invariants of motion, and the accuracy of particle drift dynamics in highly relativistic, mildly relativistic, and non-relativistic settings. The methods are compared in idealized test cases, i.e., without considering feedback on the electrodynamic fields, collisions, pair creation, or radiation. The test cases include uniform electric and magnetic fields, $\mathbf{E}\times\mathbf{B}$-fields, force-free fields, and setups relevant for high-energy astrophysics, e.g., a magnetic mirror, a magnetic dipole, and a magnetic null. These tests have direct relevance for particle acceleration in shocks and in magnetic reconnection.

We analyze particle acceleration in explosive reconnection events in magnetically dominated proton-electron plasmas. Reconnection is driven by large-scale magnetic stresses in interacting current-carrying flux tubes. Our model relies on development of current-driven instabilities on macroscopic scales. These tilt-kink instabilities develop in an initially force-free equilibrium of repelling current channels. Using MHD methods we study a 3D model of repelling and interacting flux tubes in which we simultaneously evolve test particles, guided by electromagnetic fields obtained from MHD. We identify two stages of particle acceleration; Initially particles accelerate in the current channels, after which the flux ropes start tilting and kinking and particles accelerate due to reconnection processes in the plasma. The explosive stage of reconnection produces non-thermal energy distributions with slopes that depend on plasma resistivity and the initial particle velocity. We also discuss the influence of the length of the flux ropes on particle acceleration and energy distributions. This study extends previous 2.5D results to 3D setups, providing all ingredients needed to model realistic scenarios like solar flares, black hole flares and particle acceleration in pulsar wind nebulae: formation of strong resistive electric fields, explosive reconnection and non-thermal particle distributions. By assuming initial energy equipartition between electrons and protons, applying low resistivity in accordance with solar corona conditions and limiting the flux rope length to a fraction of a solar radius we obtain realistic energy distributions for solar flares with non-thermal power law tails and maximum electron energies up to 11 MeV and maximum proton energies up to 1 GeV.

The combination of magnetohydrodynamic (MHD) simulation and multi-wavelength observations is an effective way to study mechanisms of magnetic flux rope eruption. We develop a data-driven MHD model using the zero-$\beta$ approximation. The initial condition is provided by nonlinear force-free field derived by the magneto-frictional method based on vector magnetic field observed by the Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics Observatory (SDO). The bottom boundary uses observed time series of the vector magnetic field and the vector velocity derived by the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM). We apply the data-driven model to active region 11123 observed from 06:00 UT on 2011 November 11 to about 2 hours later. The evolution of the magnetic field topology coincides with the flare ribbons observed in the 304 and 1600 \AA $ $ wavebands by the Atmospheric Imaging Assembly. The morphology, propagation path, and propagation range of the flux rope are comparable with the observations in 304 \AA . We also find that a data-constrained boundary condition, where the bottom boundary is fixed to the initial values, reproduces a similar simulation result. This model can reproduce the evolution of a magnetic flux rope in its dynamic eruptive phase.

Ultra-luminous X-ray sources (ULX) have so high X-ray luminosities that they were long thought to be accreting intermediate mass black holes. Yet, some ULX have been shown to display periodic modulations and coherent pulsations, suggestive of a neutron star in orbit around a stellar companion and accreting at super-Eddington rates. In this letter, we propose that the mass transfer in ULX could be qualitatively the same as in Supergiant X-ray binaries (SgXB), with a wind from the donor star highly beamed towards the compact object. Since the star does not fill its Roche lobe, this mass transfer mechanism known as "wind Roche lobe overflow" can remain stable even for large mass ratios. Based on realistic acceleration profiles derived from spectral observations and modeling of the stellar wind, we compute the bulk motion of the wind to evaluate the fraction of the stellar mass outflow captured by the compact object. We identify the orbital and stellar conditions for a SgXB to transfer mass at rates matching the expectations for ULX and show that the transition from SgXB to ULX luminosity levels is progressive. These results indicate that a high stellar Roche lobe filling factor is not necessary to funnel large quantities of material into the Roche lobe of the accretor. Large stellar mass loss rates such as the ones from the Wolf-Rayet star in M101 ULX-1 or the late B9 Supergiant in NGC 7793 P13 are enough to lead to a highly beamed wind and a significantly enhanced mass transfer rate.

In Supergiant X-ray binaries (SgXB), a compact object captures a fraction of the wind of an O/B supergiant on a close orbit. Proxies exist to evaluate the efficiency of mass and angular momentum accretion but they depend so dramatically on the wind speed that given the current uncertainties, they only set loose constrains. Furthermore, they often bypass the impact of orbital and shock effects on the flow structure. We study the wind dynamics and the angular momentum gained as the flow is accreted. We identify the conditions for the formation of a disc-like structure around the accretor and the observational consequences for SgXB. We use recent results on the wind launching mechanism to compute 3D streamlines, accounting for the gravitational and X-ray ionizing influence of the compact companion on the wind. Once the flow enters the Roche lobe of the accretor, we solve the hydrodynamics equations with cooling. A shocked region forms around the accretor as the flow is beamed. For wind speeds of the order of the orbital speed, the shock is highly asymmetric compared to the axisymmetric bow shock obtained for a purely planar homogeneous flow. With net radiative cooling, the flow always circularizes for wind speeds low enough. Although the donor star does not fill its Roche lobe, the wind can be significantly beamed and bent by the orbital effects. The net angular momentum of the accreted flow is then sufficient to form a persistent disc-like structure. This mechanism could explain the proposed limited outer extension of the accretion disc in Cygnus X-1 and suggests the presence of a disc at the outer rim of the neutron star magnetosphere in Vela X-1, with dramatic consequences on the spinning up of the accretor.

Hard X-ray (HXR) sources are frequently observed near the top of solar flare loops, and the emission is widely ascribed to bremsstrahlung. We here revisit an alternative scenario which stresses the importance of inverse Compton processes and the Kelvin- Helmholtz instability (KHI) proposed by Fang et al. (2016). This scenario adds a novel ingredient to the standard flare model, where evaporation flows from flare-impacted chromospheric foot-points interact with each other near the loop top and produce turbulence via KHI. The turbulence can act as a trapping region and as an efficient accelerator to provide energetic electrons, which scatter soft X-ray (SXR) photons to HXR photons via the inverse Compton mechanism. This paper focuses on the trigger of the KHI and the resulting turbulence in this new scenario. We perform a parameter survey to investigate the necessary ingredients to obtain KHI through interaction of chromospheric evaporation flows. When turbulence is produced in the loop apex, an index of -5/3 can be found in the spectra of velocity and magnetic field fluctuations. The KHI development and the generation of turbulence are controlled by the amount of energy deposited in the chromospheric foot-points and the time scale of its energy deposition, but typical values for M class flares show the KHI development routinely. Asymmetry of energy deposition determines the location where the turbulence is produced, and the synthesized SXR light curve shows a clear periodic signal related to the sloshing motion of the vortex pattern created by the KHI.

A solution to the problems of collimation and latency in relating the morphology of W 50 to SS 433

Solar prominences are subject to all kinds of perturbations during their lifetime, and frequently demonstrate oscillations. The study of prominence oscillations provides an alternative way to investigate their internal magnetic and thermal structures as the oscillation characteristics depend on their interplay with the solar corona. Prominence oscillations can be classified into longitudinal and transverse types. We perform three-dimensional ideal magnetohydrodynamic simulations of prominence oscillations along a magnetic flux rope, with the aim to compare the oscillation periods with those predicted by various simplified models and to examine the restoring force. We find that the longitudinal oscillation has a period of about 49 minutes, which is in accordance with the pendulum model where the field-ligned component of gravity serves as the restoring force. In contrast, the horizontal transverse oscillation has a period of about 10 minutes and the vertical transverse oscillation has a period of about 14 minutes, and both of them can be nicely fitted with a two-dimensional slab model. We also find that the magnetic tension force dominates most of the time in transverse oscillations, except for the first minute when magnetic pressure overwhelms.

Stratosphere circulation is important to interpret abundances of photo-chemically produced compounds like ozone that we aim to observe to assess habitability of exoplanets. We thus investigate a tidally locked ExoEarth scenario for TRAPPIST-1b, TRAPPIST-1d, Proxima Centauri~b and GJ 667 C~f with a simplified 3D atmosphere model and for different stratospheric wind breaking assumptions. These planets are representatives for different circulation regimes for orbital periods: $P_{orb}=1-100$~days. The circulation of exoplanets with $P_{orb} \leq $ 25~days can be dominated by the standing tropical Rossby wave in the troposphere and also in the stratosphere: It leads to a strong equatorial eastward wind jet and to 'Anti-Brewer-Dobson'-circulation that confines air masses to the stratospheric equatorial region. Thus, the distribution of photo-chemically produced species and aerosols may be limited to an 'equatorial transport belt'. In contrast, planets with $P_{orb}>25$~days, like GJ~667~C~f, exhibit efficient thermally driven circulation in the stratosphere that allows for a day side-wide distribution of air masses. The influence of the standing tropical Rossby waves on tidally locked ExoEarths with $P_{orb} \leq 25$~days can, however, be circumvented with deep stratospheric wind breaking alone - allowing for equator-to-pole transport like on Earth. For planets with $3 \leq P_{orb} \leq 6$~days, the extratropical Rossby wave acts as an additional safe-guard against the tropical Rossby wave in case of shallow wind breaking. Therefore, TRAPPIST-1d is less prone to have an equatorial transport belt in the stratosphere than Proxima~Centauri~b. Even our Earth model shows an equatorial wind jet, if stratosphere wind breaking is inefficient.

Supergiant X-ray Binaries (SgXB) host a compact object, often a neutron star (NS), orbiting an evolved O/B star. Mass transfer proceeds through the intense line-driven wind of the stellar donor, a fraction of which is captured by the gravitational field of the NS. The subsequent accretion process onto the NS is responsible for the abundant X-ray emission from SgXB. They also display peak-to-peak variability of the X-ray flux by a factor of a few 10 to 100, along with changes in the hardness ratios possibly due to varying absorption along the line-of-sight. We use recent radiation-hydrodynamic simulations of inhomogeneities (aka clumps) in the non-stationary wind of massive hot stars to evaluate their impact on the time-variable accretion process. For this, we run 3D hydrodynamic simulations of the wind in the vicinity of the accretor to investigate the formation of the bow shock and follow the inhomogeneous flow over several spatial orders of magnitude, down to the NS magnetosphere. In particular, we show that the impact of the wind clumps on the time-variability of the intrinsic mass accretion rate is severely tempered by the crossing of the shock, compared to the purely ballistic Bondi-Hoyle-Lyttleton estimation. We also account for the variable absorption due to clumps passing by the line-of-sight and estimate the final effective variability of the column density and mass accretion rate for different orbital separations. Finally, we compare our results to the most recent analysis of the X-ray flux and the hardness ratio in Vela X-1.

We report on the development of MPI-AMRVAC version 2.0, which is an open-source framework for parallel, grid-adaptive simulations of hydrodynamic and magnetohydrodynamic (MHD) astrophysical applications. The framework now supports radial grid stretching in combination with adaptive mesh refinement (AMR). The advantages of this combined approach are demonstrated with one-dimensional, two-dimensional and three-dimensional examples of spherically symmetric Bondi accretion, steady planar Bondi-Hoyle-Lyttleton flows, and wind accretion in Supergiant X-ray binaries. Another improvement is support for the generic splitting of any background magnetic field. We present several tests relevant for solar physics applications to demonstrate the advantages of field splitting on accuracy and robustness in extremely low plasma $\beta$ environments: a static magnetic flux rope, a magnetic null-point, and magnetic reconnection in a current sheet with either uniform or anomalous resistivity. Our implementation for treating anisotropic thermal conduction in multi-dimensional MHD applications is also described, which generalizes the original slope limited symmetric scheme from 2D to 3D. We perform ring diffusion tests that demonstrate its accuracy and robustness, and show that it prevents the unphysical thermal flux present in traditional schemes. The improved parallel scaling of the code is demonstrated with 3D AMR simulations of solar coronal rain, which show satisfactory strong scaling up to 2000 cores. Other framework improvements are also reported: the modernization and reorganization into a library, the handling of automatic regression tests, the use of inline/online Doxygen documentation, and a new future-proof data format for input/output

An exospheric kinetic solar wind model is interfaced with an observation-driven single fluid magnetohydrodynamic (MHD) model. Initially, a photospheric magnetogram serves as observational input in the fluid approach to extrapolate the heliospheric magnetic field. Then semi-empirical coronal models are used for estimating the plasma characteristics up to a heliocentric distance of 0.1AU. From there on a full MHD model which computes the three-dimensional time-dependent evolution of the solar wind macroscopic variables up to the orbit of the Earth is used. After interfacing the density and velocity at the inner MHD boundary, we compare with the results of a kinetic exospheric solar wind model based on the assumption of Maxwell and Kappa velocity distribution functions for protons and electrons respectively, as well as with \textit{in situ} observations at 1AU. This provides insight on more physically detailed processes, such as coronal heating and solar wind acceleration, that naturally arise by inclusion of suprathermal electrons in the model. We are interested in the profile of the solar wind speed and density at 1AU, in characterizing the slow and fast source regions of the wind and in comparing MHD with exospheric models in similar conditions. We calculate the energetics of both models from low to high heliocentric distances.

Several observations of astrophysical jets show evidence of a structure in the direction perpendicular to the jet axis, leading to the development of "spine & sheath" models of jets. Most studies focus on a two-component jet consisting of a highly relativistic inner jet and a slower - but still relativistic - outer jet surrounded by an unmagnetized environment. These jets are believed to be susceptible to a relativistic Rayleigh-Taylor-type instability, depending on the effective inertia ratio of the two components. We extend previous studies by taking into account the presence of a non-zero toroidal magnetic field. Different values of magnetization are examined, to detect possible differences in the evolution and stability of the jet. We find that the toroidal field, above a certain level of magnetization $\sigma$, roughly equal to 0.01, can stabilize the jet against the previously mentioned instabilities and that there is a clear trend in the behaviour of the average Lorentz factor and the effective radius of the jet when we continuously increase the magnetization. The simulations are performed using the relativistic MHD module from the open source, parallel, grid adaptive, MPI-AMRVAC code.

We intend to investigate the underlying physics for the coronal rain phenomenon in a representative bipolar magnetic field, including the formation and the dynamics of coronal rain blobs. With the MPI-AMRVAC code, we performed three dimensional radiative magnetohydrodynamic (MHD) simulation with strong heating localized on footpoints of magnetic loops after a relaxation to quiet solar atmosphere. Progressive cooling and in-situ condensation starts at the loop top due to radiative thermal instability. The first large-scale condensation on the loop top suffers Rayleigh-Taylor instability and becomes fragmented into smaller blobs. The blobs fall vertically dragging magnetic loops until they reach low beta regions and start to fall along the loops from loop top to loop footpoints. A statistic study of the coronal rain blobs finds that small blobs with masses of less than 10^10 g dominate the population. When blobs fall to lower regions along the magnetic loops, they are stretched and develop a non-uniform velocity pattern with an anti-parallel shearing pattern seen to develop along the central axis of the blobs. Synthetic images of simulated coronal rain with Solar Dynamics Observatory Atmospheric Imaging Assembly well resemble real observations presenting dark falling clumps in hot channels and bright rain blobs in a cool channel. We also find density inhomogeneities during a coronal rain "shower", which reflects the observed multi-stranded nature of coronal rain.

Magnetic reconnection and non-thermal particle distributions associated with current-driven instabilities are investigated by means of resistive magnetohydrodynamics (MHD) simulations combined with relativistic test particle methods. We propose a system with two parallel, repelling current channels in an initially force-free equilibrium, as a simplified representation of flux ropes in a stellar magnetosphere. The current channels undergo a rotation and separation on Alfv\'enic timescales, forming secondary islands and (up to tearing unstable) current sheets in which non-thermal energy distributions are expected to develop. Using the recently developed particle module of our open-source grid-adaptive MPI-AMRVAC software, we simulate MHD evolution combined with test particle treatments in MHD snapshots. We explore under which plasma-$\beta$ conditions the fastest reconnection occurs in two-and-a-half dimensional (2.5D) scenarios and in these settings test particles are evolved. We quantify energy distributions, acceleration mechanisms, relativistic corrections to the particle equations of motion and effects of resistivity in magnetically dominated proton-electron plasmas. Due to large resistive electric fields and indefinite acceleration of particles in the infinitely long current channels, hard energy spectra are found in 2.5D configurations. Solutions to these numerical artifacts are proposed for both 2.5D setups and future 3D work. We discuss the magnetohydrodynamics of an additional kink instability in 3D setups and the expected effects on energy distributions. The obtained results hold as a proof-of-principle for test particle approaches in MHD simulations, relevant to explore less idealised scenarios like solar flares and more exotic astrophysical phenomena, like black hole flares, magnetar magnetospheres and pulsar wind nebulae.

We investigate how night side cooling and surface friction impact surface temperatures and large scale circulation for tidally locked Earth-like planets. For each scenario, we vary the orbital period between $P_{rot}=1-100$~days and capture changes in climate states. We find drastic changes in climate states for different surface friction scenarios. For very efficient surface friction ($t_{s,fric}=$ 0.1 days), the simulations for short rotation periods ($P_{rot} \leq$ 10 days) show predominantly standing extra tropical Rossby waves. These waves lead to climate states with two high latitude westerly jets and unperturbed meridional direct circulation. In most other scenarios, simulations with short rotation periods exhibit instead dominance by standing tropical Rossby waves. Such climate states have a single equatorial westerly jet, which disrupts direct circulation. Experiments with weak surface friction ($t_{s,fric}=~10 -100$ days) show decoupling between surface temperatures and circulation, which leads to strong cooling of the night side. The experiment with $t_{s,fric}= 100$ days assumes climate states with easterly flow (retrograde rotation) for medium and slow planetary rotations $P_{rot}= 12 - 100$~days. We show that an increase of night side cooling efficiency by one order of magnitude compared to the nominal model leads to a cooling of the night side surface temperatures by 80-100~K. The day side surface temperatures only drop by 25~K at the same time. The increase in thermal forcing suppresses the formation of extra tropical Rossby waves on small planets ($R_P=1 R_{Earth}$) in the short rotation period regime ($P_{rot} \leq$ 10 days).

The Wolf-Rayet 98a (WR 98a) system is a prime target for interferometric surveys, since its identification as a "rotating pinwheel nebulae", where infrared images display a spiral dust lane revolving with a 1.4 year periodicity. WR 98a hosts a WC9+OB star, and the presence of dust is puzzling given the extreme luminosities of Wolf-Rayet stars. We present 3D hydrodynamic models for WR 98a, where dust creation and redistribution are self-consistently incorporated. Our grid-adaptive simulations resolve details in the wind collision region at scales below one percent of the orbital separation (~4 AU), while simulating up to 1300 AU. We cover several orbital periods under conditions where the gas component alone behaves adiabatic, or is subject to effective radiative cooling. In the adiabatic case, mixing between stellar winds is effective in a well-defined spiral pattern, where optimal conditions for dust creation are met. When radiative cooling is incorporated, the interaction gets dominated by thermal instabilities along the wind collision region, and dust concentrates in clumps and filaments in a volume-filling fashion, so WR 98a must obey close to adiabatic evolutions to demonstrate the rotating pinwheel structure. We mimic Keck, ALMA or future E-ELT observations and confront photometric long-term monitoring. We predict an asymmetry in the dust distribution between leading and trailing edge of the spiral, show that ALMA and E-ELT would be able to detect fine-structure in the spiral indicative of Kelvin-Helmholtz development, and confirm the variation in photometry due to the orientation. Historic Keck images are reproduced, but their resolution is insufficient to detect the details we predict.

Solar prominences are long-lived cool and dense plasma curtains in the hot and rarefied outer solar atmosphere or corona. The physical mechanism responsible for their formation and especially for their internal plasma circulation has been uncertain for decades. The observed ubiquitous down flows in quiescent prominences are difficult to interpret as plasma with high conductivity seems to move across horizontal magnetic field lines. Here we present three-dimensional numerical simulations of prominence formation and evolution in an elongated magnetic flux rope as a result of in-situ plasma condensations fueled by continuous plasma evaporation from the solar chromosphere. The prominence is born and maintained in a fragmented, highly dynamic state with continuous reappearance of multiple blobs and thread structures that move mainly downward dragging along mass-loaded field lines. The prominence plasma circulation is characterized by the dynamic balance between the drainage of prominence plasma back to the chromosphere and the formation of prominence plasma via continuous condensation. Plasma evaporates from the chromosphere, condenses into the prominence in the corona, and drains back to the chromosphere, establishing a stable chromosphere-corona plasma cycle. Synthetic images of the modeled prominence with the Solar Dynamics Observatory Atmospheric Imaging Assembly closely resemble actual observations, with many dynamical threads underlying an elliptical coronal cavity.

In this paper, we produce synthetic radio views on simulated flux ropes in the solar corona, where finite-beta magnetohydrodynamic (MHD) simulations serve to mimic the flux rope formation stages, as well as their stable endstates. These endstates represent twisted flux ropes where balancing Lorentz forces, gravity and pressure gradients determine the full thermodynamic variation throughout the flux rope. The obtained models are needed to quantify radiative transfer in radio bands, and allow us to contrast weak to strong magnetic field conditions. Field strengths of up to 100 G in the flux rope yield the radio views dominated by optically thin free-free emission. The forming flux rope shows clear morphological changes in its emission structure as it deforms from an arcade to a flux rope, both on disk and at the limb. For an active region filament channel with a field strength of up to 680 G in the flux rope, gyroresonance emission (from the third-fourth gyrolayers) can be detected and even dominates over free-free emission at the frequencies of up to 7 GHz. Finally, we also show synthetic views on a simulated filament embedded within a (weak-field) flux rope, resulting from an energetically consistent MHD simulation. For this filament, synthetic views at the limb show clear similarities with actual observations, and the transition from optically thick (below 10 GHz) to optically thin emission can be reproduced. On the disk, its dimension and temperature conditions are as yet not realistic enough to yield the observed radio brightness depressions.

Quasi-periodic propagating intensity disturbances have been observed in large coronal loops in EUV images over a decade, and are widely accepted to be slow magnetosonic waves. However, spectroscopic observations from Hinode/EIS revealed their association with persistent coronal upflows, making this interpretation debatable. We perform a 2.5D magnetohydrodynamic simulation to imitate the chromospheric evaporation and the following reflected patterns in a flare loop. Our model encompasses the corona, transition region, and chromosphere. We demonstrate that the quasi periodic propagating intensity variations captured by the synthesized \textit{Solar Dynamics Observatory}/Atmospheric Imaging Assembly (AIA) 131, 94~\AA~emission images match the previous observations well. With particle tracers in the simulation, we confirm that these quasi periodic propagating intensity variations consist of reflected slow mode waves and mass flows with an average speed of 310 km/s in an 80 Mm length loop with an average temperature of 9 MK. With the synthesized Doppler shift velocity and intensity maps of the \textit{Solar and Heliospheric Observatory}/Solar Ultraviolet Measurement of Emitted Radiation (SUMER) Fe XIX line emission, we confirm that these reflected slow mode waves are propagating waves.

We investigate 3D atmosphere dynamics for tidally locked terrestrial planets with an Earth-like atmosphere and irradiation for different rotation periods ($P_{rot}=1-100$ days) and planet sizes ($R_P=1-2 R_{Earth}$) with unprecedented fine detail. We could precisely identify three climate state transition regions that are associated with phase transitions in standing tropical and extra tropical Rossby waves. We confirm that the climate on fast rotating planets may assume multiple states ($P_{rot}\leq 12$ days for $R_P=2 R_{Earth}$). Our study is, however, the first to identify the type of planetary wave associated with different climate states: The first state is dominated by standing tropical Rossby waves with fast equatorial superrotation. The second state is dominated by standing extra tropical Rossby waves with high latitude westerly jets with slower wind speeds. For very fast rotations ($P_{rot}\leq 5$~days for $R_P=2 R_{Earth}$), we find another climate state transition, where the standing tropical and extra tropical Rossby wave can both fit on the planet. Thus, a third state with a mixture of the two planetary waves becomes possible that exhibits three jets. Different climate states may be observable, because the upper atmosphere's hot spot is eastward shifted with respect to the substellar point in the first state, westward shifted in the second state and the third state shows a longitudinal 'smearing' of the spot across the substellar point. We show, furthermore, that the largest fast rotating planet in our study exhibits atmosphere features known from hot Jupiters like fast equatorial superrotation and a temperature chevron in the upper atmosphere.

We extend our earlier multidimensional, magnetohydrodynamic simulations of coronal rain occurring in magnetic arcades with higher resolution, grid-adaptive computations covering a much longer ($>6$ hour) timespan. We quantify how in-situ forming blob-like condensations grow along and across field lines and show that rain showers can occur in limit cycles, here demonstrated for the first time in 2.5D setups. We discuss dynamical, multi-dimensional aspects of the rebound shocks generated by the siphon inflows and quantify the thermodynamics of a prominence-corona-transition-region like structure surrounding the blobs. We point out the correlation between condensation rates and the cross-sectional size of loop systems where catastrophic cooling takes place. We also study the variations of the typical number density, kinetic energy and temperature while blobs descend, impact and sink into the transition region. In addition, we explain the mechanisms leading to concurrent upflows while the blobs descend. As a result, there are plenty of shear flows generated with relative velocity difference around 80 km s$^{-1}$ in our simulations. These shear flows are siphon flows set up by multiple blob dynamics and they in turn affect the deformation of the falling blobs. In particular, we show how shear flows can break apart blobs into smaller fragments, within minutes.

Observations revealed rich dynamics within prominences, the cool 10,000 K, macroscopic (sizes of order 100 Mm) "clouds" in the million degree solar corona. Even quiescent prominences are continuously perturbed by hot, rising bubbles. Since prominence matter is hundredfold denser than coronal plasma, this bubbling is related to Rayleigh-Taylor instabilities. Here we report on true macroscopic simulations well into this bubbling phase, adopting a magnetohydrodynamic description from chromospheric layers up to 30 Mm height. Our virtual prominences rapidly establish fully non-linear (magneto)convective motions where hot bubbles interplay with falling pillars, with dynamical details including upwelling pillars forming within bubbles. Our simulations show impacting Rayleigh-Taylor fingers reflecting on transition region plasma, ensuring that cool, dense chromospheric material gets mixed with prominence matter up to very large heights. This offers an explanation for the return mass cycle mystery for prominence material. Synthetic views at extreme ultraviolet wavelengths show remarkable agreement with observations, with clear indications of shear-flow induced fragmentations.

We present numerical simulations in 3D settings where coronal rain phenomena take place in a magnetic configuration of a quadrupolar arcade system. Our simulation is a magnetohydrodynamic simulation including anisotropic thermal conduction, optically thin radiative losses, and parametrised heating as main thermodynamical features to construct a realistic arcade configuration from chromospheric to coronal heights. The plasma evaporation from chromospheric and transition region heights eventually causes localised runaway condensation events and we witness the formation of plasma blobs due to thermal instability, that evolve dynamically in the heated arcade part and move gradually downwards due to interchange type dynamics. Unlike earlier 2.5D simulations, in this case there is no large scale prominence formation observed, but a continuous coronal rain develops which shows clear indications of Rayleigh-Taylor or interchange instability, that causes the denser plasma located above the transition region to fall down, as the system moves towards a more stable state. Linear stability analysis is used in the non-linear regime for gaining insight and giving a prediction of the system's evolution. After the plasma blobs descend through interchange, they follow the magnetic field topology more closely in the lower coronal regions, where they are guided by the magnetic dips.

In light of the recent detection of direct evidence for the formation of Kelvin-Helmholtz instabilities in the Orion nebula, we expand upon previous modelling efforts by numerically simulating the shear-flow driven gas and dust dynamics in locations where the H$_{II}$ region and the molecular cloud interact. We aim to directly confront the simulation results with the infrared observations. Methods: To numerically model the onset and full nonlinear development of the Kelvin-Helmholtz instability we take the setup proposed to interpret the observations, and adjust it to a full 3D hydrodynamical simulation that includes the dynamics of gas as well as dust. A dust grain distribution with sizes between 5-250 nm is used, exploiting the gas+dust module of the MPI-AMRVAC code, in which the dust species are represented by several pressureless dust fluids. The evolution of the model is followed well into the nonlinear phase. The output of these simulations is then used as input for the SKIRT dust radiative transfer code to obtain infrared images at several stages of the evolution, which can be compared to the observations. Results: We confirm that a 3D Kelvin-Helmholtz instability is able to develop in the proposed setup, and that the formation of the instability is not inhibited by the addition of dust. Kelvin-Helmholtz billows form at the end of the linear phase, and synthetic observations of the billows show striking similarities to the infrared observations. It is pointed out that the high density dust regions preferentially collect on the flanks of the billows. To get agreement with the observed Kelvin-Helmholtz ripples, the assumed geometry between the background radiation, the billows and the observer is seen to be of critical importance.

Magnetohydrodynamic waves interact with structured plasmas and reveal the internal magnetic and thermal structures therein, thereby having seismological applications in the solar atmosphere. We investigate the evolution of fast magnetoacoustic pulses in randomly structured plasmas, in the context of large-scale propagating waves in the solar atmosphere. We perform one dimensional numerical simulations of fast wave pulses propagating perpendicular to a constant magnetic field in a low-$\beta$ plasma with a random density profile across the field. Both linear and nonlinear regimes are considered. We study how the evolution of the pulse amplitude and width depends on their initial values and the parameters of the random structuring. A randomly structured plasma acts as a dispersive medium for a fast magnetoacoustic pulse, causing amplitude attenuation and broadening of the pulse width. After the passage of the main pulse, secondary propagating and standing fast waves appear in the plasma. Width evolution of both linear and nonlinear pulses can be well approximated by linear functions; however, narrow pulses may have zero or negative broadening. This arises because a narrow pulse is prone to splitting, while a broad pulse usually deviates less from their initial Gaussian shape and form ripple structures on top of the main pulse. A linear pulse decays at almost a constant rate, while a nonlinear pulse decays exponentially. A pulse interacts most efficiently with a random medium which has a correlation length of about half of its initial pulse width. The development of a detailed model of a fast MHD pulse propagating in highly structured medium substantiates the interpretation of EIT waves as fast magnetoacoustic waves. Evolution of a fast pulse provides us with a novel method to diagnose the sub-resolution filamentation of the solar atmosphere.

We present a numerical study in resistive magnetohydrodynamics where the initial equilibrium configuration contains adjacent, oppositely directed, parallel current channels. Since oppositely directed current channels repel, the equilibrium is liable to an ideal magnetohydrodynamic tilt instability. This tilt evolution, previously studied in planar settings, involves two magnetic islands or fluxropes, which on Alfvenic timescales undergo a combined rotation and separation. This in turn leads to the creation of (near) singular current layers, posing severe challenges to numerical approaches. Using our open-source grid-adaptive MPI-AMRVAC software, we revisit the planar evolution case in compressible MHD, as well as its extension to 2.5D and full 3D scenarios. As long as the third dimension remains ignorable, pure tilt evolutions result which are hardly affected by out of plane magnetic field components. In all 2.5D runs, our simulations do show secondary tearing type disruptions throughout the near singular current sheets in the far nonlinear saturation regime. In full 3D, both current channels can be liable to additional ideal kink deformations. We discuss the effects of having both tilt and kink instabilities acting simultaneously in the violent, reconnection dominated evolution. In 3D, both the tilt and the kink instabilities can be stabilized by tension forces. As a concrete space plasma application, we argue that interacting tilt-kink instabilities in repelling current channels provide a novel route to initiate solar coronal mass ejections, distinctly different from currently favored pure kink or torus instability routes.

Radiative equilibrium temperatures are calculated for the troposphere of a tidally locked Super-Earth based on a simple greenhouse model, using Solar System data as a guideline. These temperatures provide in combination with a Newtonian relaxation scheme thermal forcing for a 3D atmosphere model using the dynamical core of the Massachusetts Institute of Technology global circulation model (MITgcm). Our model is of the same conceptional simplicity than the model of Held & Suarez (1994) and is thus computationally fast. Furthermore, because of the coherent, general derivation of radiative equilibrium temperatures, our model is easily adaptable for different planets and atmospheric scenarios. As a case study relevant for Super-Earths, we investigate a Gl581g-like planet with Earth-like atmosphere and irradiation and present results for two representative rotation periods of $P_{rot}$ = 10 days and $P_{rot}$ = 36.5 days. Our results provide proof of concept and highlight interesting dynamical features for the rotating regime 3 < $P_{rot}$ < 100 days, which was shown by Edson et al. (2011) to be an intermediate regime between equatorial superrotation and divergence. We confirm that the $P_{rot}$ = 10 days case is more dominated by equatorial superrotation dynamics than the $P_{rot}$ = 36.5 days case, which shows diminishing influence of standing Rossby- Kelvin waves and increasing influence of divergence at the top of the atmosphere. We argue that this dynamical regime change relates to the increase in Rossby deformation radius, in agreement with previous studies. However, we also pay attention to other features that are not or only in partial agreement with other studies, like, e.g., the number of circulation cells and their strength, the role and extent of thermal inversion layers, and the details of heat transport

Prominences in the solar corona are hundredfold cooler and denser than their surroundings, with a total mass of 1.e13 up to 1.e15 g. Here we report on the first comprehensive simulations of three-dimensional, thermally and gravitationally stratified magnetic flux ropes, where in situ condensation to a prominence happens due to radiative losses. After a gradual thermodynamic adjustment, we witness a phase where runaway cooling happens while counter-streaming shearing flows drain off mass along helical field lines. After this drainage, a prominence-like condensation resides in concave upward field regions, and this prominence retains its overall characteristics for more than two hours. While condensing, the prominence establishes a prominence-corona transition region, where magnetic field-aligned thermal conduction is operative during the runaway cooling. The prominence structure represents a force-balanced state in a helical flux rope. The simulated condensation demonstrates a right-bearing barb, as a remnant of the drainage. Synthetic images at extreme ultraviolet wavelengths follow the onset of the condensation, and confirm the appearance of horns and a three-part structure for the stable prominence state, as often seen in erupting prominences. This naturally explains recent Solar Dynamics Observatory views with the Atmospheric Imaging Assembly on prominences in coronal cavities demonstrating horns.

In this paper we present an update on the open source MPI-AMRVAC simulation toolkit where we focus on solar- and non-relativistic astrophysical magneto-fluid dynamics. We highlight recent developments in terms of physics modules such as hydrodynamics with dust coupling and the conservative implementation of Hall magnetohydrodynamics. A simple conservative high-order finite difference scheme that works in combination with all available physics modules is introduced and demonstrated at the example of monotonicity preserving fifth order reconstruction. Strong stability preserving high order Runge-Kutta time steppers are used to obtain stable evolutions in multidimensional applications realizing up to fourth order accuracy in space and time. With the new distinction between active and passive grid cells, MPI-AMRVAC is ideally suited to simulate evolutions where parts of the solution are controlled analytically, or have a tendency to progress into or out of a stationary state. Typical test problems and representative applications are discussed, with an outlook to follow-up research. Finally, we discuss the parallel scaling of the code and demonstrate excellent weak scaling up to 30 000 processors allowing to exploit modern petascale infrastructure.

In this paper we discuss the development of Rayleigh-Taylor filaments in axisymmetric simulations of Pulsar wind nebulae (PWN). High-resolution adaptive mesh refinement magnetohydrodynamic (MHD) simulations are used to resolve the non-linear evolution of the instability. The typical separation of filaments is mediated by the turbulent flow in the nebula and hierarchical growth of the filaments. The strong magnetic dissipation and field-randomization found in recent global three-dimensional simulations of PWN suggests that magnetic tension is not strong enough to suppress the growth of RT filaments, in agreement with the observations of prominent filaments in the Crab nebula. The long-term axisymmetric results presented here confirm this finding.

We present numerical simulations in 2.5D settings where large scale prominences form in situ out of coronal condensation in magnetic dips, in close agreement with early as well as recent reporting of `funnel prominences'. Our simulation uses full thermodynamic MHD with anisotropic thermal conduction, optically thin radiative losses, and parametrized heating as main ingredients to establish a realistic arcade configuration from chromosphere to corona. The chromospheric evaporation from especially transition region heights ultimately causes thermal instability and we witness the growth of a prominence suspended well above the transition region, continuously gaining mass and cross-sectional area. Several hours later, the condensation has grown into a structure connecting the prominence-corona transition region with the underlying transition region, and a continuous downward motion from the accumulated mass represents a drainage that matches observational findings. A more dynamic phase is found as well, with coronal rain, induced wave trains, and even a reconnection event when the core prominence plasma weighs down the fieldlines until a fluxrope gets formed. The upper part of the prominence is then trapped in a fluxrope structure, and we argue for its violent kink-unstable eruption as soon as the (ignored) length dimension would allow for ideal kink deformations.

Various radio galaxies show signs of having gone through episodic jet outbursts in the past. An example is the class of double-double radio galaxies (DDRGs). However, to follow the evolution of an individual source in real-time is impossible due to the large time scales involved. Numerical studies provide a powerful tool to investigate the temporal behavior of episodic jet outbursts in a (magneto-)hydrodynamical setting. We simulate the injection of two jets from active galactic nuclei (AGN), separated by a short interruption time. Three different jet models are compared. We find that an AGN jet outburst cycle can be divided into four phases. The most prominent phase occurs when the restarted jet is propagating completely inside the hot and inflated cocoon left behind by the initial jet. In that case, the jet-head advance speed of the restarted jet is significantly higher than the initial jet-head. While the head of the initial jet interacts strongly with the ambient medium, the restarted jet propagates almost unimpeded. As a result, the restarted jet maintains a strong radial integrity. Just a very small fraction of the amount of shocked jet material flows back through the cocoon compared to that of the initial jet and much weaker shocks are found at the head of the restarted jet. We find that the features of the restarted jet in this phase closely resemble the observed properties of a typical DDRG.

Dust is present in a large variety of astrophysical fluids, from tori around supermassive black holes to molecular clouds, protoplanetary discs, and cometary outflows. In many such fluids, shearing flows are present, leading to the formation of Kelvin-Helmholtz instabilities (KHI) and changing the properties and structures of the fluid through processes such as mixing and clumping of dust. We investigate how dust changes the growth rates of the KHI in 2D and 3D and how the it redistributes and clumps dust. We investigate if similarities can be found between the structures in 3D KHI and those seen in observations of molecular clouds. We do this by performing numerical hydrodynamical dust+gas simulations with in addition to gas a number of dust fluids. Each dust fluid represents a portion of the particle size-distribution. We study how dust-to-gas mass density ratios between 0.01 and 1 alter the growth rate in the linear phase of the KHI. We do this for a wide range of perturbation wavelengths, and compare these values to the analytical gas-only growth rates. As the formation of high-density dust structures is of interest in many astrophysical environments, we scale our simulations with physical quantities similar to values in molecular clouds. Large differences in dynamics are seen for different grain sizes. We demonstrate that high dust-to-gas ratios significantly reduce the growth rate of the KHI, especially for short wavelengths. We compare the dynamics in 2D and 3D simulations, where the latter demonstrates additional full 3D instabilities during the non-linear phase, leading to increased dust densities. We compare the structures formed by the KHI in 3D simulations with those in molecular clouds and see how the column density distribution of the simulation shares similarities with log-normal distributions with power-law tails sometimes seen in observations of molecular clouds.

Context. Modern high resolution observations allow to view closer into the objects powering relativistic jets. This is especially the case for SS433, an X-ray binary from which a precessing jet is observed down to the sub-parsec scale. Aims. We want to study full 3D dynamics of relativistic jets associated with AGN or XRB. We study the precessing motion of a jet as a model for the jet associated with the XRB SS433. Our study of the jet dynamics in this system focuses on the sub-parsec scales. We investigate the impact of jet precession and the variation of the Lorentz factor of the injected matter on the general 3D jet dynamics and its energy transfer to the surrounding medium. We realize synthetic radio mapping of the data, to compare our results with observations. Methods. For our study we use the code MPI-AMRVAC with SRHD model of a baryonic jet. We use a AMR scheme and an inner time-dependent boundary prescription to inject the jets. Parameters extracted from observations were used. 3D jet realizations that match the SS433 jet are intercompared. We track the energy content, as deposited in different regions of the domain affected by the jet. Our code also follows a population of particles injected with the jet. This evolving energy spectrum of accelerated electrons, allows to obtain the radio emission from our simulation. Results. we obtain meaningful observations. We find increased energy transfer for a precessing jet compared to a standing jets. We obtain synthetic radio maps for all jets. Conclusions. The synthetic radio map matches best for a model using the canonical kinematic model. Overdense precessing jets experience significant deceleration in their propagation through the ISM, and while the overall jet is of helical shape. This argument show that the kinematic model for SS433 assuming ballistic propagation has to be corrected for this decelretaion.

The magnetic configuration hosting prominences and their surrounding coro- nal structure is a key research topic in solar physics. Recent theoretical and observational studies strongly suggest that a helical magnetic flux rope is an es- sential ingredient to fulfill most of the theoretical and observational requirements for hosting prominences. To understand flux rope formation details and obtain magnetic configurations suitable for future prominence formation studies, we here report on three-dimensional isothermal magnetohydrodynamic simulations including finite gas pressure and gravity. Starting from a magnetohydrostatic corona with a linear force-free bipolar magnetic field, we follow its evolution when introducing vortex flows around the main polarities and converging flows towards the polarity inversion line near the bottom of the corona. The con- verging flows bring feet of different loops together at the polarity inversion line and magnetic reconnection and flux cancellation happens. Inflow and outflow signatures of the magnetic reconnection process are identified, and the thereby newly formed helical loops wind around pre-existing ones so that a complete flux rope grows and ascends. When a macroscopic flux rope is formed, we switch off the driving flows and find that the system relaxes to a stable state containing a helical magnetic flux rope embedded in an overlying arcade structure. A major part of the formed flux rope is threaded by dipped field lines which can stably support prominence matter, while the total mass of the flux rope is in the order of 4-5.e14 g.

The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime.

In this paper we give a detailed account of the first 3D relativistic magnetohydrodynamic (MHD) simulations of Pulsar Wind Nebulae (PWN), with parameters most suitable for the Crab Nebula. In order to clarify the new features specific to 3D models, reference 2D simulations have been carried out as well. Compared to the previous 2D simulations, we considered pulsar winds with much stronger magnetisation, up to \sigma=3, and accounted more accurately for the anticipated magnetic dissipation in the striped zone of these winds. While the 3D models preserve the separation of the post termination shock flow into the equatorial and polar components, their relative strength and significance differ. Whereas the highly magnetised 2D models produce highly coherent and well collimated polar jets capable of efficient "drilling" through the supernova shell, in the corresponding 3D models the jets are disrupted by the kink mode current driven instability and "dissolve" into the main body of PWN after propagation of several shock radii. Our results show that contrary to the expectations based on 1D analytical and semi-analytical models, our numerical solutions with highly magnetized pulsar winds still produce termination shocks with radii comparable to those deduced from the observations. We present polarized synchrotron maps constructed from the 3D simulations, showing that the wealth of morphological features observed with Hubble and Chandra can well be reproduced with 3D-MHD models.

We present the first multidimensional, magnetohydrodynamic simulations which capture the initial formation and the long-term sustainment of the enigmatic coronal rain phenomenon. We demonstrate how thermal instability can induce a spectacular display of in-situ forming blob-like condensations which then start their intimate ballet on top of initially linear force-free arcades. Our magnetic arcades host chromospheric, transition region, and coronal plasma. Following coronal rain dynamics for over 80 minutes physical time, we collect enough statistics to quantify blob widths, lengths, velocity distributions, and other characteristics which directly match with modern observational knowledge. Our virtual coronal rain displays the deformation of blobs into $V$-shaped like features, interactions of blobs due to mostly pressure-mediated levitations, and gives the first views on blobs which evaporate in situ, or get siphoned over the apex of the background arcade. Our simulations pave the way for systematic surveys of coronal rain showers in true multidimensional settings, to connect parametrized heating prescriptions with rain statistics, ultimately allowing to quantify the coronal heating input.

Current observations have shown that astrophysical jets reveal strong signs of radial structure. They suggest that the inner region of the jet, the jet spine, consists of a low-density, fast-moving gas, while the outer region of the jet consists of a more dense and slower moving gas, called the jet sheath. Moreover, if jets carry angular momentum, the resultant centrifugal forces lead to a radial stratification. Current observations are not able to fully resolve the radial structure, so little is known about its actual profile. We present three AGN jet models in $2.5D$ of which two have been given a radial structure. The first model is a homogeneous jet, the only model that doesn't carry angular momentum; the second model is a spine-sheath jet with an isothermal equation of state; and the third jet model is a (piecewise) isochoric spine-sheath jet, with constant but different densities for jet spine and jet sheath. In this paper, we look at the effects of radial stratification on jet integrity, mixing between the different jet components and global morphology of the jet-head and surrounding cocoon.

We present in this paper both a linear study and numerical relativistic MHD simulations of the non-resonant streaming instability occurring in the precursor of relativistic shocks. In the shock front restframe, we perform a linear analysis of this instability in a likely configuration for ultra-relativistic shock precursors. This considers magneto-acoustic waves having a wave vector perpendicular to the shock front and the large scale magnetic field. Our linear analysis is achieved without any assumption on the shock velocity and is thus valid for all velocity regimes. In order to check our calculation, we also perform relativistic MHD simulations describing the propagation of the aforementioned magneto-acoustic waves through the shock precursor. The numerical calculations confirm our linear analysis, which predicts that the growth rate of the instability is maximal for ultra-relativistic shocks and exhibits a wavenumber dependence $\propto k_x^{1/2}$. Our numerical simulations also depict the saturation regime of the instability where we show that the magnetic amplification is moderate but nevertheless significant ($\delta B/B\leq 10$). This latter fact may explain the presence of strong turbulence in the vicinity of relativistic magnetized shocks. Our numerical approach also introduces a convenient means to handle isothermal (ultra-)relativistic MHD conditions.

It is found that both microflare-sized impulsive heating at one leg of the loop and a suddenly imposed velocity perturbation can propel the prominence to oscillate along the magnetic dip. An extensive parameter survey results in a scaling law, showing that the period of the oscillation, which weakly depends on the length and height of the prominence, and the amplitude of the perturbations, scales with $\sqrt{R/g_\odot}$, where $R$ represents the curvature radius of the dip, and $g_\odot$ is the gravitational acceleration of the Sun. This is consistent with the linear theory of a pendulum, which implies that the field-aligned component of gravity is the main restoring force for the prominence longitudinal oscillations, as confirmed by the force analysis. However, the gas pressure gradient becomes non-negligible for short prominences. The oscillation damps with time in the presence of non-adiabatic processes. Compared to heat conduction, the radiative cooling is the dominant factor leading to the damping. A scaling law for the damping timescale is derived, i.e., $\tau\sim l^{1.63} D^{0.66}w^{-1.21}v_{0}^{-0.30}$, showing strong dependence on the prominence length $l$, the geometry of the magnetic dip (characterized by the depth $D$ and the width $w$), and the velocity perturbation amplitude $v_0$. The larger the amplitude, the faster the oscillation damps. It is also found that mass drainage significantly reduces the damping timescale when the perturbation is too strong.

We present first results of three dimensional relativistic magnetohydrodynamical simulations of Pulsar Wind Nebulae. They show that the kink instability and magnetic dissipation inside these nebulae may be the key processes allowing to reconcile their observations with the theory of pulsar winds. In particular, the size of the termination shock, obtained in the simulations, agrees very well with the observations even for Poynting-dominated pulsar winds. Due to magnetic dissipation the total pressure in the simulated nebulae is particle-dominated and more or less uniform. While in the main body of the simulated nebulae the magnetic field becomes rather randomized, close to the termination shock, it is dominated by the regular toroidal field freshly injected by the pulsar wind. This field is responsible for driving polar outflows and may explain the high polarization observed in pulsar wind nebulae.

Context. Relativistic jets emerging from AGN cores transfer energy from the core to their surrounding ISM/IGM. Because jets are observed to have finite opening angles, one needs to quantify the role of conical versus cylindrical jet propagation in this energy transfer. Aims. We use FR-II AGN jets parameter with finite opening angles. We study the effect of the variation of the opening angle on the dynamics and energy transfer of the jet. We also point out how the characteristics of this external medium, such as its density profile, play a role in the dynamics. Methods. This study exploits our parallel AMR code MPI-AMRVAC with its special relativistic hydrodynamic model, incorporating an equation of state with varying effective polytropic index. We studied mildly under-dense jets up to opening angles of 10 degrees, at Lorentz factors of about 10, inspired by observations. Instantaneous quantification of the various ISM volumes and their energy content allows one to quantify the role of mixing versus shock-heated cocoon regions over the time intervals. Results. We show that a wider opening angle jet results in a faster deceleration of the jet and leads to a wider cocoon dominated by Kelvin-Helmholtz and Rayleigh-Taylor instabilities. The energy transfer mainly occurs in the shocked ISM region by both the frontal bow shock and cocoon-traversing shock waves, in a roughly 3 to 1 ratio to the energy transfer of the mixing zone, for a 5 degree opening angle jet. A rarefaction wave induces a dynamically formed layered structure of the jet beam. Conclusions. Finite opening angle jets can efficiently transfer significant fractions (25 % up to 70 %) of their injected energy over a growing region of shocked ISM matter. The role of the ISM stratification is prominent for determining the overall volume that is affected by relativistic jet injection.

Massive stars shape their surrounding medium through the force of their stellar winds, which collide with the circumstellar medium. Because the characteristics of these stellar winds vary over the course of the evolution of the star, the circumstellar matter becomes a reflection of the stellar evolution and can be used to determine the characteristics of the progenitor star. In particular, whenever a fast wind phase follows a slow wind phase, the fast wind sweeps up its predecessor in a shell, which is observed as a circumstellar nebula. We make 2-D and 3-D numerical simulations of fast stellar winds sweeping up their slow predecessors to investigate whether numerical models of these shells have to be 3-D, or whether 2-D models are sufficient to reproduce the shells correctly. We focus on those situations where a fast Wolf-Rayet (WR) star wind sweeps up the slower wind emitted by its predecessor, being either a red supergiant or a luminous blue variable. As the fast WR wind expands, it creates a dense shell of swept up material that expands outward, driven by the high pressure of the shocked WR wind. These shells are subject to a fair variety of hydrodynamic-radiative instabilities. If the WR wind is expanding into the wind of a luminous blue variable phase, the instabilities will tend to form a fairly small-scale, regular filamentary lattice with thin filaments connecting knotty features. If the WR wind is sweeping up a red supergiant wind, the instabilities will form larger interconnected structures with less regularity. Our results show that 3-D models, when translated to observed morphologies, give realistic results that can be compared directly to observations. The 3-D structure of the nebula will help to distinguish different progenitor scenarios.

In the context of planet formation, anticyclonic vortices have recently received lots of attention for the role they can play in planetesimals formation. Radial migration of intermediate size solids toward the central star may prevent their growth to larger solid grains. On the other hand, vortices can trap the dust and accelerate this growth, counteracting fast radial transport. Multiple effects have been shown to affect this scenario, such as vortex migration or decay. The aim of this paper is to study the formation of vortices by the Rossby wave instability and their long term evolution in a full three dimensional protoplanetary disc. We use a robust numerical scheme combined with adaptive mesh refinement in cylindrical coordinates, allowing to affordably compute long term 3D evolutions. We consider a full disc stratified both radially and vertically that is prone to formation of vortices by the Rossby wave instability. We show that the 3D Rossby vortices grow and survive over hundreds of years without migration. The localized overdensity which initiated the instability and vortex formation survives the growth of the Rossby wave instability for very long times. When the vortices are no longer sustained by the Rossby wave instability, their shape changes toward more elliptical vortices. This allows them to survive shear-driven destruction, but they may be prone to elliptical instability and slow decay. When the conditions for growing Rossby wave-related instabilities are maintained in the disc, large-scale vortices can survive over very long timescales and may be able to concentrate solids.

Filament longitudinal oscillations have been observed on the solar disk in H$\alpha$. We intend to find an example of the longitudinal oscillations of a prominence, where the magnetic dip can be seen directly, and examine what is the restoring force of such kind of oscillations. We carry out a multiwavelength data analysis of the active region prominence oscillations above the western limb on 2007 February 8. Besides, we perform a one-dimensional hydrodynamic simulation of the longitudinal oscillations. The high-resolution observations by Hinode/SOT indicate that the prominence, seen as a concave-inward shape in lower-resolution Extreme Ultraviolet (EUV) images, actually consists of many concave-outward threads, which is indicative of the existence of magnetic dips. After being injected into the dip region, a bulk of prominence material started to oscillate for more than 3.5 hours, with the period being 52 min. The oscillation decayed with time, with the decay timescale being 133 min. Our hydrodynamic simulation can well reproduce the oscillation period, but the damping timescale in the simulation is 1.5 times as long as the observations. The results clearly show the prominence longitudinal oscillations around the dip of the prominence and our study suggests that the restoring force of the longitudinal oscillations might be the gravity. Radiation and heat conduction are insufficient to explain the decay of the oscillations. Other mechanisms, such as wave leakage and mass accretion, have to be considered. The possible relation between the longitudinal oscillations and the later eruption of a prominence thread, as well as a coronal mass ejection (CME), is also discussed.

Starting from a realistically sheared magnetic arcade connecting chromospheric, transition region to coronal plasma, we simulate the in-situ formation and sustained growth of a quiescent prominence in the solar corona. Contrary to previous works, our model captures all phases of the prominence formation, including the loss of thermal equilibrium, its successive growth in height and width to macroscopic dimensions, and the gradual bending of the arched loops into dipped loops, as a result of the mass accumulation. Our 2.5-dimensional, fully thermodynamically and magnetohydrodynamically consistent model mimics the magnetic topology of normal-polarity prominences above a photospheric neutral line, and results in a curtain-like prominence above the neutral line through which the ultimately dipped magnetic field lines protrude at a finite angle. The formation results from concentrated heating in the chromosphere, followed by plasma evaporation and later rapid condensation in the corona due to thermal instability, as verified by linear instability criteria. Concentrated heating in the lower atmosphere evaporates plasma from below to accumulate at the top of coronal loops and supply mass to the later prominence constantly. This is the first evaporation-condensation model study where we can demonstrate how the formed prominence stays in a force balanced state, which can be compared to the Kippenhahn-Schluter type magnetohydrostatic model, all in a finite low-beta corona.

We present numerical simulations of the hydrodynamical interactions that produce circumstellar shells. These simulations include several scenarios, such as wind-wind interaction and wind-ISM collisions. In our calculations we have taken into account the presence of dust in the stellar wind. Our results show that, while small dust grains tend to be strongly coupled to the gas, large dust grains are only weakly coupled. As a result, the distribution of the large dust grains is not representative of the gas distribution. Combining these results with observations may give us a new way of validating hydrodynamical models of the circumstellar medium.

Massive stars lose a large fraction of their original mass over the course of their evolution. These stellar winds shape the surrounding medium according to parameters that are the result of the characteristics of the stars, varying over time as the stars evolve, leading to both permanent and temporary features that can be used to constrain the evolution of the progenitor star. Because long Gamma-Ray Bursts (GRBs) are thought to originate from massive stars, the characteristics of the circumstellar medium (CSM) should be observable in the signal of GRBs. This can occur directly, as the characteristics of the GRB-jet are changed by the medium it collides with, and indirectly because the GRB can only be observed through the extended circumstellar bubble that surrounds each massive star. We use computer simulations to describe the circumstellar features that can be found in the vicinity of massive stars and discuss if, and how, they may appear in GRB observations. Specifically, we make hydrodynamical models of the circumstellar environment of a rapidly rotating, chemically near-homogeneous star, which is represents a GRB progenitor candidate model. The simulations show that the star creates a large scale bubble of shocked wind material, which sweeps up the interstellar medium in an expanding shell. Within this bubble, temporary circumstellar shells, clumps and voids are created as a result of changes in the stellar wind. Most of these temporary features have disappeared by the time the star reaches the end of its life, leaving a highly turbulent circumstellar bubble behind. Placing the same star in a high density environment simplifies the evolution of the CSM as the more confined bubble prohibits the formation of some of the temporary structures.

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.

Starting from accurate MHD flux rope equilibria containing prominence condensations, we initiate a systematic survey of their linear eigenoscillations. To quantify the full spectrum of linear MHD eigenmodes, we require knowledge of all flux-surface localized modes, charting out the continuous parts of the MHD spectrum. We combine analytical and numerical findings for the continuous spectrum for realistic prominence configurations. The equations governing all eigenmodes for translationally symmetric, gravitating equilibria containing an axial shear flow, are analyzed, along with their flux-surface localized limit. The analysis is valid for general 2.5D equilibria, where either density, entropy, or temperature vary from one flux surface to another. We analyze the mode couplings caused by the poloidal variation in the flux rope equilibria, by performing a small gravity parameter expansion. We contrast the analytical results with continuous spectra obtained numerically. For equilibria where the density is a flux function, we show that continuum modes can be overstable, and we present the stability criterion for these convective continuum instabilities. Furthermore, for all equilibria, a four-mode coupling scheme between an Alfvenic mode of poloidal mode number m and three neighboring (m-1, m, m+1) slow modes is identified, occurring in the vicinity of rational flux surfaces. For realistically prominence equilibria, this coupling is shown to play an important role, from weak to stronger gravity parameter g values. The analytic predictions for small g are compared with numerical spectra, and progressive deviations for larger g are identified. The unstable continuum modes could be relevant for short-lived prominence configurations. The gaps created by poloidal mode coupling in the continuous spectrum need further analysis, as they form preferred frequency ranges for global eigenoscillations.

It has been established that cold plasma condensations can form in a magnetic loop subject to localized heating of the footpoints. In this paper, we use grid-adaptive numerical simulations of the radiative hydrodynamic equations to parametrically investigate the filament formation process in a pre-shaped loop with both steady and finite-time chromospheric heating. Compared to previous works, we consider low-lying loops with shallow dips, and use a more realistic description for the radiative losses. We demonstrate for the first time that the onset of thermal instability satisfies the linear instability criterion. The onset time of the condensation is roughly \sim 2 hr or more after the localized heating at the footpoint is effective, and the growth rate of the thread length varies from 800 km hr-1 to 4000 km hr-1, depending on the amplitude and the decay length scale characterizing this localized chromospheric heating. We show how single or multiple condensation segments may form in the coronal portion. In the asymmetric heating case, when two segments form, they approach and coalesce, and the coalesced condensation later drains down into the chromosphere. With a steady heating, this process repeats with a periodicity of several hours. While our parametric survey confirms and augments earlier findings, we also point out that steady heating is not necessary to sustain the condensation. Once the condensation is formed, it can keep growing also when the localized heating ceases. Finally, we show that the condensation can survive continuous buffeting by perturbations resulting from the photospheric p-mode waves.

We study the hydrodynamical behavior occurring in the turbulent interaction zone of a fast moving red supergiant star, where the circumstellar and interstellar material collide. In this wind-interstellar medium collision, the familiar bow shock, contact discontinuity, and wind termination shock morphology forms, with localized instability development. Our model includes a detailed treatment of dust grains in the stellar wind, and takes into account the drag forces between dust and gas. The dust is treated as pressureless gas components binned per grainsize, for which we use ten representative grainsize bins. Our simulations allow to deduce how dust grains of varying sizes become distributed throughout the circumstellar medium. We show that smaller dust grains (radius <0.045 micro-meters) tend to be strongly bound to the gas and therefore follow the gas density distribution closely, with intricate finestructure due to essentially hydrodynamical instabilities at the wind-related contact discontinuity. Larger grains which are more resistant to drag forces are shown to have their own unique dust distribution, with progressive deviations from the gas morphology. Specifically, small dust grains stay entirely within the zone bound by shocked wind material. The large grains are capable of leaving the shocked wind layer, and can penetrate into the shocked or even unshocked interstellar medium. Depending on how the number of dust grains varies with grainsize, this should leave a clear imprint in infrared observations of bowshocks of red supergiants and other evolved stars.

Strong optical and radio flares often appear in the afterglow phase of Gamma-Ray Bursts (GRBs). It has been proposed that colliding ultra-relativistic shells can produce these flares. Such consecutive shells can be formed due to the variability in the central source of a GRB. We perform high resolution 1D numerical simulations of late collisions between two ultra-relativistic shells in order to explore these events. We examine the case where a cold uniform shell collides with a self-similar Blandford and McKee shell in a constant density environment and consider cases with different Lorentz factor and energy for the uniform shell. We produce the corresponding on-axis light curves and emission images for the afterglow phase and examine the occurrence of optical and radio flares assuming a spherical explosion and a hard-edged jet scenario. For our simulations we use the Adaptive Mesh Refinement version of the Versatile Advection Code (AMRVAC) coupled to a linear radiative transfer code to calculate synchrotron emission. We find steeply rising flare like behavior for small jet opening angles and more gradual rebrightenings for large opening angles. Synchrotron self-absorption is found to strongly influence the onset and shape of the radio flare.

As massive stars evolve, their winds change. This causes a series of hydrodynamical interactions in the surrounding medium. Whenever a fast wind follows a slow wind phase, the fast wind sweeps up the slow wind in a shell, which can be observed as a circumstellar nebula. One of the most striking examples of such an interaction is when a massive star changes from a red supergiant into a Wolf-Rayet star. Nebulae resulting from such a transition have been observed around many Wolf-Rayet stars and show detailed, complicated structures owing to local instabilities in the swept-up shells. Shells also form in the case of massive binary stars, where the winds of two stars collide with one another. Along the collision front gas piles up, forming a shell that rotates along with the orbital motion of the binary stars. In this case the shell follows the surface along which the ram pressure of the two colliding winds is in balance. Using the MPI-AMRVAC hydrodynamics code we have made multi-dimensional simulations of these interactions in order to model the formation and evolution of these circumstellar nebulae and explore whether full 3D simulations are necessary to obtain accurate models of such nebulae.

Energy loss through optically thin radiative cooling plays an important part in the evolution of astrophysical gas dynamics and should therefore be considered a necessary element in any numerical simulation. Although the addition of this physical process to the equations of hydrodynamics is straightforward, it does create numerical challenges that have to be overcome in order to ensure the physical correctness of the simulation. First, the cooling has to be treated (semi-)implicitly, owing to the discrepancies between the cooling timescale and the typical timesteps of the simulation. Secondly, because of its dependence on a tabulated cooling curve, the introduction of radiative cooling creates the necessity for an interpolation scheme. In particular, we will argue that the addition of radiative cooling to a numerical simulation creates the need for extremely high resolution, which can only be fully met through the use of adaptive mesh refinement.

We investigate the morphology of the collision front between the stellar winds of binary components in two long-period binary systems, one consisting of a hydrogen rich Wolf-Rayet star (WNL) and an O-star and the other of a Luminous Blue Variable (LBV) and an O-star. Specifically, we follow the development and evolution of instabilities that form in such a shell, if it is sufficiently compressed, due to both the wind interaction and the orbital motion. We use MPI-AMRVAC to time-integrate the equations of hydrodynamics, combined with optically thin radiative cooling, on an adaptive mesh 3D grid. Using parameters for generic binary systems, we simulate the interaction between the winds of the two stars. The WNL+O star binary shows a typical example of an adiabatic wind collision. The resulting shell is thick and smooth, showing no instabilities. On the other hand, the shell created by the collision of the O star wind with the LBV wind, combined with the orbital motion of the binary components, is susceptible to thin shell instabilities, which create a highly structured morphology. We identify the nature of the instabilities as both linear and non-linear thin-shell instabilities, with distinct differences between the leading and the trailing parts of the collision front. We also find that for binaries containing a star with a (relatively) slow wind, the global shape of the shell is determined more by the slow wind velocity and the orbital motion of the binary, than the ram pressure balance between the two winds. The interaction between massive binary winds needs further parametric exploration, to identify the role and dynamical importance of multiple instabilities at the collision front, as shown here for an LBV+O star system.

In gamma-ray-bursts (GRB), ultra-relativistic blast waves are ejected into the circumburst medium. We analyse in unprecedented detail the deceleration of a self-similar Blandford-McKee blast wave from a Lorentz factor 25 to the nonrelativistic Sedov phase. Our goal is to determine the stability properties of its frontal shock. We carried out a grid-adaptive relativistic 2D hydro-simulation at extreme resolving power, following the GRB jet during the entire afterglow phase. We investigate the effect of the finite initial jet opening angle on the deceleration of the blast wave, and identify the growth of various instabilities throughout the coasting shock front. We find that during the relativistic phase, the blast wave is subject to pressure-ram pressure instabilities that ripple and fragment the frontal shock. These instabilities manifest themselves in the ultra-relativistic phase alone, remain in full agreement with causality arguments, and decay slowly to finally disappear in the near-Newtonian phase as the shell Lorentz factor drops below 3. From then on, the compression rate decreases to levels predicted to be stable by a linear analysis of the Sedov phase. Our simulations confirm previous findings that the shell also spreads laterally because a rarefaction wave slowly propagates to the jet axis, inducing a clear shell deformation from its initial spherical shape. The blast front becomes meridionally stratified, with decreasing speed from axis to jet edge.

We present a study of the intermediate regime between ultra-relativistic and nonrelativistic flow for gamma-ray burst afterglows. The hydrodynamics of spherically symmetric blast waves is numerically calculated using the AMRVAC adaptive mesh refinement code. Spectra and light curves are calculated using a separate radiation code that, for the first time, links a parametrisation of the microphysics of shock acceleration, synchrotron self-absorption and electron cooling to a high-performance hydrodynamics simulation. For the dynamics we find that the transition to the nonrelativistic regime generally occurs later than expected, that the Sedov-Taylor solution overpredicts the late time blast wave radius and that the analytical formula for the blast wave velocity from Huang (1999) overpredicts the late time velocity by a factor 4/3. For the radiation we find that the flux may differ up to an order of magnitude depending on the equation of state that is used for the fluid and that the counterjet leads to a clear rebrightening at late times for hard-edged jets. Simulating GRB030329 using predictions for its physical parameters from the literature leads to spectra and light curves that may differ significantly from the actual data, emphasizing the need for very accurate modelling. Predicted light curves at low radio frequencies for a hard-edged jet model of GRB030329 with opening angle 22 degrees show typically two distinct peaks, due to the combined effect of jet break, non relativistic break and counterjet.

The conventional derivation of the gamma-ray burst afterglow jet break time uses only the blast wave fluid Lorentz factor and therefore leads to an achromatic break. We show that in general gamma-ray burst afterglow jet breaks are chromatic across the self-absorption break. Depending on circumstances, the radio jet break may be postponed significantly. Using high-accuracy adaptive mesh fluid simulations in one dimension, coupled to a detailed synchrotron radiation code, we demonstrate that this is true even for the standard fireball model and hard-edged jets. We confirm these effects with a simulation in two dimensions. The frequency dependence of the jet break is a result of the angle dependence of the emission, the changing optical depth in the self-absorbed regime and the shape of the synchrotron spectrum in general. In the optically thin case the conventional analysis systematically overestimates the jet break time, leading to inferred opening angles that are underestimated by a factor 1.32 and explosion energies that are underestimated by a factor 1.73, for explosions in a homogeneous environment.

We simulate time-dependent particle acceleration in the blast wave of a young supernova remnant (SNR), using a Monte Carlo approach for the diffusion and acceleration of the particles, coupled to an MHD code. We calculate the distribution function of the cosmic rays concurrently with the hydrodynamic evolution of the SNR, and compare the results with those obtained using simple steady-state models. The surrounding medium into which the supernova remnant evolves turns out to be of great influence on the maximum energy to which particles are accelerated. In particular, a shock going through a $\rho \propto r^{-2}$ density profile causes acceleration to typically much higher energies than a shock going through a medium with a homogeneous density profile. We find systematic differences between steady-state analytical models and our time-dependent calculation in terms of spectral slope, maximum energy, and the shape of the cut-off of the particle spectrum at the highest energies. We also find that, provided that the magnetic field at the reverse shock is sufficiently strong to confine particles, cosmic rays can be easily re-accelerated at the reverse shock.

This work presents a new plasma cooling curve that is calculated using the SPEX package. We compare our cooling rates to those in previous works, and implement the new cooling function in the grid-adaptive framework `AMRVAC'. Contributions to the cooling rate by the individual elements are given, to allow for the creation of cooling curves tailored to specific abundance requirements. In some situations, it is important to be able to include radiative losses in the hydrodynamics. The enhanced compression ratio can trigger instabilities (such as the Vishniac thin-shell instability) that would otherwise be absent. For gas with temperatures below 10,000 K, the cooling time becomes very long and does not affect the gas on the timescales that are generally of interest for hydrodynamical simulations of circumstellar plasmas. However, above this temperature, a significant fraction of the elements is ionised, and the cooling rate increases by a factor 1000 relative to lower temperature plasmas.

Transverse stratification is a common intrinsic feature of astrophysical jets. There is growing evidence that jets in radio galaxies consist of a fast low density outflow at the jet axis, surrounded by a slower, denser, extended jet. The inner and outer jet components then have a different origin and launching mechanism, making their effective inertia, magnetization, associated energy flux and angular momentum content different as well. Their interface will develop differential rotation, where disruptions may occur. We here investigate the stability of rotating, two-component relativistic outflows typical for jets in radio galaxies. For this purpose, we parametrically explore the long term evolution of a transverse cross-section of radially stratified jets numerically With grid-adaptive relativistic MHD simulations, augmented with approximate linear stability analysis. We study the influence of dynamically important poloidal magnetic fields, with varying contributions of the inner component jet to the total kinetic energy flux of the jet, on their non-linear azimuthal stability. We demonstrate that two-component jets with high kinetic energy flux, and an inner jet effective inertia which is higher than the outer jet effective inertia are subject to the development of a relativistically enhanced, rotation-induced Rayleigh-Taylor type instability. This instability plays a major role in decelerating the inner jet and the overall jet decollimation. This novel deceleration scenario can partly explain the radio source dichotomy, relating it directly to the efficiency of the central engine in launching the inner jet component. The FRII/FRI transition could then occur when the relative kinetic energy flux of the inner to the outer jet grows beyond a certain treshold.

Gamma-ray burst afterglow flares and rebrightenings of the optical and X-ray light curve have been attributed to both late time inner engine activity and density changes in the medium surrounding the burster. To test the latter, we study the encounter between the relativistic blast wave from a gamma-ray burster and a stellar wind termination shock. The blast wave is simulated using a high performance adaptive mesh relativistic hydrodynamics code, AMRVAC, and the synchrotron emission is analyzed in detail with a separate radiation code. We find no bump in the resulting light curve, not even for very high density jumps. Furthermore, by analyzing the contributions from the different shock wave regions we are able to establish that it is essential to resolve the blast wave structure in order to make qualitatively correct predictions on the observed output and that the contribution from the reverse shock region will not stand out, even when the magnetic field is increased in this region by repeated shocks. This study resolves a controversy in recent literature.

Observations show that the magnetic field in young supernova remnants (SNRs) is significantly stronger than can be expected from the compression of the circumstellar medium (CSM) by a factor of four expected for strong blast waves. Additionally, the polarization is mainly radial, which is also contrary to expectation from compression of the CSM magnetic field. Cosmic rays (CRs) may help to explain these two observed features. They can increase the compression ratio to factors well over those of regular strong shocks by adding a relativistic plasma component to the pressure, and by draining the shock of energy when CRs escape from the region. The higher compression ratio will also allow for the contact discontinuity, which is subject to the Rayleigh-Taylor (R-T) instability, to reach much further out to the forward shock. This could create a preferred radial polarization of the magnetic field. With an adaptive mesh refinement MHD code (AMRVAC), we simulate the evolution of SNRs with three different configurations of the initial CSM magnetic field, and look at two different equations of state in order to look at the possible influence of a CR plasma component. The spectrum of CRs can be simulated using test particles, of which we also show some preliminary results that agree well with available analytical solutions.

Supernova remnants (SNR) are now widely believed to be a source of cosmic rays (CRs) up to an energy of 1 PeV. The magnetic fields required to accelerate CRs to sufficiently high energies need to be much higher than can result from compression of the circumstellar medium (CSM) by a factor 4, as is the case in strong shocks. Non-thermal synchrotron maps of these regions indicate that indeed the magnetic field is much stronger, and for young SNRs has a dominant radial component while for old SNRs it is mainly toroidal. How these magnetic fields get enhanced, or why the field orientation is mainly radial for young remnants, is not yet fully understood. We use an adaptive mesh refinement MHD code, AMRVAC, to simulate the evolution of supernova remnants and to see if we can reproduce a mainly radial magnetic field in early stages of evolution. We follow the evolution of the SNR with three different configurations of the initial magnetic field in the CSM: an initially mainly toroidal field, a turbulent magnetic field, and a field parallel to the symmetry axis. Although for the latter two topologies a significant radial field component arises at the contact discontinuity due to the Rayleigh-Taylor instability, no radial component can be seen out to the forward shock. Ideal MHD appears not sufficient to explain observations. Possibly a higher compression ratio and additional turbulence due to dominant presence of CRs can help us to better reproduce the observations in future studies.

The properties of linear Alfv\'en, slow, and fast magnetoacoustic waves for uniform plasmas in relativistic magnetohydrodynamics (MHD) are discussed, augmenting the well-known expressions for their phase speeds with knowledge on the group speed. A 3+1 formalism is purposely adopted to make direct comparison with the Newtonian MHD limits easier and to stress the graphical representation of their anisotropic linear wave properties using the phase and group speed diagrams. By drawing these for both the fluid rest frame and for a laboratory Lorentzian frame which sees the plasma move with a three-velocity having an arbitrary orientation with respect to the magnetic field, a graphical view of the relativistic aberration effects is obtained for all three MHD wave families. Moreover, it is confirmed that the classical Huygens construction relates the phase and group speed diagram in the usual way, even for the lab frame viewpoint. Since the group speed diagrams correspond to exact solutions for initial conditions corresponding to a localized point perturbation, their formulae and geometrical construction can serve to benchmark current high-resolution algorithms for numerical relativistic MHD.

The deceleration mechanisms for relativistic jets in active galactic nuclei remain an open question, and in this paper we propose a model which could explain sudden jet deceleration, invoking density discontinuities. This is particularly motivated by recent indications from HYMORS. Exploiting high resolution, numerical simulations, we demonstrate that for both high and low energy jets, always at high Lorentz factor, a transition to a higher density environment can cause a significant fraction of the directed jet energy to be lost on reflection. This can explain how one-sided jet deceleration and a transition to FR I type can occur in HYMORS, which start as FR II (and remain so on the other side). For that purpose, we implemented in the relativistic hydrodynamic grid-adaptive AMRVAC code, the Synge-type equation of state introduced in the general polytropic case by Meliani et al. (2004). We present results for 10 model computations, varying the inlet Lorentz factor from 10 to 20, including uniform or decreasing density profiles, and allowing for cylindrical versus conical jet models. As long as the jet propagates through uniform media, we find that the density contrast sets most of the propagation characteristics, fully consistent with previous modeling efforts. When the jet runs into a denser medium, we find a clear distinction in the decelaration of high energy jets depending on the encountered density jump. For fairly high density contrast, the jet becomes destabilised and compressed, decelerates strongly (up to subrelativistic speeds) and can form knots. We point out differences that are found between cylindrical and conical jet models, together with dynamical details like the Richtmyer-Meshkov instabilities developing at the original contact interface.

A robust second order, shock-capturing numerical scheme for multi-dimensional special relativistic magnetohydrodynamics on computational domains with adaptive mesh refinement is presented. The base solver is a total variation diminishing Lax-Friedrichs scheme in a finite volume setting and is combined with a diffusive approach for controlling magnetic monopole errors. The consistency between the primitive and conservative variables is ensured at all limited reconstructions and the spatial part of the four velocity is used as a primitive variable. Demonstrative relativistic examples are shown to validate the implementation. We recover known exact solutions to relativistic MHD Riemann problems, and simulate the shock-dominated long term evolution of Lorentz factor 7 vortical flows distorting magnetic island chains.

Extragalactic jets are inferred to harbor dynamically important, organized magnetic fields which presumably aid in the collimation of the relativistic jet flows. We here explore by means of grid-adaptive, high resolution numerical simulations the morphology of AGN jets pervaded by helical field and flow topologies. We concentrate on morphological features of the bow shock and the jet beam behind the Mach disk, for various jet Lorentz factors and magnetic field helicities. We investigate the influence of helical magnetic fields on jet beam propagation in overdense external medium. We use the AMRVAC code, employing a novel hybrid block-based AMR strategy, to compute ideal plasma dynamics in special relativity. The helicity of the beam magnetic field is effectively transported down the beam, with compression zones in between diagonal internal cross-shocks showing stronger toroidal field regions. In comparison with equivalent low-relativistic jets which get surrounded by cocoons with vortical backflows filled by mainly toroidal field, the high speed jets demonstrate only localized, strong toroidal field zones within the backflow vortical structures. We find evidence for a more poloidal, straight field layer, compressed between jet beam and backflows. This layer decreases the destabilizing influence of the backflow on the jet beam. In all cases, the jet beam contains rich cross-shock patterns, across which part of the kinetic energy gets transferred. For the high speed reference jet considered here, significant jet deceleration only occurs beyond distances exceeding ${\cal O}(100 R_j)$, as the axial flow can reaccelerate downstream to the internal cross-shocks. This reacceleration is magnetically aided, due to field compression across the internal shocks which pinch the flow.

Aims : We re-examine the conditions required to steadily deviate an accretion flow from a circumstellar disc into a magnetospheric funnel flow onto a slow rotating young forming star. Methods : New analytical constraints on the formation of accretion funnels flows due to the presence of a dipolar stellar magnetic field disrupting the disc are derived. The Versatile Advection Code is used to confirm these constraints numerically. Axisymmetric MHD simulations are performed, where a stellar dipole field enters the resistive accretion disc, whose structure is self-consistently computed. Results : The analytical criterion derived allows to predict a priori the position of the truncation radius from a non perturbative accretion disc model. Accretion funnels are found to be robust features which occur below the co-rotation radius, where the stellar poloidal magnetic pressure becomes both at equipartition with the disc thermal pressure and is comparable to the disc poloidal ram pressure. We confirm the results of Romanova et al. 2002 and find accretion funnels for stellar dipole fields as low as 140 G in the low accretion rate limit of $10^{-9} M_\odot.yr^{-1}$. With our present numerical setup with no disc magnetic field, we found no evidence of winds, neither disc driven nor X-winds, and the star is only spun up by its interaction with the disc. Conclusions : Weak dipole fields, similar in magnitude to those observed, lead to the development of accretion funnel flows in weakly accreting T Tauri stars. However, the higher accretion observed for most T Tauri stars (${\dot M} \sim 10^{-8} M_\odot.yr^{-1}$) requires either larger stellar field strength and/or different magnetic topologies to allow for magnetospheric accretion.

Context: Astrophysical jets from various sources seem to be stratified, with a fast inner jet and a slower outer jet. As it is likely that the launching mechanism for each component is different, their interface will develop differential rotation, while the outer jet radius represents a second interface where disruptions may occur. Aims: We explore the stability of stratified, rotating, relativistic two-component jets, in turn embedded in static interstellar medium. Methods: In a grid-adaptive relativistic hydrodynamic simulation with the AMRVAC code, the non-linear azimuthal stability of two-component relativistic jets is investigated. We simulate until multiple inner jet rotations have been completed. Results: We find evidence for the development of an extended shear flow layer between the two jet components, resulting from the growth of a body mode in the inner jet, Kelvin-Helmholtz surface modes at their original interface, and their nonlinear interaction. Both wave modes are excited by acoustic waves which are reflected between the symmetry axis and the interface of the two jet components. Their interaction induces the growth of near stationary, counterrotating vortices at the outer edge of the shear flow layer. The presence of a heavy external jet allows to slow down their further development, and maintain a collimated flow. At the outer jet boundary, small-scale Rayleigh-Taylor instabilities develop, without disrupting the jet configuration. Conclusion: We demonstrate that the cross-section of two-component relativistic jets, with a heavy, cold outer jet, is non-linearly stable.

We present numerical magnetohydrostatic solutions describing the gravitationally stratified, bulk equilibrium of cool, dense prominence plasma embedded in a near-potential coronal field. These solutions are calculated using the FINESSE magnetohydrodynamics equilibrium solver and describe the morphologies of magnetic field distributions in and around prominences and the cool prominence plasma that these fields support. The equilibrium condition for this class of problem is usually different in distinct subdomains, separated by free boundaries, across which solutions are matched by suitable continuity or jump conditions describing force balance. We employ our precise finite element elliptic solver to calculate solutions not accessible by previous analytical techniques with temperature or entropy prescribed as free functions of the magnetic flux function, including a range of values of the polytropic index, temperature variations mainly across magnetic field lines and photospheric field profiles sheared close to the polarity inversion line. Out of the many examples computed here, perhaps the most noteworthy is one which reproduces precisely the three-part structure often encountered in observations: a cool dense prominence within a cavity/flux rope embedded in a hot corona. The stability properties of these new equilibria, which may be relevant to solar eruptions, can be determined in the form of a full resistive MHD spectrum using a companion hyperbolic stability solver.

Context:A significant fraction of progenitors for long gamma-ray bursts (GRBs) are believed to be massive stars. The investigation of long GRBs therefore requires modeling the propagation of ultra-relativistic blastwaves through the circumburst medium surrounding massive stars. We simulate the expansion of an isotropic, adiabatic relativistic fireball into the wind-shaped medium around a massive GRB progenitor. The circumburst medium is composed of a realistically stratified stellar wind zone up to its termination shock, followed by a region of shocked wind characterized by a constant density. Aims: We followed the evolution of the blastwave through all its stages, including the extremely rapid acceleration up to a Lorentz factor 75 flow, its deceleration by interaction with stellar wind, its passage of the wind termination shock, until its propagation through shocked wind. Methods: We used the adaptive mesh refinement versatile advection code to follow the evolution of the fireball. Results: We show that the acceleration from purely thermal to ultra-relativistic kinetic regimes is abrupt and produces an internally structured blastwave. We resolved the structure of this ultra-relativistic shell in all stages, thanks to the adaptive mesh. We comment on the dynamical roles played by forward and reverse shock pairs in the phase of interaction with the free stellar wind and clearly identify the complex shock-dominated structure created when the shell crosses the terminal shock. Conclusion: We show that in our model where the terminal shock is taken relatively close to the massive star, the phase of self-similar deceleration of Blandford-McKee type can only be produced in the constant density, shocked wind zone.

We present a detailed study of localised magnetohydrodynamical (MHD) instabilities occuring in two--dimensional magnetized accretion disks. We model axisymmetric MHD disk tori, and solve the equations governing a two--dimensional magnetized accretion disk equilibrium and linear wave modes about this equilibrium. We show the existence of novel MHD instabilities in these two--dimensional equilibria which do not occur in an accretion disk in the cylindrical limit. The disk equilibria are numerically computed by the FINESSE code. The stability of accretion disks is investigated analytically as well as numerically. We use the PHOENIX code to compute all the waves and instabilities accessible to the computed disk equilibrium. We concentrate on strongly magnetized disks and sub--Keplerian rotation in a large part of the disk. These disk equilibria show that the thermal pressure of the disk can only decrease outwards if there is a strong gravitational potential. Our theoretical stability analysis shows that convective continuum instabilities can only appear if the density contours coincide with the poloidal magnetic flux contours. Our numerical results confirm and complement this theoretical analysis. Furthermore, these results show that the influence of gravity can either be stabilizing or destabilizing on this new kind of MHD instability. In the likely case of a non--constant density, the height of the disk should exceed a threshold before this type of instability can play a role. This localised MHD instability provides an ideal, linear route to MHD turbulence in strongly magnetized accretion disk tori.

We apply a novel adaptive mesh refinement code, AMRVAC, to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold Interstellar Medium (ISM). We do this for both 1D isotropic as well as full 2D jetlike fireball models. This is relevant for Gamma Ray Bursts, and we demonstrate that, thanks to the AMR strategy, we resolve the internal structure of the shocked shell-ISM matter, which will leave its imprint on the GRB afterglow. We determine the deceleration from an initial Lorentz factor $\gamma=100$ up to the almost Newtonian $\gamma\sim{\cal O}(2)$ phase of the flow. We present axisymmetric 2D shell evolutions, with the 2D extent characterized by their initial opening angle. In such jetlike GRB models, we discuss the differences with the 1D isotropic GRB equivalents. These are mainly due to thermally induced sideways expansions of both the shocked shell and shocked ISM regions. We found that the propagating 2D ultrarelativistic shell does not accrete all the surrounding medium located within its initial opening angle. Part of this ISM matter gets pushed away laterally and forms a wide bow-shock configuration with swirling flow patterns trailing the thin shell. The resulting shell deceleration is quite different from that found in isotropic GRB models. As long as the lateral shell expansion is merely due to ballistic spreading of the shell, isotropic and 2D models agree perfectly. As thermally induced expansions eventually lead to significantly higher lateral speeds, the 2D shell interacts with comparably more ISM matter and decelerates earlier than its isotropic counterpart.

We present analytical and numerical studies of magnetorotational instabilities occuring in magnetized accretion disks. In these studies we make use of the linearised compressible MHD equations. These calculations are performed for general radially stratified disks in the cylindrical limit. In particular, we investigate the influence of nonvanishing toroidal magnetic field component on the growth rate and oscillation frequency of magnetorotational instabilities in Keplerian disks. We find the persistence of these instabilities in accretion disks close to equipartition. Our calculations show that these eigenmodes become overstable (complex eigenvalue), due to the presence of a toroidal magnetic field component, while their growth rate reduces slightly. Furthermore, we demonstrate the presence of magneto-rotational overstabilities in weakly magnetized sub-Keplerian rotating disks. We show that the growth rate scales with the rotation frequency of the disk. These eigenmodes also have a nonzero oscillation frequency, due to the presence of the dominant toroidal magnetic field component. The overstable character of the MRI increases as the rotation frequency of the disk decreases.

We present a study of instabilities occuring in thick magnetized accretion disks. We calculate the growth rates of these instabilities and characterise precisely the contribution of the magneto-rotational and the convective mechanism. All our calculations are performed in radially stratified disks in the cylindrical limit. The numerical calculations are performed using the appropriate local dispersion equation solver discussed in Blokland et al. (2005). A comparison with recent results by Narayan et al. (2002) shows excellent agreement with their approximate growth rates only if the disks are weakly magnetized. However, for disks close to equipartition, the dispersion equation from Narayan et al. (2002) loses its validity. Our calculations allow for a quantitative determination of the increase of the growth rate due to the magneto-rotational mechanism. We find that the increase of the growth rate for long wavelength convective modes caused by this mechanism is almost neglible. On the other hand, the growth rate of short wavelength instabilities can be significantly increased by this mechanism, reaching values up to 60%.

In two previous publications$^{1,2}$, we have demonstrated that stationary rotation of magnetized plasma about a compact central object permits an enormous number of different MHD instabilities, with the well-known magneto-rotational instability as just one of them. We here concentrate on the new instabilities found that are driven by transonic transitions of the poloidal flow. A particularly promising class of instabilities, from the point of view of MHD turbulence in accretion disks, is the class of {\em trans-slow Alfven continuum modes}, that occur when the poloidal flow exceeds a critical value of the slow magnetosonic speed. When this happens, virtually every magnetic/flow surface of the disk becomes unstable with respect to highly localized modes of the continuous spectrum. The mode structures rotate, in turn, about the rotating disk. These structure lock and become explosively unstable when the mass of the central object is increased beyond a certain critical value. Their growth rates then become huge, of the order of the Alfven transit time. These instabilities appear to have all requisite properties to facilitate accretion flows across magnetic surfaces and jet formation.[1] R. Keppens, F. Casse, J.P. Goedbloed, "Waves and instabilities in accretion disks: Magnetohydrodynamic spectroscopic analysis", Astrophys. J. {\bf 569}, L121--L126 (2002).[2] J.P. Goedbloed, A.J.C. Belien, B. van der Holst, R. Keppens, "Unstable continuous spectra of transonic axisymmetric plasmas", Phys. Plasmas {\bf 11}, 28--54 (2004).

We present a numerical study of an idealized magnetohydrodynamic (MHD) configuration consisting of a planar wake flow embedded into a three-dimensional (3D) sheared magnetic field. Our simulations investigate the possibility for in-situ development of large-scale compressive disturbances at cospatial current sheet -- velocity shear regions in the heliosphere. Using a linear MHD solver, we first systematically chart the destabilized wavenumbers, corresponding growth rates, and physical parameter ranges for dominant 3D sinuous-type instabilities in an equilibrium wake--current sheet system. Wakes bounded by sufficiently supersonic (Mach number $M_s > 2.6$) flow streams are found to support dominant fully 3D sinuous instabilities when the plasma beta is of order unity. Fully nonlinear, compressible 2.5D and 3D MHD simulations show the self-consistent formation of shock fronts of fast magnetosonic type. They carry density perturbations far away from the wake's center. Shock formation conditions are identified in sonic and Alfv\'enic Mach number parameter space. Depending on the wake velocity contrast and magnetic field magnitude, as well as on the initial perturbation, the emerging shock patterns can be plane-parallel as well as fully three-dimensionally structured. Similar large-scale transients could therefore originate at distances far above coronal helmet streamers or at the location of the ecliptic current sheet.

Obtainable computational efficiency is evaluated when using an Adaptive Mesh Refinement (AMR) strategy in time accurate simulations governed by sets of conservation laws. For a variety of 1D, 2D, and 3D hydro- and magnetohydrodynamic simulations, AMR is used in combination with several shock-capturing, conservative discretization schemes. Solution accuracy and execution times are compared with static grid simulations at the corresponding high resolution and time spent on AMR overhead is reported. Our examples reach corresponding efficiencies of 5 to 20 in multidimensional calculations and only 1.5 -- 8 % overhead is observed. For AMR calculations of multi-dimensional magnetohydrodynamic problems, several strategies for controlling the $\nabla \cdot \BB =0$ constraint are examined. Three source term approaches suitable for cell-centered $\BB$ representations are shown to be effective. For 2D and 3D calculations where a transition to a more globally turbulent state takes place, it is advocated to use an approximate Riemann solver based discretization at the highest allowed level(s), in combination with the robust Total Variation Diminishing Lax-Friedrichs method on the coarser levels. This level-dependent use of the spatial discretization acts as a computationally efficient, hybrid scheme.

The Kelvin-Helmholtz (KH) instability occurring in a single shear flow configuration that is embedded in a uniform flow-aligned magnetic field, is revisited by means of high resolution two-dimensional (2D) magnetohydrodynamic (MHD) simulations. First, the calculations extend previous studies of magnetized shear flows to a higher compressibility regime. The nonlinear evolution of an isolated KH billow emerging from the fastest growing linear mode for a convective sonic Mach number $M_{cs}=0.7$ layer is in many respects similar to its less compressible counterpart (Mach $M_{cs}=0.5$). In particular, the disruptive regime where locally amplified, initially weak magnetic fields, control the nonlinear saturation process is found for Alfv\'en Mach numbers $4\simlt M_A \simlt 30$. The most notable difference between $M_{cs}=0.7$ versus $M_{cs}=0.5$ layers is that higher density contrasts and fast magnetosonic shocklet structures are observed. Second, the use of adaptive mesh refinement allows to parametrically explore much larger computational domains, including up to 22 wavelengths of the linearly dominant mode. A strong process of large-scale coalescence is found, whatever the magnetic field regime. It proceeds through continuous pairing/merging events between adjacent vortices up to the point where the final large-scale vortical structure reaches the domain dimensions. This pairing/merging process is attributed to the growth of subharmonic modes and is mainly controlled by relative phase differences between them. These grid-adaptive simulations demonstrate that even in very weak magnetic field regimes ($M_A \simeq 30$), the large-scale KH coalescence process can trigger tearing-type reconnection events previously identified in cospatial current-vortex sheets.

We present magnetohydrodynamic simulations of a resistive accretion disk continuously launching transmagnetosonic, collimated jets. We time-evolve the full set of magnetohydrodynamic equations, but neglect radiative losses in the energetics (radiatively inefficient). Our calculations demonstrate that a jet is self-consistently produced by the interaction of an accretion disk with an open, initially bent large-scale magnetic field. A constant fraction of heated disk material is launched in the inner equipartition disk regions, leading to the formation of a hot corona and a bright collimated, super-fastmagnetosonic jet. We illustrate the complete dynamics of the ``hot'' near steady-state outflow (where thermal pressure $\simeq$ magnetic pressure) by showing force balance, energy budget and current circuits. The evolution to this near stationary state is analyzed in terms of the temporal variation of energy fluxes controlling the energetics of the accretion disk. We find that unlike advection-dominated accretion flow, the energy released by accretion is mainly sent into the jet rather than transformed into disk enthalpy. These magnetized, radiatively inefficient accretion-ejection structures can account for under-luminous thin disks supporting bright fast collimated jets as seen in many systems displaying jets (for instance M87).

We present numerical magnetohydrodynamic (MHD) simulations of a magnetized accretion disk launching trans-Alfvenic jets. These simulations, performed in a 2.5 dimensional time-dependent polytropic resistive MHD framework, model a resistive accretion disk threaded by an initial vertical magnetic field. The resistivity is only important inside the disk, and is prescribed as eta = alpha_m V_AH exp(-2Z^2/H^2), where V_A stands for Alfven speed, H is the disk scale height and the coefficient alpha_m is smaller than unity. By performing the simulations over several tens of dynamical disk timescales, we show that the launching of a collimated outflow occurs self-consistently and the ejection of matter is continuous and quasi-stationary. These are the first ever simulations of resistive accretion disks launching non-transient ideal MHD jets. Roughly 15% of accreted mass is persistently ejected. This outflow is safely characterized as a jet since the flow becomes super-fastmagnetosonic, well-collimated and reaches a quasi-stationary state. We present a complete illustration and explanation of the `accretion-ejection' mechanism that leads to jet formation from a magnetized accretion disk. In particular, the magnetic torque inside the disk brakes the matter azimuthally and allows for accretion, while it is responsible for an effective magneto-centrifugal acceleration in the jet. As such, the magnetic field channels the disk angular momentum and powers the jet acceleration and collimation. The jet originates from the inner disk region where equipartition between thermal and magnetic forces is achieved. A hollow, super-fastmagnetosonic shell of dense material is the natural outcome of the inwards advection of a primordial field.

A complete analytical and numerical treatment of all magnetohydrodynamic waves and instabilities for radially stratified, magnetized accretion disks is presented. The instabilities are a possible source of anomalous transport. While recovering results on known hydrodynamicand both weak- and strong-field magnetohydrodynamic perturbations, the full magnetohydrodynamic spectra for a realistic accretion disk model demonstrates a much richer variety of instabilities accessible to the plasma than previously realized. We show that both weakly and strongly magnetized accretion disks are prone to strong non-axisymmetric instabilities.The ability to characterize all waves arising in accretion disks holds great promise for magnetohydrodynamic spectroscopic analysis.

Hydrodynamical simulations are presented of a pulsar wind emitted by a supersonically moving pulsar. The pulsar moves through the interstellar medium or, in the more interesting case, through the supernova remnant created at its birth event. In both cases there exists a three-fold structure consisting of the wind termination shock, contact discontinuity and a bow shock bounding the pulsar wind nebula. Using hydrodynamical simulations we study the behaviour of the pulsar wind nebula inside a supernova remnant, and in particular the interaction with the outer shell of swept up interstellar matter and the blast wave surrounding the remnant. This interaction occurs when the pulsar breaks out of the supernova remnant. We assume the remnant is in the Sedov stage of its evolution. Just before break-through, the Mach number associated with the pulsar motion equals ${\cal M}_{\rm psr} = 7/\sqrt{5}$, {\em independent} of the supernova explosion energy and pulsar velocity. The bow shock structure is shown to survive this break-through event.

Axisymmetric stellar wind solutions are presented, obtained by numerically solving the ideal magnetohydrodynamic (MHD) equations. Stationary solutions are critically analysed using the knowledge of the flux functions. These flux functions enter in the general variational principle governing all axisymmetric stationary ideal MHD equilibria. The magnetized wind solutions for (differentially) rotating stars contain both a `wind' and a `dead' zone. We illustrate the influence of the magnetic field topology on the wind acceleration pattern, by varying the coronal field strength and the extent of the dead zone. This is evident from the resulting variations in the location and appearance of the critical curves where the wind speed equals the slow, Alfven, and fast speed. Larger dead zones cause effective, fairly isotropic acceleration to super-Alfvenic velocities as the polar, open field lines are forced to fan out rapidly with radial distance. A higher field strength moves the Alfven transition outwards. In the ecliptic, the wind outflow is clearly modulated by the extent of the dead zone. The combined effect of a fast stellar rotation and an equatorial `dead' zone in a bipolar field configuration can lead to efficient thermo-centrifugal equatorial winds. Such winds show both a strong poleward collimation and some equatorward streamline bending due to significant toroidal field pressure at mid-latitudes. We discuss how coronal mass ejections are then simulated on top of the transonic outflows.

A numerical study of the Kelvin-Helmholtz instability in compressible magnetohydrodynamics is presented. The three-dimensional simulations consider shear flow in a cylindrical jet configuration, embedded in a uniform magnetic field directed along the jet axis. The growth of linear perturbations at specified poloidal and axial mode numbers demonstrate intricate non-linear coupling effects. The physical mechanims leading to induced secondary Kelvin-Helmholtz instabilities at higher mode numbers are identified. The initially weak magnetic field becomes locally dominant in the non-linear dynamics before and during saturation. Thereby, it controls the jet deformation and eventual breakup. The results are obtained using the Versatile Advection Code [G. Toth, Astrophys. Lett. Comm. 34, 245 (1996)], a software package designed to solve general systems of conservation laws. An independent calculation of the same Kelvin-Helmholtz unstable jet configuration using a three-dimensional pseudo-spectral code gives important insights into the coupling and excitation events of the various linear mode numbers.

We discuss steady-state transonic outflows obtained by direct numerical solution of the hydrodynamic and magnetohydrodynamic equations. We make use of the Versatile Advection Code, a software package for solving systems of (hyperbolic) partial differential equations. We proceed stepwise from a spherically symmetric, isothermal, unmagnetized, non-rotating Parker wind to arrive at axisymmetric, polytropic, magnetized, rotating models. These represent 2D generalisations of the analytical 1D Weber-Davis wind solution, which we obtain in the process. Axisymmetric wind solutions containing both a `wind' and a `dead' zone are presented. Since we are solving for steady-state solutions, we efficiently exploit fully implicit time stepping. The method allows us to model thermally and/or magneto-centrifugally driven stellar outflows. We particularly emphasize the boundary conditions imposed at the stellar surface. For these axisymmetric, steady-state solutions, we can use the knowledge of the flux functions to verify the physical correctness of the numerical solutions.

We investigate the Kelvin-Helmholtz instability occuring at the interface of a shear flow configuration in 2D compressible magnetohydrodynamics (MHD). The linear growth and the subsequent non-linear saturation of the instability are studied numerically. We consider an initial magnetic field aligned with the shear flow, and analyze the differences between cases where the initial field is unidirectional everywhere (uniform case), and where the field changes sign at the interface (reversed case). We recover and extend known results for pure hydrodynamic and MHD cases with a discussion of the dependence of the non-linear saturation on the wavenumber, the sound Mach number, and the Alfvenic Mach number for the MHD case. A reversed field acts to destabilize the linear phase of the Kelvin-Helmholtz instability compared to the pure hydrodynamic case, while a uniform field suppresses its growth. In resistive MHD, reconnection events almost instantly accelerate the buildup of a global plasma circulation. They play an important role throughout the further non-linear evolution as well, since the initial current sheet gets amplified by the vortex flow and can become unstable to tearing instabilities forming magnetic islands. As a result, the saturation behaviour and the overall evolution of the density and the magnetic field is markedly different for the uniform versus the reversed field case.