BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Rony Keppens

Rony Keppens contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRphysics.plasm-phastro-ph.HEastro-ph.IMphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 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

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