Researcher profile

Allard Jan van Marle

Allard Jan van Marle contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

9 published item(s)

preprint2022arXiv

Diffusive shock acceleration at oblique high Mach number shocks

The current paradigm of cosmic ray (CR) origin states that the most part of galactic CRs is produced by supernova remnants. The interaction of supernova ejecta with the interstellar medium after supernova's explosions results in shocks responsible for CR acceleration via diffusive shock acceleration (DSA). We use particle-in-cell (PIC) simulations and a combined PIC-magnetohydrodynamic (PIC-MHD) technique to investigate whether DSA can occur in oblique high Mach number shocks. Using the PIC method, we follow the formation of the shock and determine the fraction of the particles that gets involved in DSA. Then, with this result, we use PIC-MHD simulations to model the large-scale structure of the plasma and the magnetic field surrounding the shock and find out whether or not the reflected particles can generate the upstream turbulence and trigger DSA. We find that the feasibility of this process in oblique shocks depends strongly on the Alfvenic Mach number, and the DSA process is more likely triggered at high Mach number shocks.

preprint2020arXiv

On the influence of supra-thermal particle acceleration on the morphology of low-Mach, high-$β$ shocks

When two galaxy clusters encounter each other, the interaction results in a collisionless shock that is characterized by a low (1-4) sonic Mach number, and a high Alfv{é}nic Mach number. Our goal is to determine if, and to what extent, such shocks can accelerate particles to sufficient velocities that they can contribute to the cosmic ray spectrum. We combine two different computational methods, magnetohydrodynamics (MHD) and particle-in-cell (PIC) into a single code that allows us to take advantage of the high computational efficiency of MHD while maintaining the ability to model the behaviour of individual non-thermal particles. Using this method, we perform a series of simulations covering the expected parameter space of galaxy cluster collision shocks. Our results show that for shocks with a sonic Mach number below 2.25 no diffusive shock acceleration can take place because of a lack of instabilities in the magnetic field, whereas for shocks with a sonic Mach number $\geq\,3$ the acceleration is efficient and can accelerate particles to relativistic speeds. In the regime between these two extremes, diffusive shock acceleration can occur but is relatively inefficient because of the time- and space-dependent nature of the instabilities. For those shocks that show efficient acceleration, the instabilities in the upstream gas increase to the point where they change the nature of the shock, which, in turn, will influence the particle injection process.

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.

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.

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.

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.

preprint2010arXiv

Numerical models of collisions between core-collapse supernovae and circumstellar shells

Recent observations of luminous Type IIn supernovae (SNe) provide compelling evidence that massive circumstellar shells surround their progenitors. In this paper we investigate how the properties of such shells influence the SN lightcurve by conducting numerical simulations of the interaction between an expanding SN and a circumstellar shell ejected a few years prior to core collapse. Our parameter study explores how the emergent luminosity depends on a range of circumstellar shell masses, velocities, geometries, and wind mass-loss rates, as well as variations in the SN mass and energy. We find that the shell mass is the most important parameter, in the sense that higher shell masses (or higher ratios of M_shell/M_SN) lead to higher peak luminosities and higher efficiencies in converting shock energy into visual light. Lower mass shells can also cause high peak luminosities if the shell is slow or if the SN ejecta are very fast, but only for a short time. Sustaining a high luminosity for durations of more than 100 days requires massive circumstellar shells of order 10 M_sun or more. This reaffirms previous comparisons between pre-SN shells and shells produced by giant eruptions of luminous blue variables (LBVs), although the physical mechanism responsible for these outbursts remains uncertain. The lightcurve shape and observed shell velocity can help diagnose the approximate size and density of the circumstellar shell, and it may be possible to distinguish between spherical and bipolar shells with multi-wavelength lightcurves. These models are merely illustrative. One can, of course, achieve even higher luminosities and longer duration light curves from interaction by increasing the explosion energy and shell mass beyond values adopted here.

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.

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.