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

R. Keppens

R. Keppens contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRastro-ph.GAastro-ph.HEastro-ph.IMphysics.comp-phphysics.flu-dynphysics.plasm-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
7topics
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

7 published item(s)

preprint2022arXiv

BxC: a swift generator for 3D magnetohydrodynamic turbulence

Magnetohydrodynamic turbulence is central to laboratory and astrophysical plasmas, and is invoked for interpreting many observed scalings. Verifying predicted scaling law behaviour requires extreme-resolution direct numerical simulations (DNS), with needed computing resources excluding systematic parameter surveys. We here present an analytic generator of realistically looking turbulent magnetic fields, that computes 3D ${\cal{O}}(1000^3)$ solenoidal vector fields in minutes to hours on desktops. Our model is inspired by recent developments in 3D incompressible fluid turbulence theory, where a Gaussian white noise vector subjected to a non-linear transformation results in an intermittent, multifractal random field. Our $B\times C$ model has only few parameters that have clear geometric interpretations. We directly compare a (costly) DNS with a swiftly $B\times C$-generated realization, in terms of its (i) characteristic sheet-like structures of current density, (ii) volume-filling aspects across current intensity, (iii) power-spectral behaviour, (iv) probability distribution functions of increments for magnetic field and current density, structure functions, spectra of exponents, and (v) partial variance of increments. The model even allows to mimic time-evolving magnetic and current density distributions and can be used for synthetic observations on 3D turbulent data cubes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Two-fluid implementation in MPI-AMRVAC, with applications in the solar chromosphere

The chromosphere is a partially ionized layer of the solar atmosphere, the transition between the photosphere where the gas is almost neutral and the fully ionized corona. As the collisional coupling between neutral and charged particles decreases in the upper part of the chromosphere, the hydrodynamical timescales may become comparable to the collisional timescale, and a two-fluid model is needed. In this paper we describe the implementation and validation of a two-fluid model which simultaneously evolves charges and neutrals, coupled by collisions. The two-fluid equations are implemented in the fully open-source MPI-AMRVAC code. In the photosphere and the lower part of the solar atmosphere, where collisions between charged and neutral particles are very frequent, an explicit time-marching would be too restrictive, since for stability the timestep needs to be proportional to the inverse of the collision frequency. This is overcome by evaluating the collisional terms implicitly using an explicit-implicit (IMEX) scheme. The cases presented cover very different collisional regimes and our results are fully consistent with related literature findings. If collisional time and length scales are smaller than the hydrodynamical scales usually considered in the solar chromosphere, density structures seen in the neutral and charged fluids are similar, with the effect of elastic collisions between charges and neutrals being similar to diffusivity. Otherwise, density structures are different and the decoupling in velocity between the two species increases. The use of IMEX schemes efficiently avoids the small timestep constraints of fully explicit implementations in strongly collisional regimes. Adaptive Mesh Refinement (AMR) greatly decreases the computational cost, compared to uniform grid runs at the same effective resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Multi-threaded prominence oscillations triggered by a coronal shock wave

In this work, we study the causal relations between a localised energy release and a remote prominence oscillation, where the prominence has a realistic thread-like structure. We used an open source magnetohydrodynamic (MHD) code known as MPI-AMRVAC to create a multithreaded prominence body. We introduced an additional energy source from which a shock wave originates, thereby inducing prominence oscillation. We studied two cases with different source amplitudes to analyze its effect on the oscillations. Our results show that the frequently used pendulum model does not suffice to fully estimate the period of the prominence oscillation, in addition to showing that the influence of the source and the thread-like prominence structure needs to be taken into account. Repeated reflections and transmissions of the initial shock wave occur at the specific locations of multiple high-temperature and high-density gradients in the domain. This includes the left and right transition region (TR) located at the footpoints of the magnetic arcade, as well as the various transition regions between the prominence and the corona (PCTR). This results in numerous interferences of compressional waves. They contribute to the restoring forces of the oscillation, causing the period to deviate from the expected pendulum model, in addition to leading to differences in attributed damping or even growth in amplitude between the various threads. Along with the global longitudinal motion that result from the shock impact, small-scale transverse oscillations are also evident. Multiple high-frequency oscillations represent the propagation of magnetoacoustic waves. The damping we see is linked to the conversion of energy and its exchange with the surrounding corona. Our simulations demonstrate the exchange of energy between different threads and their different modes of oscillation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Radiation-Hydrodynamics with MPI-AMRVAC: Flux-Limited Diffusion

