TY - JOUR A1 - Thomas, Timon A1 - Pfrommer, Christoph A1 - Pakmor, Rüdiger T1 - A finite volume method for two-moment cosmic ray hydrodynamics on a moving mesh JF - Monthly notices of the Royal Astronomical Society N2 - 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. KW - hydrodynamics KW - MHD KW - methods: numerical KW - cosmic rays Y1 - 2021 U6 - https://doi.org/10.1093/mnras/stab397 SN - 0035-8711 SN - 1365-2966 VL - 503 IS - 2 SP - 2242 EP - 2264 PB - Oxford University Press CY - Oxford ER - TY - JOUR A1 - Krtička, Jiří A1 - Feldmeier, Achim T1 - Stochastic light variations in hot stars from wind instability BT - finding photometric signatures and testing against the TESS data JF - Astronomy and astrophysics : an international weekly journal / European Southern Observatory (ESO) N2 - Context Line-driven wind instability is expected to cause small-scale wind inhomogeneities, X-ray emission, and wind line profile variability. The instability can already develop around the sonic point if it is initiated close to the photosphere due to stochastic turbulent motions. In such cases, it may leave its imprint on the light curve as a result of wind blanketing. Aims We study the photometric signatures of the line-driven wind instability. Methods We used line-driven wind instability simulations to determine the wind variability close to the star. We applied two types of boundary perturbations: a sinusoidal one that enables us to study in detail the development of the instability and a stochastic one given by a Langevin process that provides a more realistic boundary perturbation. We estimated the photometric variability from the resulting mass-flux variations. The variability was simulated assuming that the wind consists of a large number of independent conical wind sectors. We compared the simulated light curves with TESS light curves of OB stars that show stochastic variability. Results We find two typical signatures of line-driven wind instability in photometric data: a knee in the power spectrum of magnitude fluctuations, which appears due to engulfment of small-scale structure by larger structures, and a negative skewness of the distribution of fluctuations, which is the result of spatial dominance of rarefied regions. These features endure even when combining the light curves from independent wind sectors. Conclusions The stochastic photometric variability of OB stars bears certain signatures of the line-driven wind instability. The distribution function of observed photometric data shows negative skewness and the power spectra of a fraction of light curves exhibit a knee. This can be explained as a result of the line-driven wind instability triggered by stochastic base perturbations. KW - stars: winds KW - outflows KW - stars: mass-loss KW - stars: early-type KW - hydrodynamics KW - instabilities KW - stars: variables: general Y1 - 2021 U6 - https://doi.org/10.1051/0004-6361/202040148 SN - 1432-0746 VL - 648 PB - EDP Sciences CY - Les Ulis ER - TY - THES A1 - Thomas, Timon T1 - Cosmic-ray hydrodynamics: theory, numerics, applications T1 - Hydrodynamik der kosmischen Strahlung: Theorie, Numerik, Anwendungen N2 - 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. N2 - Kosmische Strahlung (CR) ist ein allgegenwärtiger und wichtiger Bestandteil astrophysikalischer Umgebungen wie des interstellaren Mediums (ISM) und des Intracluster-Mediums (ICM). Ihre plasmaphysikalischen Wechselwirkungen mit elektromagnetischen Feldern beeinflussen ihre Transporteigenschaften weitgehend. Effektive Modelle, die die Mikrophysik des CR-Transports einbeziehen, sind erforderlich, um die Auswirkungen von CRs auf die sie umgebenden makrophysikalischen Medien zu untersuchen. Die Entwicklung solcher Modelle ist eine Herausforderung, aufgrund der konzeptionellen, Längenskalen-, und Zeitskalen-Unterschiede zwischen den Mikroskalen der Plasmaphysik und den Makroskalen der Umgebung. Hydrodynamische Theorien des CR-Transports erreichen dies, indem sie die Entwicklung der CR-Population in Form von statistischen Momenten erfassen. Im etablierten hydrodynamischen Ein-Moment Modell für den CR-Transport wird die Dynamik der gesamten CR-Population durch eine einzige statistische Größe wie der häufig verwendeten CR-Energiedichte beschrieben. In dieser Arbeit entwickle ich eine neue hydrodynamische Zwei-Momenten Theorie für den CR-Transport, die das etablierte hydrodynamische Modell um den CR-Energiefluss als zweite unabhängige hydrodynamische Größe erweitert. Ich erläutere, wie dieses Modell die Wechselwirkung zwischen CRs und gyroresonanten Alfvén-Wellen berücksichtigt. Die mit diesen Alfvén-Wellen verbundenen kleinskaligen Magnetfelder streuen die CRs, was den CR-Transport entlang großskaligen Magnetfeldlinien grundlegend verändert. Dies führt zu den CR-Strömungs-und Diffusioneffekten, welche beide in der neu vorgestellten hydrodynamischen Theorie erfasst werden. Ich verwende eine adaptierte Eddington Näherung, um die hydrodynamischen Gleichungen zu schließen und untersuche die Genauigkeit dieser Näherung, indem ich sie mit Näherungen höherer Ordnung für den CR-Transport vergleiche. Darüber hinaus entwickle ich ein Finite-Volumen-Schema für das neue hydrodynamische Modell und passe es an den mitbewegten Gitter Code Arepo an. Dieses Schema wird mittels einer Simulation eines CR-getriebenen galaktischen Windes angewendet. Ich untersuche, wie CRs den Wind beschleunigen und führe eine statistische Analyse der CR-Transporteigenschaften innerhalb des simulierten zirkumgalaktischen Mediums (CGM) durch. Ich zeige, dass das neue hydrodynamische Modell das morphologische Erscheinungsbild eines neu-entdeckten bestimmten Typs von filamentartigen Radiostrukturen, welcher in der zentralen molekularen Zone (CMZ) auffindbar ist, erklären kann. Ich schlage vor, dass es sich bei diesen harfenartigen Strukturen um synchrotronstrahlende CRs handelt, die zuvor von einer punktförmigen Quelle wie dem stellaren Wind eines massereichen Sterns oder eines Pulsars in geflochtene Magnetfeldlinien injiziert wurden. Schließlich stelle ich den Finite-Volumen-Code Blinc vor, der adaptive Gitterverfeinerungstechniken (AMR) verwendet, um Simulationen von Strahlungs-und Magnetohydrodynamik (MHD) durchzuführen. Das Gitter von Blinc ist blockstrukturiert und wird im Computerspe-icher mittels eines graphbasierten Ansatzes dargestellt. Ich beschreibe die Implementierung des Gittergraphen und wie ein Diffusionsprozess eingesetzt wird, um einen Lastausgleich in parallelen Rechenumgebungen zu erreichen. Verschiedene Testprobleme werden verwendet, um die Genauigkeit und Robustheit der verwendeten numerischen Algorithmen zu überprüfen. KW - cosmic rays KW - hydrodynamics KW - radiative transfer KW - methods: analytical KW - methods: numerical KW - Galactic center KW - Non-thermal radiation sources KW - galaktisches Zentrum KW - Quellen nichtthermischer Strahlung KW - kosmische Strahlung KW - Hydrodynamik KW - Methoden: analytisch KW - Methoden: numerisch KW - Strahlungstransport Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-563843 ER - TY - JOUR A1 - Kurfürst, P. A1 - Feldmeier, Achim A1 - Krticka, Jiri T1 - Two-dimensional modeling of density and thermal structure of dense circumstellar outflowing disks JF - Astronomy and astrophysics : an international weekly journal N2 - Context. Evolution of massive stars is affected by a significant loss of mass either via (nearly) spherically symmetric stellar winds or by aspherical mass-loss mechanisms, namely the outflowing equatorial disks. However, the scenario that leads to the formation of a disk or rings of gas and dust around massive stars is still under debate. It is also unclear how various forming physical mechanisms of the circumstellar environment affect its shape and density, as well as its kinematic and thermal structure. Results. Our models show the geometric distribution and contribution of viscous heating that begins to dominate in the central part of the disk for mass-loss rates higher than (M) over dot greater than or similar to 10(-10) M-circle dot yr(-1). In the models of dense viscous disks with (M) over dot > 10(-8) M-circle dot yr(-1), the viscosity increases the central temperature up to several tens of thousands of Kelvins, however the temperature rapidly drops with radius and with distance from the disk midplane. The high mass-loss rates and high viscosity lead to instabilities with significant waves or bumps in density and temperature in the very inner disk region. Conclusions. The two-dimensional radial-vertical models of dense outflowing disks including the full Navier-Stokes viscosity terms show very high temperatures that are however limited to only the central disk cores inside the optically thick area, while near the edge of the optically thick region the temperature may be low enough for the existence of neutral hydrogen, for example. KW - stars: massive KW - stars: mass-loss KW - stars: winds, outflows KW - stars: evolution KW - stars: rotation KW - hydrodynamics Y1 - 2018 U6 - https://doi.org/10.1051/0004-6361/201731300 SN - 1432-0746 VL - 613 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - del Valle, Maria Victoria A1 - Pohl, Martin T1 - Nonthermal emission from Stellar Bow Shocks JF - The astrophysical journal : an international review of spectroscopy and astronomical physics N2 - Since the detection of nonthermal radio emission from the bow shock of the massive runaway star BD +43 degrees 3654, simple models have predicted high-energy emission, at X-rays and gamma-rays, from these Galactic sources. Observational searches for this emission so far give no conclusive evidence but a few candidates at gamma-rays. In this work we aim at developing a more sophisticated model for the nonthermal emission from massive runaway star bow shocks. The main goal is to establish whether these systems are efficient nonthermal emitters, even if they are not strong enough yet to be detected. For modeling the collision between the stellar wind and the interstellar medium we use 2D hydrodynamic simulations. We then adopt the flow profile of the wind and the ambient medium obtained with the simulation as the plasma state for solving the transport of energetic particles injected in the system, as well as the nonthermal emission they produce. For this purpose we solve a 3D (two spatial vertical bar energy) advection-diffusion equation in the test-particle approximation. We find that a massive runaway star with a powerful wind converts 0.16%-0.4% of the power injected in electrons into nonthermal emission, mostly produced by inverse Compton scattering of dust-emitted photons by relativistic electrons, and second by synchrotron radiation. This represents a fraction of similar to 10(-5) to 10(-4) of the wind kinetic power. Given the better sensibility of current instruments at radio wavelengths, these systems are more prone to be detected at radio through the synchrotron emission they produce rather than at gamma energies. KW - gamma-rays: stars KW - hydrodynamics KW - radiation mechanisms: nonthermal KW - stars: winds, outflows Y1 - 2018 U6 - https://doi.org/10.3847/1538-4357/aad333 SN - 0004-637X SN - 1538-4357 VL - 864 IS - 1 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Krticka, Jiri A1 - Feldmeier, Achim T1 - Light variations due to the line-driven wind instability and wind blanketing in O stars JF - Astronomy and astrophysics : an international weekly journal N2 - A small fraction of the radiative flux emitted by hot stars is absorbed by their winds and redistributed towards longer wavelengths. This effect, which leads also to the heating of the stellar photosphere, is termed wind blanketing. For stars with variable winds, the effect of wind blanketing may lead to the photometric variability. We have studied the consequences of line driven wind instability and wind blanketing for the light variability of O stars. We combined the results of wind hydrodynamic simulations and of global wind models to predict the light variability of hot stars due to the wind blanketing and instability. The wind instability causes stochastic light variability with amplitude of the order of tens of millimagnitudes and a typical timescale of the order of hours for spatially coherent wind structure. The amplitude is of the order of millimagnitudes when assuming that the wind consists of large number of independent concentric cones. The variability with such amplitude is observable using present space borne photometers. We show that the simulated light curve is similar to the light curves of O stars obtained using BRITE and CoRoT satellites. KW - stars: winds, outflows KW - stars: mass-loss KW - stars: early-type KW - stars: variables: general KW - hydrodynamics Y1 - 2018 U6 - https://doi.org/10.1051/0004-6361/201731614 SN - 1432-0746 VL - 617 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Seiß, Martin A1 - Albers, Nicole A1 - Sremčević, Miodrag A1 - Schmidt, Jürgen A1 - Salo, Heikki A1 - Seiler, Michael A1 - Hoffmann, Holger A1 - Spahn, Frank T1 - Hydrodynamic Simulations of Moonlet-induced Propellers in Saturn's Rings BT - Application to Bleriot JF - The astronomical journal N2 - 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. KW - diffusion KW - hydrodynamics KW - planets and satellites: rings Y1 - 2018 U6 - https://doi.org/10.3847/1538-3881/aaed44 SN - 0004-6256 SN - 1538-3881 VL - 157 IS - 1 PB - IOP Publishing Ltd. CY - Bristol ER - TY - JOUR A1 - Grätz, Fabio M. A1 - Seiß, Martin A1 - Schmidt, Jürgen A1 - Colwell, Joshua A1 - Spahn, Frank T1 - Sharp Gap Edges in Dense Planetary Rings BT - an Axisymmetric Diffusion Model JF - The astrophysical journal : an international review of spectroscopy and astronomical physics N2 - One of the most intriguing facets of Saturn's rings are the sharp edges of gaps in the rings where the surface density abruptly drops to zero. This is despite of the fact that the range over which a moon transfers angular momentum onto the ring material is much larger. Recent UVIS-scans of the edges of the Encke and Keeler gap show that this drop occurs over a range approximately equal to the rings' thickness. Borderies et al. show that this striking feature is likely related to the local reversal of the usually outward directed viscous transport of angular momentum in strongly perturbed regions. In this article we revise the Borderies et al. model using a granular flow model to define the shear and bulk viscosities, ν and ζ, and incorporate the angular momentum flux reversal effect into the axisymmetric diffusion model we developed for gaps in dense planetary rings. Finally, we apply our model to the Encke and Keeler division in order to estimate the shear and bulk viscosities in the vicinity of both gaps KW - celestial mechanics KW - diffusion KW - hydrodynamics KW - planets and satellites: rings KW - scattering Y1 - 2019 U6 - https://doi.org/10.3847/1538-4357/ab007e SN - 0004-637X SN - 1538-4357 VL - 872 IS - 2 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Thomas, T. A1 - Pfrommer, Christoph T1 - Cosmic-ray hydrodynamics BT - alfvén-wave regulated transport of cosmic rays JF - Monthly notices of the Royal Astronomical Society N2 - 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. KW - hydrodynamics KW - radiative transfer KW - methods: analytical KW - methods: numerical KW - cosmic rays Y1 - 2019 U6 - https://doi.org/10.1093/mnras/stz263 SN - 0035-8711 SN - 1365-2966 VL - 485 IS - 3 SP - 2977 EP - 3008 PB - Oxford Univ. Press CY - Oxford ER - TY - THES A1 - Grätz, Fabio M. T1 - Nonlinear diffusion in granular gases and dense planetary rings N2 - Small moonlets or moons embedded in dense planetary rings create S-shaped density modulations called propellers if their masses are smaller than a certain threshold, alternatively they create a circumferential gap in the disk if the embedded body’s mass exceeds this threshold (Spahn and Sremčević, 2000). The gravitational perturber scatters the ring particles, depletes the disk’s density, and, thus, clears a gap, whereas counteracting viscous diffusion of the ring material has the tendency to close the created gap, thereby forming a propeller. Propeller objects were predicted by Spahn and Sremčević (2000) and Sremčević et al. (2002) and were later discovered by the Cassini space probe (Tiscareno et al., 2006, Sremčević et al., 2007, Tiscareno et al., 2008, and Tiscareno et al., 2010). The ring moons Pan and Daphnis are massive enough to maintain the circumferential Encke and Keeler gaps in Saturn’s A ring and were detected by Showalter (1991) and Porco (2005) in Voyager and Cassini images, respectively. In this thesis, a nonlinear axisymmetric diffusion model is developed to describe radial density profiles of circumferential gaps in planetary rings created by embedded moons (Grätz et al., 2018). The model accounts for the gravitational scattering of the ring particles by the embedded moon and for the counteracting viscous diffusion of the ring matter back into the gap. With test particle simulations it is shown that the scattering of the ring particles passing the moon is larger for small impact parameters than estimated by Goldreich and Tremaine (1980). This is especially significant for the modeling of the Keeler gap. The model is applied to the Encke and Keeler gaps with the aim to estimate the shear viscosity of the ring in their vicinities. In addition, the model is used to analyze whether tiny icy moons whose dimensions lie below Cassini’s resolution capabilities would be able to cause the poorly understood gap structure of the C ring and the Cassini Division. One of the most intriguing facets of Saturn’s rings are the extremely sharp edges of the Encke and Keeler gaps: UVIS-scans of their gap edges show that the optical depth drops from order unity to zero over a range of far less than 100 m, a spatial scale comparable to the ring’s vertical extent. This occurs despite the fact that the range over which a moon transfers angular momentum onto the ring material is much larger. Borderies et al. (1982, 1989) have shown that this striking feature is likely related to the local reversal of the usually outward-directed viscous transport of angular momentum in strongly perturbed regions. We have revised the Borderies et al. (1989) model using a granular flow model to define the shear and bulk viscosities, ν and ζ, in order to incorporate the angular momentum flux reversal effect into the axisymmetric diffusion model for circumferential gaps presented in this thesis (Grätz et al., 2019). The sharp Encke and Keeler gap edges are modeled and conclusions regarding the shear and bulk viscosities of the ring are discussed. Finally, we explore the question of whether the radial density profile of the central and outer A ring, recently measured by Tiscareno and Harris (2018) in the highest resolution to date, and in particular, the sharp outer A ring edge can be modeled consistently from the balance of gravitational scattering by several outer moons and the mass and momentum transport. To this aim, the developed model is extended to account for the inward drifts caused by multiple discrete and overlapping resonances with multiple outer satellites and is then used to hydrodynamically simulate the normalized surface mass density profile of the A ring. This section of the thesis is based on studies by Tajeddine et al. (2017a) who recently discussed the common misconception that the 7:6 resonance with Janus alone maintains the outer A ring edge, showing that the combined effort of several resonances with several outer moons is required to confine the A ring as observed by the Cassini spacecraft. KW - celestial mechanics KW - diffusion KW - hydrodynamics KW - planets and satellites: rings KW - scattering Y1 - 2020 ER - TY - JOUR A1 - Sandin, C. A1 - Steffen, M. A1 - Schoenberner, D. A1 - Rühling, Ute T1 - Hot bubbles of planetary nebulae with hydrogen-deficient winds I. Heat conduction in a chemically stratified plasma JF - Frontiers in psychology N2 - 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. KW - conduction KW - hydrodynamics KW - planetary nebulae: general KW - stars: AGB and post-AGB KW - stars: Wolf-Rayet KW - X-rays: stars Y1 - 2016 U6 - https://doi.org/10.1051/0004-6361/201527357 SN - 1432-0746 VL - 586 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Thomas, Timon A1 - Feldmeier, Achim T1 - Radiative waves in stellar winds with line scattering JF - Monthly notices of the Royal Astronomical Society N2 - 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. KW - hydrodynamics KW - radiative transfer KW - waves KW - stars: winds KW - outflows Y1 - 2016 U6 - https://doi.org/10.1093/mnras/stw1008 SN - 0035-8711 SN - 1365-2966 VL - 460 SP - 1923 EP - 1933 PB - Oxford Univ. Press CY - Oxford ER - TY - GEN A1 - Sandin, Christer A1 - Steffen, Matthias A1 - Jacob, Ralf A1 - Schönberner, Detlef A1 - Rühling, Ute A1 - Hamann, Wolf-Rainer A1 - Todt, Helge Tobias T1 - The role of heat conduction to the formation of [WC]-type planetary nebulae T2 - Proceedings of the International Astronomical Union N2 - 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 582 KW - conduction KW - hydrodynamics KW - planetary nebulae: general Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-413702 SN - 1866-8372 IS - 582 SP - 494 EP - 495 ER - TY - GEN A1 - Seiß, Martin A1 - Spahn, Frank T1 - Hydrodynamics of Saturn’s dense rings T2 - Postprints der Universität Potsdam : Postprint Mathematisch Naturwissenschaftliche Reihe N2 - 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 574 KW - granular gas KW - instabilities KW - hydrodynamics KW - planetary rings Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-413139 SP - 191 EP - 218 ER - TY - JOUR A1 - Kurfuerst, P. A1 - Feldmeier, Achim A1 - Krticka, Jiri T1 - Time-dependent modeling of extended thin decretion disks of critically rotating stars JF - Astronomy and astrophysics : an international weekly journal N2 - Context. During their evolution massive stars can reach the phase of critical rotation when a further increase in rotational speed is no longer possible. Direct centrifugal ejection from a critically or near-critically rotating surface forms a gaseous equatorial decretion disk. Anomalous viscosity provides the efficient mechanism for transporting the angular momentum outwards. The outer part of the disk can extend up to a very large distance from the parent star. Aims. We study the evolution of density, radial and azimuthal velocity, and angular momentum loss rate of equatorial decretion disks out to very distant regions. We investigate how the physical characteristics of the disk depend on the distribution of temperature and viscosity. Methods. We calculated stationary models using the Newton-Raphson method. For time-dependent hydrodynamic modeling we developed the numerical code based on an explicit finite difference scheme on an Eulerian grid including full Navier-Stokes shear viscosity. Results. The sonic point distance and the maximum angular momentum loss rate strongly depend on the temperature profile and are almost independent of viscosity. The rotational velocity at large radii rapidly drops accordingly to temperature and viscosity distribution. The total amount of disk mass and the disk angular momentum increase with decreasing temperature and viscosity. Conclusions. The time-dependent one-dimensional models basically confirm the results obtained in the stationary models as well as the assumptions of the analytical approximations. Including full Navier-Stokes viscosity we systematically avoid the rotational velocity sign change at large radii. The unphysical drop of the rotational velocity and angular momentum loss at large radii (present in some models) can be avoided in the models with decreasing temperature and viscosity. KW - stars: mass-loss KW - stars: evolution KW - stars: rotation KW - hydrodynamics Y1 - 2014 U6 - https://doi.org/10.1051/0004-6361/201424272 SN - 0004-6361 SN - 1432-0746 VL - 569 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Seiss, Martin A1 - Spahn, Frank T1 - Hydrodynamics of saturn's dense rings JF - Mathematical modelling of natural phenomena N2 - 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. KW - granular gas KW - instabilities KW - hydrodynamics KW - planetary rings Y1 - 2011 U6 - https://doi.org/10.1051/mmnp/20116409 SN - 0973-5348 SN - 1760-6101 VL - 6 IS - 4 SP - 191 EP - 218 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Parkin, E. R. A1 - Broos, Patrick S. A1 - Townsley, L. K. A1 - Pittard, J. M. A1 - Moffat, Anthony F. J. A1 - Naze, Y. A1 - Rauw, G. A1 - Oskinova, Lida A1 - Waldron, W. L. T1 - X-RAY EMISSION FROM THE DOUBLE-BINARY OB-STAR SYSTEM QZ CAR (HD 93206) JF - ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES N2 - 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. KW - hydrodynamics KW - stars: early-type KW - stars: individual (QZ Carinae) KW - stars: massive KW - stars: winds, outflows KW - X-rays: stars Y1 - 2011 U6 - https://doi.org/10.1088/0067-0049/194/1/8 SN - 0067-0049 VL - 194 IS - 1 PB - IOP PUBLISHING LTD CY - BRISTOL ER - TY - THES A1 - Klar, Jochen T1 - A detailed view of filaments and sheets of the warm-hot intergalactic medium T1 - Eine detaillierte Ansicht der Filamente und Ebenen des warm-heißen intergalaktischen Mediums N2 - In the context of cosmological structure formation sheets, filaments and eventually halos form due to gravitational instabilities. It is noteworthy, that at all times, the majority of the baryons in the universe does not reside in the dense halos but in the filaments and the sheets of the intergalactic medium. While at higher redshifts of z > 2, these baryons can be detected via the absorption of light (originating from more distant sources) by neutral hydrogen at temperatures of T ~ 10^4 K (the Lyman-alpha forest), at lower redshifts only about 20 % can be found in this state. The remain (about 50 to 70 % of the total baryons mass) is unaccounted for by observational means. Numerical simulations predict that these missing baryons could reside in the filaments and sheets of the cosmic web at high temperatures of T = 10^4.5 - 10^7 K, but only at low to intermediate densities, and constitutes the warm-hot intergalactic medium (WHIM). The high temperatures of the WHIM are caused by the formation of shocks and the subsequent shock-heating of the gas. This results in a high degree of ionization and renders the reliable detection of the WHIM a challenging task. Recent high-resolution hydrodynamical simulations indicate that, at redshifts of z ~ 2, filaments are able to provide very massive galaxies with a significant amount of cool gas at temperatures of T ~ 10^4 K. This could have an important impact on the star-formation in those galaxies. It is therefore of principle importance to investigate the particular hydro- and thermodynamical conditions of these large filament structures. Density and temperature profiles, and velocity fields, are expected to leave their special imprint on spectroscopic observations. A potential multiphase structure may act as tracer in observational studies of the WHIM. In the context of cold streams, it is important to explore the processes, which regulate the amount of gas transported by the streams. This includes the time evolution of filaments, as well as possible quenching mechanisms. In this context, the halo mass range in which cold stream accretion occurs is of particular interest. In order to address these questions, we perform particular hydrodynamical simulations of very high resolution, and investigate the formation and evolution of prototype structures representing the typical filaments and sheets of the WHIM. We start with a comprehensive study of the one-dimensional collapse of a sinusoidal density perturbation (pancake formation) and examine the influence of radiative cooling, heating due to an UV background, thermal conduction, and the effect of small-scale perturbations given by the cosmological power spectrum. We use a set of simulations, parametrized by the wave length of the initial perturbation L. For L ~ 2 Mpc/h the collapse leads to shock-confined structures. As a result of radiative cooling and of heating due to an UV background, a relatively cold and dense core forms. With increasing L the core becomes denser and more concentrated. Thermal conduction enhances this trend and may lead to an evaporation of the core at very large L ~ 30 Mpc/h. When extending our simulations into three dimensions, instead of a pancake structure, we obtain a configuration consisting of well-defined sheets, filaments, and a gaseous halo. For L > 4 Mpc/h filaments form, which are fully confined by an accretion shock. As with the one-dimensional pancakes, they exhibit an isothermal core. Thus, our results confirm a multiphase structure, which may generate particular spectral tracers. We find that, after its formation, the core becomes shielded against further infall of gas onto the filament, and its mass content decreases with time. In the vicinity of the halo, the filament's core can be attributed to the cold streams found in other studies. We show, that the basic structure of these cold streams exists from the very beginning of the collapse process. Further on, the cross section of the streams is constricted by the outwards moving accretion shock of the halo. Thermal conduction leads to a complete evaporation of the cold stream for L > 6 Mpc/h. This corresponds to halos with a total mass higher than M_halo = 10^13 M_sun, and predicts that in more massive halos star-formation can not be sustained by cold streams. Far away from the gaseous halo, the temperature gradients in the filament are not sufficiently strong for thermal conduction to be effective. N2 - Im Rahmen der kosmologischen Strukturbildung entstehen durch Gravitationsinstabilitäten Flächen, Filamente und schließlich Halos. Interessanterweise befinden sich zu jedem Zeitpunkt der kosmologischen Entwicklung der Großteil der Baryonen nicht in den Halos, sondern in den Filamenten und Ebenen des intergalaktischen Mediums. Während diese Baryonen bei höheren Rotverschiebungen (z ~ 2) noch in Form durch die Absorbtion von Licht (von weit entfernteren Quellen) durch neutralen Wasserstoff bei einer Temperatur von T ~ 10^4 K beobachtbar sind (Lyman-Alpha Wald), gilt dies bei niedrigeren Rotverschiebungen für nur noch ca. 20 % der Baryonen. Der überwiegende Teil (ca. 50-70 % der gesamten baryonischen Masse) sind bisher noch nicht direkt beobachtbar. Numerische Simulationen sagen jedoch voraus, das sich diese Baryonen in den Filamenten und Flächen des kosmischen Netzes befinden. Die entsprechende Gasverteilung zeichnet sich durch hohe Temperaturen T = 10^5 - 10^7 K und geringe bis mittlere Dichten aus und wird als warm-heißes intergalaktisches Medium (WHIM) bezeichnet. Die hohen Temperaturen entstehen in Folge der Bildung von Stoßwellen und der darauf folgenden Erhitzung des Gases (shock-heating). Das WHIM ist daher hochgradig ionisiert und sein verlässlicher Nachweis stellt eine große Herausforderung für die beobachtende Kosmologie dar. Neuere hydrodynamische Simulationen zeigen, dass sich bei höheren Rotverschiebungen von z ~ 2 Gasströmungen entlang der Filamente bilden, die massive Galaxien mit erheblichen Mengen an relativ kaltem Gas (T ~ 10^4 K) versorgen können. Dies hätte einen erheblichen Einfluss auf die Sternentstehung in diesen Galaxien. Es ist daher von grundsätzlichem Interesse, die spezifischen hydro- und thermodynamischen Bedingungen in den Strukturen des WHIM zu untersuchen. Sowohl Dichte- und Temperaturprofile als auch Geschwindigkeitsfelder prägen spektroskopische Beobachtungen. Eine mögliche Mehrphasenstruktur des WHIM könnte daher als Indikator in beobachtenden Studien dienen. Im Zusammenhang mit den kalten Strömen ist es besonders interessant, Prozesse zu untersuchen die den Zufluss von kaltem Gas zu den Galaxien regulieren. Dies umfasst die Zeitentwicklung des Anteils an kaltem Gas in den Filamenten, sowie mögliche Mechanismen, die zum Versiegen des Zuflusses von kaltem Gas auf die Galaxienscheibe führen. Um diese Zusammenhänge zu erforschen, führen wir spezielle hydrodynamische Simulationen mit sehr hoher Auflösung durch, die zu ausgewählten, wohldefinierten Strukturen führen, die das WHIM charakterisieren. Wir beginnen mit einer ausführlichen Untersuchung des eindimensionalen Kollaps einer sinusförmigen Störung (pancake formation). Hierbei untersuchen wir den Einfluss von Strahlungkühlung, Heizung durch den intergalaktischen UV Hintergrund, Wärmeleitung, sowie von kleinskaligen Störungen, welche dem kosmologischen Störungsspektrum folgen. Wir benutzen hierbei eine Reihe von Simulationen, welche die Längenskala der anfänglichen Störung L als Parameter verwenden. Für L ~ 2 Mpc/h führt der Kollaps zur Ausbildung einer Stoßwelle. Zusätzlich entsteht als Folge der Strahlungskühlung und der Heizung durch den UV Hintergrund ein relativ dichter und kalter isothermer Kern. Mit ansteigendem L wird dieser Kern dichter und kompakter. Durch Wärmeleitung reduziert sich die räumliche Ausdehnung des Kerns. Für L ~ 30 Mpc/h führt dies zu einem Verschwinden des Kerns. Mit der Erweiterung unserer Methodik auf dreidimensionale Simulationen, entsteht nun eine Konfiguration, welche aus wohldefinierten Flächen, Filamenten und einem gasförmigen Halo besteht. Für L > 4 Mpc/h, erhalten wir Filamente, die vollständig durch Akkretionsschocks begrenzt sind. Wie in unseren eindimensionalen Simulationen weisen auch sie einen isothermen Kern auf. Dies legt nahe, dass das WHIM eine Mehrphasenstruktur besitzt und mögliche Spektralsignaturen erzeugen kann. Nach seiner Entstehung ist der Kern gegen weiteren Zufluss von Gas abgeschirmt und seine Masse reduziert sich mit der Zeit. In der direkten Umgebung des Halos entspricht der Kern des Filamentes den oben angesprochenen kalten Strömen. Unsere Untersuchung zeigt, dass diese während der gesamten Entwicklung des Halos existent sind. In der weiteren Entwicklung werden sie durch den expandierenden Akkretionsschock des Halos verengt. Ab einer Skala von L > 6 Mpc/h kann Wärmeleitung zu einem Verschwinden des Zustroms von kaltem Gas führen. Diese Skala entspricht Halos mit einer Gesamtmasse von M_halo = 10^13 M_sun. Galaxien, die sich in noch massiveren Halos bilden, können daher nicht durch kalte Ströme mit Gas für die Sternentstehung versorgt werden. Im Filament, weit außerhalb des gasförmigen Halos, sind die Temperaturgradienten zu klein, um effiziente Wärmeleitung zu ermöglichen. KW - Kosmologie KW - Hydrodynamik KW - Intergalaktisches Medium KW - cosmology KW - hydrodynamics KW - intergalactic medium Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-58038 ER - TY - THES A1 - Feldmeier, Achim T1 - Hydrodynamics of astrophysical winds driven by scattering in spectral lines N2 - Liniengetriebene Winde werden durch Impulsübertrag von Photonen auf ein Plasma bei Absorption oder Streuung in zahlreichen Spektrallinien beschleunigt. Dieser Prozess ist besonders effizient für ultraviolette Strahlung und Plasmatemperaturen zwischen 10^4 K und 10^5 K. Zu den astronomischen Objekten mit liniengetriebenen Winden gehören Sterne der Spektraltypen O, B und A, Wolf-Rayet-Sterne sowie Akkretionsscheiben verschiedenster Größenordnung, von Scheiben um junge Sterne und in kataklysmischen Veränderlichen bis zu Quasarscheiben. Es ist bislang nicht möglich, das vollständige Windproblem numerisch zu lösen, also die Hydrodynamik, den Strahlungstransport und das statistische Gleichgewicht dieser Strömungen gleichzeitig zu behandeln. Die Betonung liegt in dieser Arbeit auf der Windhydrodynamik, mit starken Vereinfachungen in den beiden anderen Gebieten. Wegen persönlicher Beteiligung betrachte ich drei Themen im Detail. 1. Windinstabilität durch Dopplerde-shadowing des Gases. Die Instabilität bewirkt, dass Windgas in dichte Schalen komprimiert wird, die von starken Stoßfronten begrenzt sind. Schnelle Wolken entstehen im Raum zwischen den Schalen und stoßen mit diesen zusammen. Dies erzeugt Röntgenflashes, die die beobachtete Röntgenstrahlung heißer Sterne erklären können. 2. Wind runway durch radiative Wellen. Der runaway zeigt, warum beobachtete liniengetriebene Winde schnelle, kritische Lösungen anstelle von Brisenlösungen (oder shallow solutions) annehmen. Unter bestimmten Bedingungen stabilisiert der Wind sich auf masseüberladenen Lösungen, mit einem breiten, abbremsenden Bereich und Knicken im Geschwindigkeitsfeld. 3. Magnetische Winde von Akkretionsscheiben um Sterne oder in aktiven Galaxienzentren. Die Linienbeschleunigung wird hier durch die Zentrifugalkraft entlang korotierender poloidaler Magnetfelder und die Lorentzkraft aufgrund von Gradienten im toroidalen Feld unterstützt. Ein Wirbelblatt, das am inneren Scheibenrand beginnt, kann zu stark erhöhten Massenverlustraten führen. N2 - Line driven winds are accelerated by the momentum transfer from photons to a plasma, by absorption and scattering in numerous spectral lines. Line driving is most efficient for ultraviolet radiation, and at plasma temperatures from 10^4 K to 10^5 K. Astronomical objects which show line driven winds include stars of spectral type O, B, and A, Wolf-Rayet stars, and accretion disks over a wide range of scales, from disks in young stellar objects and cataclysmic variables to quasar disks. It is not yet possible to solve the full wind problem numerically, and treat the combined hydrodynamics, radiative transfer, and statistical equilibrium of these flows. The emphasis in the present writing is on wind hydrodynamics, with severe simplifications in the other two areas. I consider three topics in some detail, for reasons of personal involvement. 1. Wind instability, as caused by Doppler de-shadowing of gas parcels. The instability causes the wind gas to be compressed into dense shells enclosed by strong shocks. Fast clouds occur in the space between shells, and collide with the latter. This leads to X-ray flashes which may explain the observed X-ray emission from hot stars. 2. Wind runaway, as caused by a new type of radiative waves. The runaway may explain why observed line driven winds adopt fast, critical solutions instead of shallow (or breeze) solutions. Under certain conditions the wind settles on overloaded solutions, which show a broad deceleration region and kinks in their velocity law. 3. Magnetized winds, as launched from accretion disks around stars or in active galactic nuclei. Line driving is assisted by centrifugal forces along co-rotating poloidal magnetic field lines, and by Lorentz forces due to toroidal field gradients. A vortex sheet starting at the inner disk rim can lead to highly enhanced mass loss rates. KW - Hydrodynamik KW - Strahlungstransport KW - Sternwinde KW - Akkretionsscheiben KW - hydrodynamics KW - radiative transfer KW - stellar winds KW - accretion disks Y1 - 2001 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-0000388 ER -