Search

Oceans play a big role in the nature of our planet, about $ 70 \% $ of our earth is covered by water. Strong currents are transporting warm water around the world making life possible, and allowing us to harvest its power producing energy. Yet, oceans also carry a much more deadly side. Floods and tsunamis can easily annihilate whole cities and destroy life in seconds. The earth's climate system is also very much linked to the currents in the ocean due to its large coverage of the earth's surface, thus, gaining scientific insights into the mechanisms and effects through simulations is of high importance. Deep ocean currents can be simulated by means of wall-bounded turbulent flow simulations. To support these very large scale numerical simulations and enable the scientists to interpret their output, we deploy an interactive visualization framework to study sheared thermal convection. The visualizations are based on volume rendering of the temperature field. To address the needs of supercomputer users with different hardware and software resources, we evaluate different volume rendering implementations supported in the ParaView environment: two GPU-based solutions with Kitware's native volume mapper or NVIDIA's IndeX library, and a CPU-only Intel OSPRay-based implementation.

Collisionless magnetic reconnection is a prime candidate to account for flare-like or steady emission, outflow launching, or plasma heating, in a variety of high-energy astrophysical objects, including ones with relativistic ion-electron plasmas. But the fate of the initial magnetic energy in a reconnection event remains poorly known: what is the amount given to kinetic energy, the ion/electron repartition, and the hardness of the particle distributions? We explore these questions with 2D particle-in-cell simulations of ion-electron plasmas. We find that 45 to 75% of the total initial magnetic energy ends up in kinetic energy, this fraction increasing with the inflow magnetization. Depending on the guide field strength, ions get from 30 to 60% of the total kinetic energy. Particles can be separated into two populations that only weakly mix: (i) particles initially in the current sheet, heated by its initial tearing and subsequent contraction of the islands; and (ii) particles from the background plasma that primarily gain energy via the reconnection electric field when passing near the X-point. Particles (ii) tend to form a power-law with an index $p=-d\log n(\gamma)/d\log\gamma$, that depends mostly on the inflow Alfv\'en speed $V_A$ and magnetization $\sigma_s$ of species $s$, with for electrons $p=5$ to $1.2$ for increasing $\sigma_e$. The highest particle Lorentz factor, for ions or electrons, increases roughly linearly with time for all the relativistic simulations. This is faster, and the spectra can be harder, than for collisionless shock acceleration. We discuss applications to microquasar and AGN coronae, to extragalactic jets, and to radio lobes. We point out situations where effects such as Compton drag or pair creation are important., Comment: 15 pages, submitted to A&A

We present first results of a long-term project which aims at multi-scale, multi-physics simulations of wind accretion in microquasars and high-mass X-ray binaries. The 3D hydrodynamical simulations cover all scales, from the circum-binary environment down to the immediate vicinity of the black hole. We first introduce the numerical method and parallelization strategy of the AMR A-MAZE code. We then discuss some preliminary results of how, and on what scales, an accretion disk is formed around the black hole. We finally present some characteristics of this disk, which is far from Keplerian. We emphasize that on all scales shocks play a decisive role for the accretion process and the process of structure formation -- for the formation of the large scale, nearly coherent structure of the disk, but also for the formation of turbulent fluctuations., Comment: Invited contribution to ASTRONUM 2013

Magnetic reconnection is a leading mechanism for magnetic energy conversion and high-energy non-thermal particle production in a variety of high-energy astrophysical objects, including ones with relativistic ion-electron plasmas (e.g., microquasars or AGNs) - a regime where first principle studies are scarce. We present 2D particle-in-cell (PIC) simulations of low $\beta$ ion-electron plasmas under relativistic conditions, i.e., with inflow magnetic energy exceeding the plasma rest-mass energy. We identify outstanding properties: (i) For relativistic inflow magnetizations (here $10 < \sigma_e < 360$), the reconnection outflows are dominated by thermal agitation instead of bulk kinetic energy. (ii) At large inflow electron magnetization ($\sigma_e > 80$), the reconnection electric field is sustained more by bulk inertia than by thermal inertia. It challenges the thermal-inertia-paradigm and its implications. (iii) The inflows feature sharp transitions at the entrance of the diffusion zones. These are not shocks but results from particle ballistic motions, all bouncing at the same location, provided that the thermal velocity in the inflow is far smaller than the inflow E cross B bulk velocity. (iv) Island centers are magnetically isolated from the rest of the flow, and can present a density depletion at their center. (v) The reconnection rates are slightly larger than in non-relativistic studies. They are best normalized by the inflow relativistic Alfv\'en speed projected in the outflow direction, which then leads to rates in a close range (0.14-0.25) thus allowing for an easy estimation of the reconnection electric field., Comment: Submitted to A&A

Colliding supersonic bulk flows shape observable properties and internal physics of various astrophysical objects, like O-star winds, molecular clouds, galactic sheets, binaries, or gamma-ray bursts. Using numerical simulations, we show that the bulk flows leave a clear imprint on the collision zone, its mean properties and the turbulence it naturally develops. Our model setup consists of 3D head-on colliding isothermal hydrodynamical flows with Mach numbers between 2 and 43. Simulation results are in line with expectations from self-similarity: root mean square Mach numbers (Mrms) scale linearly with upstream Mach numbers, mean densities remain limited to a few times the upstream density. The density PDF is not log-normal. The turbulence is inhomogeneous: weaker in the zone center than close to the confining shocks. It is anisotropic: while Mrms is generally supersonic, Mrms transverse to the upstream flow is always subsonic. We argue that uniform, isothermal, head-on colliding flows generally disfavor isotropic, supersonic turbulence. The anisotropy carries over to other quantities like the density variance - Mach number relation. Structure functions differ depending on whether they are computed along a line-of-sight perpendicular or parallel to the upstream flow. We suggest that such line-of-sight effects should be kept in mind when interpreting turbulence characteristics derived from observations., Comment: 20 pages, 14 figures, 4 tables, accepted by Astronomy and Astrophysics

