@phdthesis{Klar2012, author = {Klar, Jochen}, title = {A detailed view of filaments and sheets of the warm-hot intergalactic medium}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-58038}, school = {Universit{\"a}t Potsdam}, year = {2012}, abstract = {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.}, language = {en} } @article{ThomasPfrommerPakmor2021, author = {Thomas, Timon and Pfrommer, Christoph and Pakmor, R{\"u}diger}, title = {A finite volume method for two-moment cosmic ray hydrodynamics on a moving mesh}, series = {Monthly notices of the Royal Astronomical Society}, volume = {503}, journal = {Monthly notices of the Royal Astronomical Society}, number = {2}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stab397}, pages = {2242 -- 2264}, year = {2021}, abstract = {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.}, language = {en} } @article{ThomasPfrommer2019, author = {Thomas, T. and Pfrommer, Christoph}, title = {Cosmic-ray hydrodynamics}, series = {Monthly notices of the Royal Astronomical Society}, volume = {485}, journal = {Monthly notices of the Royal Astronomical Society}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stz263}, pages = {2977 -- 3008}, year = {2019}, abstract = {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.}, language = {en} } @phdthesis{Thomas2022, author = {Thomas, Timon}, title = {Cosmic-ray hydrodynamics: theory, numerics, applications}, doi = {10.25932/publishup-56384}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-563843}, school = {Universit{\"a}t Potsdam}, pages = {334}, year = {2022}, abstract = {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{\´e}n waves. The small-scale magnetic fields associated with these Alfv{\´e}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.}, language = {en} } @article{SandinSteffenSchoenberneretal.2016, author = {Sandin, C. and Steffen, M. and Schoenberner, D. and R{\"u}hling, Ute}, title = {Hot bubbles of planetary nebulae with hydrogen-deficient winds I. Heat conduction in a chemically stratified plasma}, series = {Frontiers in psychology}, volume = {586}, journal = {Frontiers in psychology}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201527357}, pages = {11}, year = {2016}, abstract = {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.}, language = {en} } @article{SeissAlbersSremčevićetal.2019, author = {Seiß, Martin and Albers, Nicole and Sremčević, Miodrag and Schmidt, J{\"u}rgen and Salo, Heikki and Seiler, Michael and Hoffmann, Holger and Spahn, Frank}, title = {Hydrodynamic Simulations of Moonlet-induced Propellers in Saturn's Rings}, series = {The astronomical journal}, volume = {157}, journal = {The astronomical journal}, number = {1}, publisher = {IOP Publishing Ltd.}, address = {Bristol}, issn = {0004-6256}, doi = {10.3847/1538-3881/aaed44}, pages = {11}, year = {2019}, abstract = {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{\´e}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{\´e}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.}, language = {en} } @phdthesis{Feldmeier2001, author = {Feldmeier, Achim}, title = {Hydrodynamics of astrophysical winds driven by scattering in spectral lines}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-0000388}, school = {Universit{\"a}t Potsdam}, year = {2001}, abstract = {Liniengetriebene Winde werden durch Impuls{\"u}bertrag von Photonen auf ein Plasma bei Absorption oder Streuung in zahlreichen Spektrallinien beschleunigt. Dieser Prozess ist besonders effizient f{\"u}r ultraviolette Strahlung und Plasmatemperaturen zwischen 10^4 K und 10^5 K. Zu den astronomischen Objekten mit liniengetriebenen Winden geh{\"o}ren Sterne der Spektraltypen O, B und A, Wolf-Rayet-Sterne sowie Akkretionsscheiben verschiedenster Gr{\"o}ßenordnung, von Scheiben um junge Sterne und in kataklysmischen Ver{\"a}nderlichen bis zu Quasarscheiben. Es ist bislang nicht m{\"o}glich, das vollst{\"a}ndige Windproblem numerisch zu l{\"o}sen, also die Hydrodynamik, den Strahlungstransport und das statistische Gleichgewicht dieser Str{\"o}mungen gleichzeitig zu behandeln. Die Betonung liegt in dieser Arbeit auf der Windhydrodynamik, mit starken Vereinfachungen in den beiden anderen Gebieten. Wegen pers{\"o}nlicher Beteiligung betrachte ich drei Themen im Detail. 1. Windinstabilit{\"a}t durch Dopplerde-shadowing des Gases. Die Instabilit{\"a}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{\"o}ntgenflashes, die die beobachtete R{\"o}ntgenstrahlung heißer Sterne erkl{\"a}ren k{\"o}nnen. 2. Wind runway durch radiative Wellen. Der runaway zeigt, warum beobachtete liniengetriebene Winde schnelle, kritische L{\"o}sungen anstelle von Brisenl{\"o}sungen (oder shallow solutions) annehmen. Unter bestimmten Bedingungen stabilisiert der Wind sich auf masse{\"u}berladenen L{\"o}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{\"u}tzt. Ein Wirbelblatt, das am inneren Scheibenrand beginnt, kann zu stark erh{\"o}hten Massenverlustraten f{\"u}hren.}, language = {en} } @article{SeissSpahn2011, author = {Seiss, Martin and Spahn, Frank}, title = {Hydrodynamics of saturn's dense rings}, series = {Mathematical modelling of natural phenomena}, volume = {6}, journal = {Mathematical modelling of natural phenomena}, number = {4}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0973-5348}, doi = {10.1051/mmnp/20116409}, pages = {191 -- 218}, year = {2011}, abstract = {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.}, language = {en} } @misc{SeissSpahn2011, author = {Seiß, Martin and Spahn, Frank}, title = {Hydrodynamics of Saturn's dense rings}, series = {Postprints der Universit{\"a}t Potsdam : Postprint Mathematisch Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Postprint Mathematisch Naturwissenschaftliche Reihe}, doi = {10.25932/publishup-41313}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-413139}, pages = {191 -- 218}, year = {2011}, abstract = {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.}, language = {en} } @article{KrtickaFeldmeier2018, author = {Krticka, Jiri and Feldmeier, Achim}, title = {Light variations due to the line-driven wind instability and wind blanketing in O stars}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {617}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201731614}, pages = {7}, year = {2018}, abstract = {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.}, language = {en} }