Refine
Year of publication
Document Type
- Article (13)
- Doctoral Thesis (3)
- Postprint (2)
- Habilitation Thesis (1)
Language
- English (19)
Is part of the Bibliography
- yes (19)
Keywords
- hydrodynamics (19) (remove)
Institute
We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh AREPO code. Therein, the anisotropic transport of cosmic rays (CRs) along magnetic field lines is discretized using a path-conservative finite volume method on the unstructured time-dependent Voronoi mesh of AREPO. The interaction of CRs and gyroresonant Alfven waves is described by short time-scale source terms in the CRHD equations. We employ a custom-made semi-implicit adaptive time stepping source term integrator to accurately integrate this interaction on the small light-crossing time of the anisotropic transport step. Both the transport and the source term integration step are separated from the evolution of the magnetohydrodynamical equations using an operator split approach. The new algorithm is tested with a variety of test problems, including shock tubes, a perpendicular magnetized discontinuity, the hydrodynamic response to a CR overpressure, CR acceleration of a warm cloud, and a CR blast wave, which demonstrate that the coupling between CR and magnetohydrodynamics is robust and accurate. We demonstrate the numerical convergence of the presented scheme using new linear and non-linear analytic solutions.
Photospheric radiation can drive winds from hot, massive stars by direct momentum transfer through scattering in bound-bound transitions of atmospheric ions. The line radiation force should cause a new radiative wave mode. The dispersion relation from perturbations of the line force was analysed so far either in Sobolev approximation or for pure line absorption. The former does not include the line-driven instability, and the latter cannot account for upstream propagating, radiative waves. We consider a non-Sobolev line force that includes scattering in a simplified way, accounting however for the important line-drag effect. We derive a new dispersion relation for radiative waves, and analyse wave propagation using Fourier methods, and by numerical solution of an integro-differential equation. The existence of an upstream propagating, dispersive radiative wave mode is demonstrated.
Cosmic rays (CRs) are a ubiquitous and an important component of astrophysical environments such as the interstellar medium (ISM) and intracluster medium (ICM). Their plasma physical interactions with electromagnetic fields strongly influence their transport properties. Effective models which incorporate the microphysics of CR transport are needed to study the effects of CRs on their surrounding macrophysical media. Developing such models is challenging because of the conceptional, length-scale, and time-scale separation between the microscales of plasma physics and the macroscales of the environment. Hydrodynamical theories of CR transport achieve this by capturing the evolution of CR population in terms of statistical moments. In the well-established one-moment hydrodynamical model for CR transport, the dynamics of the entire CR population are described by a single statistical quantity such as the commonly used CR energy density. In this work, I develop a new hydrodynamical two-moment theory for CR transport that expands the well-established hydrodynamical model by including the CR energy flux as a second independent hydrodynamical quantity. I detail how this model accounts for the interaction between CRs and gyroresonant Alfvén waves. The small-scale magnetic fields associated with these Alfvén waves scatter CRs which fundamentally alters CR transport along large-scale magnetic field lines. This leads to the effects of CR streaming and diffusion which are both captured within the presented hydrodynamical theory. I use an Eddington-like approximation to close the hydrodynamical equations and investigate the accuracy of this closure-relation by comparing it to high-order approximations of CR transport. In addition, I develop a finite-volume scheme for the new hydrodynamical model and adapt it to the moving-mesh code Arepo. This scheme is applied using a simulation of a CR-driven galactic wind. I investigate how CRs launch the wind and perform a statistical analysis of CR transport properties inside the simulated circumgalactic medium (CGM). I show that the new hydrodynamical model can be used to explain the morphological appearance of a particular type of radio filamentary structures found inside the central molecular zone (CMZ). I argue that these harp-like features are synchrotron-radiating CRs which are injected into braided magnetic field lines by a point-like source such as a stellar wind of a massive star or a pulsar. Lastly, I present the finite-volume code Blinc that uses adaptive mesh refinement (AMR) techniques to perform simulations of radiation and magnetohydrodynamics (MHD). The mesh of Blinc is block-structured and represented in computer memory using a graph-based approach. I describe the implementation of the mesh graph and how a diffusion process is employed to achieve load balancing in parallel computing environments. Various test problems are used to verify the accuracy and robustness of the employed numerical algorithms.
Cosmic-ray hydrodynamics
(2019)
Star formation in galaxies appears to be self-regulated by energetic feedback processes. Among the most promising agents of feedback are cosmic rays (CRs), the relativistic ion population of interstellar and intergalactic plasmas. In these environments, energetic CRs are virtually collisionless and interact via collective phenomena mediated by kinetic-scale plasma waves and large-scale magnetic fields. The enormous separation of kinetic and global astrophysical scales requires a hydrodynamic description. Here, we develop a new macroscopic theory for CR transport in the self-confinement picture, which includes CR diffusion and streaming. The interaction between CRs and electromagnetic fields of Alfvenic turbulence provides the main source of CR scattering, and causes CRs to stream along the magnetic field with the Alfven velocity if resonant waves are sufficiently energetic. However, numerical simulations struggle to capture this effect with current transport formalisms and adopt regularization schemes to ensure numerical stability. We extent the theory by deriving an equation for the CRmomentum density along the mean magnetic field and include a transport equation for the Alfven-wave energy. We account for energy exchange of CRs and Alfven waves via the gyroresonant instability and include other wave damping mechanisms. Using numerical simulations, we demonstrate that our new theory enables stable, self-regulated CR transport. The theory is coupled to magnetohydrodynamics, conserves the total energy and momentum, and correctly recovers previous macroscopic CR transport formalisms in the steady-state flux limit. Because it is free of tunable parameters, it holds the promise to provide predictable simulations of CR feedback in galaxy formation.
The space missions Voyager and Cassini together with earthbound observations re-vealed a wealth of structures in Saturn’s rings. There are, for example, waves being excited at ring positions which are in orbital resonance with Saturn’s moons. Other structures can be assigned to embedded moons like empty gaps, moon induced wakes or S-shaped propeller features. Further-more, irregular radial structures are observed in the range from 10 meters until kilometers. Here some of these structures will be discussed in the frame of hydrodynamical modeling of Saturn’s dense rings. For this purpose we will characterize the physical properties of the ring particle ensemble by mean field quantities and point to the special behavior of the transport coefficients. We show that unperturbed rings can become unstable and how diffusion acts in the rings. Additionally, the alternative streamline formalism is introduced to describe perturbed regions of dense rings with applications to the wake damping and the dispersion relation of the density waves.
One of the biggest successes of the Cassini mission is the detection of small moons (moonlets) embedded in Saturns rings that cause S-shaped density structures in their close vicinity, called propellers. Here, we present isothermal hydrodynamic simulations of moonlet-induced propellers in Saturn's A ring that denote a further development of the original model. We find excellent agreement between these new hydrodynamic and corresponding N-body simulations. Furthermore, the hydrodynamic simulations confirm the predicted scaling laws and the analytical solution for the density in the propeller gaps. Finally, this mean field approach allows us to simulate the pattern of the giant propeller Blériot, which is too large to be modeled by direct N-body simulations. Our results are compared to two stellar occultation observations by the Cassini Ultraviolet Imaging Spectrometer (UVIS), which intersect the propeller Blériot. Best fits to the UVIS optical depth profiles are achieved for a Hill radius of 590 m, which implies a moonlet diameter of about 860 m. Furthermore, the model favors a kinematic shear viscosity of the surrounding ring material of ν0 = 340 cm2 s−1, a dispersion velocity in the range of 0.3 cm s−1 < c0 < 1.5 cm s−1, and a fairly high bulk viscosity 7 < ξ0/ν0 < 17. These large transport values might be overestimated by our isothermal ring model and should be reviewed by an extended model including thermal fluctuations.
The space missions Voyager and Cassini together with earthbound observations revealed a wealth of structures in Saturn's rings. There are, for example, waves being excited at ring positions which are in orbital resonance with Saturn's moons. Other structures can be assigned to embedded moons like empty gaps, moon induced wakes or S-shaped propeller features. Furthermore, irregular radial structures are observed in the range from 10 meters until kilometers. Here some of these structures will be discussed in the frame of hydrodynamical modeling of Saturn's dense rings. For this purpose we will characterize the physical properties of the ring particle ensemble by mean field quantities and point to the special behavior of the transport coefficients. We show that unperturbed rings can become unstable and how diffusion acts in the rings. Additionally, the alternative streamline formalism is introduced to describe perturbed regions of dense rings with applications to the wake damping and the dispersion relation of the density waves.
X-ray observations of young Planetary Nebulæ (PNe) have revealed diffuse emission in extended regions around both H-rich and H-deficient central stars. In order to also repro-duce physical properties of H-deficient objects, we have, at first, extended our time-dependent radiation-hydrodynamic models with heat conduction for such conditions. Here we present some of the important physical concepts, which determine how and when a hot wind-blown bubble forms. In this study we have had to consider the, largely unknown, evolution of the CSPN, the slow (AGB) wind, the fast hot-CSPN wind, and the chemical composition. The main conclusion of our work is that heat conduction is needed to explain X-ray properties of wind-blown bubbles also in H-deficient objects.
Heat conduction has been found a plausible solution to explain discrepancies between expected and measured temperatures in hot bubbles of planetary nebulae (PNe). While the heat conduction process depends on the chemical composition, to date it has been exclusively studied for pure hydrogen plasmas in PNe. A smaller population of PNe show hydrogen-deficient and helium-and carbon-enriched surfaces surrounded by bubbles of the same composition; considerable differences are expected in physical properties of these objects in comparison to the pure hydrogen case. The aim of this study is to explore how a chemistry-dependent formulation of the heat conduction affects physical properties and how it affects the X-ray emission from PN bubbles of hydrogen-deficient stars. We extend the description of heat conduction in our radiation hydrodynamics code to work with any chemical composition. We then compare the bubble-formation process with a representative PN model using both the new and the old descriptions. We also compare differences in the resulting X-ray temperature and luminosity observables of the two descriptions. The improved equations show that the heat conduction in our representative model of a hydrogen-deficient PN is nearly as efficient with the chemistry-dependent description; a lower value on the diffusion coefficient is compensated by a slightly steeper temperature gradient. The bubble becomes somewhat hotter with the improved equations, but differences are otherwise minute. The observable properties of the bubble in terms of the X-ray temperature and luminosity are seemingly unaffected.
X-ray observations of the double-binary OB-star system QZ Car (HD 93206) obtained with the Chandra X-ray Observatory over a period of roughly 2 years are presented. The respective orbits of systems A (O9.7 I+b2 v, P-A = 21 days) and B (O8 III+o9 v, P-B = 6 days) are reasonably well sampled by the observations, allowing the origin of the X-ray emission to be examined in detail. The X-ray spectra can be well fitted by an attenuated three-temperature thermal plasma model, characterized by cool, moderate, and hot plasma components at kT similar or equal to 0.2, 0.7, and 2 keV, respectively, and a circumstellar absorption of similar or equal to 0.2 x 10(22) cm(-2). Although the hot plasma component could be indicating the presence of wind-wind collision shocks in the system, the model fluxes calculated from spectral fits, with an average value of similar or equal to 7x10(-13) erg s(-1) cm(-2), do not show a clear correlation with the orbits of the two constituent binaries. A semi-analytical model of QZ Car reveals that a stable momentum balance may not be established in either system A or B. Yet, despite this, system B is expected to produce an observed X-ray flux well in excess of the observations. If one considers the wind of the O8 III star to be disrupted by mass transfer, the model and observations are in far better agreement, which lends support to the previous suggestion of mass transfer in the O8 III+o9 v binary. We conclude that the X-ray emission from QZ Car can be reasonably well accounted for by a combination of contributions mainly from the single stars and the mutual wind-wind collision between systems A and B.