BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Rony Keppens

Rony Keppens appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.HEastro-ph.GAastro-ph.IMphysics.plasm-phastro-ph.EPphysics.comp-ph

Catalog footprint

What is connected

34works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2022arXiv

Legolas: magnetohydrodynamic spectroscopy with viscosity and Hall current

Many linear stability aspects in plasmas are heavily influenced by non-ideal effects beyond the basic ideal magnetohydrodynamics (MHD) description. Here, the extension of the modern open-source MHD spectroscopy code Legolas with viscosity and the Hall current is highlighted and benchmarked on a stringent set of historic and recent findings. The viscosity extension is demonstrated in a cylindrical setup featuring Taylor-Couette flow and in a viscoresistive plasma slab with a tearing mode. For the Hall extension, we show how the full eigenmode spectrum relates to the analytic dispersion relation in an infinite homogeneous medium. We quantify the Hall term influence on the resistive tearing mode in a Harris current sheet, including the effect of compressibility, which is absent in earlier studies. Furthermore, we illustrate how Legolas mimics the incompressible limit easily to compare to literature results. Going beyond published findings, we emphasise the importance of computing the full eigenmode spectrum, and how elements of the spectrum are modified by compressibility. These extensions allow for future stability studies with Legolas that are relevant to ongoing dynamo experiments, protoplanetary disks, or magnetic reconnection.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Plasmoid-fed prominence formation (PF$^2$) during flux rope eruption

We report a new, plasmoid-fed scenario for the formation of an eruptive prominence (PF$^2$), involving reconnection and condensation. We use grid-adaptive resistive two-and-a-half-dimensional magnetohydrodynamic (MHD) simulations in a chromosphere-to-corona setup to resolve this plasmoid-fed scenario. We study a pre-existing flux rope (FR) in the low corona that suddenly erupts due to catastrophe, which also drives a fast shock above the erupting FR. A current sheet (CS) forms underneath the erupting FR, with chromospheric matter squeezed into it. The plasmoid instability occurs and multiple magnetic islands appear in the CS once the Lundquist number reaches $\sim 3.5\times 10^{4}$. The remnant chromospheric matter in the CS is then transferred to the FR by these newly formed magnetic islands. The dense and cool mass transported by the islands accumulates in the bottom of the FR, thereby forming a prominence during the eruption phase. More coronal plasma continuously condenses into the prominence due to the thermal instability as the FR rises. Due to the fine structure brought in by the PF$^2$ process, the model naturally forms filament threads, aligned above the polarity inversion line. Synthetic views at our resolution of $15 \mathrm{km}$ show many details that may be verified in future high-resolution observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Super-Alfvénic Rotational Instability in accretion disks about black holes