Aims. Observations show nebulae around some massive stars but not around others. If observed, their chemical composition is far from homogeneous. Our goal is to put these observational features into the context of the evolution of massive stars and their circumstellar medium (CSM) and, more generally, to quantify the role of massive stars for the chemical and dynamical evolution of the ISM. Methods. Using the A-MAZE code, we perform 2d-axisymmetric hydrodynamical simulations of the evolution of the CSM, shaped by stellar winds, for a whole grid of massive stellar models from 15 to 120 Msun and following the stellar evolution from the zero-age main-sequence to the time of supernova explosion. In addition to the usual quantities, we also follow five chemical species: H, He, C, N, and O. Results. We show how various quantities evolve as a function of time: size of the bubble, position of the wind termination shock, chemical composition of the bubble, etc. The chemical composition of the bubble changes considerably compared to the initial composition, particularly during the red-supergiant (RSG) and Wolf-Rayet (WR) phases. In some extreme cases, the inner region of the bubble can be completely depleted in hydrogen and nitrogen, and is mainly composed of carbon, helium and oxygen. We argue why the bubble typically expands at a lower rate than predicted by self-similarity theory. In particular, the size of the bubble is very sensitive to the density of the ISM, decreasing by a factor of around 2.5 for each additional dex in ISM density. The bubble size also decreases with the metallicity of the central star, as low-metallicity stars have weaker winds. Our models qualitatively fit the observations of WR ejecta nebulae., Comment: 13 pages, 14 figures, accepted for publication in A&A. Corrected after language editing

We present the parallel particle-in-cell (PIC) code Apar-T and, more importantly, address the fundamental question of the relations between the PIC model, the Vlasov-Maxwell theory, and real plasmas. First, we present four validation tests: spectra from simulations of thermal plasmas, linear growth rates of the relativistic tearing instability and of the filamentation instability, and non-linear filamentation merging phase. For the filamentation instability we show that the effective growth rates measured on the total energy can differ by more than 50% from the linear cold predictions and from the fastest modes of the simulation. Second, we detail a new method for initial loading of Maxwell-J\"uttner particle distributions with relativistic bulk velocity and relativistic temperature, and explain why the traditional method with individual particle boosting fails. Third, we scrutinize the question of what description of physical plasmas is obtained by PIC models. These models rely on two building blocks: coarse-graining, i.e., grouping of the order of p~10^10 real particles into a single computer superparticle, and field storage on a grid with its subsequent finite superparticle size. We introduce the notion of coarse-graining dependent quantities, i.e., quantities depending on p. They derive from the PIC plasma parameter Lambda^{PIC}, which we show to scale as 1/p. We explore two implications. One is that PIC collision- and fluctuation-induced thermalization times are expected to scale with the number of superparticles per grid cell, and thus to be a factor p~10^10 smaller than in real plasmas. The other is that the level of electric field fluctuations scales as 1/Lambda^{PIC} ~ p. We provide a corresponding exact expression. Fourth, we compare the Vlasov-Maxwell theory, which describes a phase-space fluid with infinite Lambda, to the PIC model and its relatively small Lambda., Comment: 24 pages, 14 figures, accepted in Astronomy & Astrophysics

We present 3D hydrodynamical simulations of the separated binary RS Ophiuchi (RS Oph), a recurrent nova and potential progenitor of a SN Ia. RS Oph is composed of a red giant (RG) and a white dwarf (WD) whose mass is close to the Chandrasekhar limit. In an isothermal scenario, the WD accrets about 10% of a 20 km/s RG wind by a non-Keplerian accretion disk with strong spiral shocks, and about 2% of a 60 km/s RG wind by what we term a 'turbulent accretion ball'. A significantly larger impact have the thermodynamics. In an adiabatic scenario only about 0.7% of the 20 km/s RG wind is accreted. The rate of change of the system separation due to mass and angular momentum loss out of the system is negative in all three cases studied, but is ten times smaller for a fast RG wind (60 km/s) than for a slow RG wind (20 km/s). The results demonstrate that existing nova models and observed recurrence times fit well together with 3D wind accretion and that RS Oph is one of the most promising systems to become an SN Ia., Comment: 6 pages, 2 figures; ASTRONUM-2009 "Numerical Modeling of Space Plasma Flows", Chamonix, France, July 2009, ASP Conf. Proc., in press

Shock bound interaction zones (SBIZs) are ubiquitous in astrophysics. We present numerical results for 2D and 3D, plane-parallel, infinitely extended SBIZs. Isothermal settings and parameterized cooling are considered. We highlight and compare characteristic of such zones. We emphasize the mutual coupling between the turbulence within the SBIZ and the confining shocks, point out potential differences to 3D periodic box studies of supersonic turbulence, and contemplate on possible effects on the X-ray emission of such zones., Comment: 7 pages, 4 figures; ASTRONUM-2009 "Numerical Modeling of Space Plasma Flows", Chamonix, France, July 2009, ASP Conf. Proc., in press