Radiation controls the dynamics and energetics of many astrophysical environments. To capture the coupling between the radiation and matter, however, is often a physically complex and computationally expensive endeavour. We develop a numerical tool to perform radiation-hydrodynamics simulations in various configurations at an affordable cost. We build upon the finite volume code MPI-AMRVAC to solve the equations of hydrodynamics on multi-dimensional adaptive meshes and introduce a new module to handle the coupling with radiation. A non-equilibrium, flux-limiting diffusion approximation is used to close the radiation momentum and energy equations. The time-dependent radiation energy equation is then solved within a flexible framework, accounting fully for radiation forces and work terms and further allowing the user to adopt a variety of descriptions for the radiation-matter interaction terms (the 'opacities'). We validate the radiation module on a set of standard testcases for which different terms of the radiative energy equation predominate. As a preliminary application to a scientific case, we calculate spherically symmetric models of the radiation-driven and optically thick supersonic outflows from massive Wolf-Rayet stars. This also demonstrates our code's flexibility, as the illustrated simulation combines opacities typically used in static stellar structure models with a parametrised form for the enhanced line-opacity expected in supersonic flows. This new module provides a convenient and versatile tool to perform multi-dimensional and high resolution radiative-hydrodynamics simulations in optically thick environments with the MPI-AMRVAC code. The code is ready to be used for a variety of astrophysical applications, where a first target for us will be multi-dimensional simulations of stellar outflows from Wolf-Rayet stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetohydrodynamic Non-linearities in Sunspot Atmospheres: Chromospheric Detections of Intermediate Shocks

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én shock solutions to the theoretical Rankine-Hugoniot relations. Employing over 200,000 independent measurements, we reveal that the Alfvén (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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Relativistic AGN jets III. Synthesis of synchrotron emission from Double-Double Radio Galaxies

The class of Double-Double Radio Galaxies (DDRGs) relates to episodic jet outbursts. How various regions and components add to the total intensity in radio images is less well known. In this paper we synthesize synchrotron images for DDRGs based on special relativistic hydrodynamic simulations, making advanced approximations for the magnetic fields. We study the synchrotron images for: Three different radial jet profiles; Ordered, entangled or mixed magnetic fields; Spectral ageing from synchrotron cooling; The contribution from different jet components; The viewing angle and Doppler (de-)boosting; The various epochs of the evolution of the DDRG. To link our results to observational data, we adopt to J1835+6204 as a reference source. In all cases the synthesized synchrotron images show two clear pairs of hotspots, in the inner and outer lobes. The best resemblance is obtained for the piecewise isochoric jet model, for a viewing angle of approximately $\vartheta \sim -71^{\circ}$, i.e. inclined with the lower jet towards the observer, with predominantly entangled ($\gtrsim 70$ per cent of the magnetic pressure) in turbulent, rather than ordered fields. The effects of spectral ageing become significant when the ratio of observation frequencies and cut-off frequency $ν_{\rm obs}/ν_{\infty,0} \gtrsim 10^{-3}$, corresponding to $\sim 3 \cdot 10^2$ MHz. For viewing angles $\vartheta \lesssim -30^{\circ}$, a DDRG morphology can no longer be recognized. The second jets must be injected within $\lesssim$ 4 per cent of the lifetime of the first jets for a DDRG structure to emerge, which is relevant for Active Galactic Nuclei feedback constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Wind morphology around cool evolved stars in binaries: the case of slowly accelerating oxygen-rich outflows

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.

0 citations0 reviewsNo code linkedOpen access