The theory of instability of accretion disks about black holes, neutron stars or proto-planets, is revisited by means of the recent method of the Spectral Web. The cylindrical accretion disk differential equation is shown to be governed by the forward and backward Doppler-shifted continuous Alfvén spectra $Ω_{\rm A}^\pm \equiv m Ω\pm ω_{\rm A}$, where $ω_{\rm A}$ is the static Alfvén frequency. It is crucial to take non-axisymmetry ($m \ne 0$) and super-Alfvénic rotation of the Doppler frames ($|mΩ| \gg |ω_{\rm A}|$) into account. The continua $Ω_{\rm A}^+$ and $Ω_{\rm A}^-$ then overlap, ejecting a plethora of Super-Alfvénic Rotational Instabilities (SARIs). In-depth analysis for small inhomogeneity shows that the two Alfvén singularities reduce the extent of the modes to sizes much smaller than the width of the accretion disk. Generalization for large inhomogeneity leads to the completely unprecedented result that, for mode numbers $|k| \gg |m|$, any complex $ω$ in a wide neighborhood of the real axis is an approximate `eigenvalue'. The difference with genuine eigenmodes is that the amount of complementary energy to excite the modes is tiny, $|W_{\rm com}| \le c$, with $c$ the machine accuracy of the computation. This yields a multitude of two-dimensional continua of quasi-discrete modes: quasi-continuum SARIs. We conjecture that the onset of 3D turbulence in magnetized accretion disks is governed, not by the excitation of discrete axisymmetric Magneto-Rotational Instabilities, but by the excitation of modes from these two-dimensional continua of quasi-discrete non-axisymmetric Super-Alfvénic Rotational Instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Three-dimensional MHD wave propagation near a coronal null point: a new wave mode decomposition approach

We present a new MHD wave decomposition method that overcomes the limitations of existing wave identification methods. Our method allows to investigate the energy fluxes in different MHD modes at different locations of the solar atmosphere as waves generated by vortex flows travel through the solar atmosphere and pass near the magnetic null. We simulate wave dynamics through a coronal null configuration and apply a rotational wave driver at our bottom photospheric boundary. To identify the wave energy fluxes associated with different MHD wave modes, we employ a wave-decomposition method that is able to uniquely distinguish different MHD modes. Our proposed method utilizes the geometry of an individual magnetic field-line in 3D space to separate out velocity perturbations associated with the three fundamental MHD waves. Our method for wave identification is consistent with previous flux-surface-based methods and gives expected results in terms of wave energy fluxes at various locations of the null configuration. We show that ubiquitous vortex flows excite MHD waves that contribute significantly to the Poynting flux in the solar corona. Alfvén wave energy flux accumulates on the fan surface and fast wave energy flux accumulates near the null point. There is a strong current density buildup at the spine and fan surface.The proposed method has advantages over previously utilized wave decomposition methods, since it may be employed in realistic simulations or magnetic extrapolations, as well as in real solar observations, whenever the 3D fieldline shape is known. The enhancement in energy flux associated with magneto-acoustic waves near nulls may have important implications in the formation of jets and impulsive plasma flows.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Coronal rain in randomly heated arcades

Adopting the MPI-AMRVAC code, we present a 2.5-dimensional magnetohydrodynamic (MHD) simulation, which includes thermal conduction and radiative cooling, to investigate the formation and evolution of the coronal rain phenomenon. We perform the simulation in initially linear force-free magnetic fields which host chromospheric, transition region, and coronal plasma, with turbulent heating localized on their footpoints. Due to thermal instability, condensations start to occur at the loop top, and rebound shocks are generated by the siphon inflows. Condensations fragment into smaller blobs moving downwards and as they hit the lower atmosphere, concurrent upflows are triggered. Larger clumps show us clear "coronal rain showers" as dark structures in synthetic EUV hot channels and bright blobs with cool cores in the 304 Å channel, well resembling real observations. Following coronal rain dynamics for more than 10 hours, we carry out a statistical study of all coronal rain blobs to quantify their widths, lengths, areas, velocity distributions, and other properties. The coronal rain shows us continuous heating-condensation cycles, as well as cycles in EUV emissions. Compared to previous studies adopting steady heating, the rain happens faster and in more erratic cycles. Although most blobs are falling downward, upward-moving blobs exist at basically every moment. We also track the movement of individual blobs to study their dynamics and the forces driving their movements. The blobs have a prominence-corona transition-region-like structure surrounding them, and their movements are dominated by the pressure evolution in the very dynamic loop system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A fully self-consistent model for solar flares

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

MPI-AMRVAC: a parallel, grid-adaptive PDE toolkit

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Prominence formation by levitation-condensation at extreme resolutions

Prominences in the solar atmosphere represent an intriguing and delicate balance of forces and thermodynamics in an evolving magnetic topology. How this relatively cool material comes to reside at coronal heights, and what drives its evolution prior to, during, and after its appearance remains an area full of open questions. We deliberately focus on the levitation-condensation scenario, where a coronal flux rope forms and eventually demonstrates in-situ condensations, revisiting it at extreme resolutions down to order 6 km in scale. We perform grid-adaptive numerical simulations in a 2.5D translationally invariant setup, where we can study the distribution of all metrics involved in advanced magnetohydrodynamic stability theory for nested flux rope equilibria. We quantify in particular Convective Continuum Instability (CCI), Thermal Instability (TI), baroclinicity, and mass-slipping metrics within a series of numerical simulations of prominences formed via levitation-condensation. Overall, we find that the formation and evolution of prominence condensations happens in a clearly defined sequence regardless of resolution or background field strength between 3 and 10 Gauss. The CCI governs the slow evolution of plasma prior to the formation of distinct condensations found to be driven by the TI. Evolution of the condensations towards the topological dips of the magnetic flux rope is a consequence of these condensations forming initially outside of pressure balance with their surroundings. From the baroclinicity distributions, smaller-scale rotational motions are inferred within forming and evolving condensations. Upon the complete condensation of a prominence `monolith', the slow descent of this plasma towards lower heights appears consistent with the mass-slippage mechanism driven by episodes of both local current diffusion and magnetic reconnection.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermal instabilities: Fragmentation and field misalignment of filament fine structure

Prominences show a surprising amount of fine structure and it is widely believed that their threads, as seen in H$α$ 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$α$ observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation and plasma circulation of solar prominences

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pinwheels in the sky, with dust: 3D modeling of the Wolf-Rayet 98a environment

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Synthetic radio views on simulated solar flux ropes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Connecting the dots II: Phase changes in the climate dynamics of tidally locked terrestrial exoplanets

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solar prominences: 'double, double ... boil and bubble'

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Connecting the dots: A versatile model for the atmospheres of tidally locked Super-Earths

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effect of dust on Kelvin-Helmholtz instabilities

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interacting tilt and kink instabilities in repelling current channels

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rayleigh-Taylor instability in Magnetohydrodynamic Simulations of the Crab Nebula

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relativistic AGN jets II. Jet properties and mixing effects for episodic jet activity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating the in situ condensation process of solar prominences

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Dynamics of Funnel Prominences

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

3D Prominence-hosting Magnetic Configurations: Creating a Helical Magnetic Flux Rope

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relativistic 3D jet simulations for the X-ray binary SS433

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Three-Dimensional Magnetohydrodynamic Simulations of the Crab Nebula

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 σ=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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effect of angular opening on the dynamics of relativistic hydro jets

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multi-dimensional models of circumstellar shells around evolved massive stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Observation and Simulation of Longitudinal Oscillations of an Active Region Prominence

Filament longitudinal oscillations have been observed on the solar disk in H$α$. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations of Prominence Formation in the Magnetized Solar Corona by Chromospheric Heating

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

3-D simulations of shells around massive stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Computing the dust distribution in the bowshock of a fast moving, evolved star

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dust distribution in circumstellar shells

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the circumstellar medium of massive stars and how it may appear in GRB observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radiative cooling in numerical astrophysics: the need for adaptive mesh refinement

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thin shell morphology in the circumstellar medium of massive binaries

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.

0 citations0 reviewsNo code linkedOpen access