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.

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

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

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.